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Bcthodology  to  a  aacrotconoattric  aodal  to  facilitata  battar  analyiii 
Nith  tha  aodal.  First  and  second  dagraa  polynoaial  rasponsa  surfaca 
equations  express  endogenous  variables  as  functions  of  selected  exoge¬ 
nous  variables  in  the  Klein-Goldberger  econoaetric  aodal. 

Second  degree  polynoaial  response  functions  are  used  to  derive 
function  aultipliers.  The  function  aultipliers  shOM  that  the  iapact  of 
changes  in  exogenous  variables  on  endogenous  variables  depends  on  the 
levels  of  one  or  aore  other  exogenous  variables.  The  function  eulti- 
pliers  are  used  to  conduct  policy  analysis  and  assess  factor  iaportanca. 
As  an  extension,  first  degree  polynoaial  response  functions  are  used  in 
an  example  problem  to  maximize  gross  national  product  subject  to  con¬ 
straints  on  unemployment,  inflation,  and  ranges  of  fiscal  policy 
variables.  The  example  problem  demonstrates  the  flexibility  and  value  of 
developing  a  response  surface  equation  for  complex  macroeconoaetric 
models.  _ _ 


The  study  concludes  that  a  response  surface  can  capture  the 
complexity  of  macroeconometr i c  models  such  as  the  Klein-Goldbergcr 
model.  Results  also  shOH  that  the  assumptions  of  linearity  for 
developing  multipliers  can  result  in  misleading  values  when  nonlinearity 
IS  present.  Recommendations  for  further  research  include  fitting  a  more 
nonlinear  model  with  response  surfaces,  and  including  time  as  an 
independent  variable  in  the  response  surface  equations. 
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AN  APPLICATION  OF  RESPONSE  SURFACE  NETH0D0L06Y 


TO  A  NACROECONOHIC  MODEL 


!•  Ibf  BflSSCSb  0C5bifi! 


Introduction 

This  thesis  dtsls  with  hoM  ascroiconoaic  aodils  and  response 
surface  nethodoloqy  (RSM)  can  be  brought  together  to  provide  better 
analysis  of  a  national  econoey.  Its  purpose  is  to  deeonstrate  that  RSM 
can  reduce  coeplex  relationships  enbodied  in  eacroeconoeetric  eodels  to 
simple  equations.  The  thesis  also  interprets  the  simple  equations  and 
shows  how  they  can  be  used  for  practical  applications.  This  chapter 
puts  the  research  effort  into  perspective  by  briefly  describing  macro- 
econometric  models  and  response  surface  methodology  suggesting  possible 
wavs  to  combine  the  two.  The  chapter  then  outlines  the  research  plan 
including  the  research  problem,  research  questions,  research  objectives, 
scooe,  and  general  methodology  for  attacking  the  problem. 
b§£i:5Scgngmic  Mgdfls 

Macroeconomic  models  are  a  set  of  economic  relationships  expressed 
in  mathematical  equations  which  allow  economists  to  predict  the  perfor¬ 
mance  of  a  national  economy.  Economists  have  developed  several  types  of 
macroeconomic  models.  One  type  uses  certain  economic  indicators  which 
have  historically  led  cyclical  changes  in  the  economy.  Another  type 
uses  consumer  attitudes  and  buying  plans  to  predict  economic  perfor¬ 
mance.  The  type  of  particular  interest  to  this  thesis  is  the  econome¬ 
tric  model.  Econometric  models  are  systems  of  statistically  derived 
simultaneous  equations  based  on  theory  and  historical  data  which  predict 


AF1T/60R/ENS/850-4 


AN  APPLICATION  OF  RESPONSE  SURFACE  NETHODOLOSY 
TO  A  NACROECONOMIC  MODEL 

THESIS 

Priaintid  to  tht  Faculty  to  tho  School  of  Enginitring 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fullfillnent  of  the  Requireeents  for  the  Degree  of 
Master  of  Science  in  Operations  Research 


Janes  L.  Donovan,  B.  S. 
Captain,  USAF 

DecenOer  1985 


Approved  for  public  release:  distribution  unlimited 


r»i.Tf  irTW”3r".rFT-',»-  r'.'^ 


AcknoM^edggAents 

First,  I  Mish  to  thank  God,  Mho  mss  always  there  to  help  during  the 
very  rare  occasions  when  ny  thesis  advisors  were  unavailable.  I  sin¬ 
cerely  thank  ny  two  thesis  advisors  LtCol  Palner  W.  Smith  and  Dr  Robert 
F,  Allen  for  their  many  hours  of  help  and  constant  encouragement.  I  have 
learned  much  from  their  technical  expertise  and  personal  examples. 

I  sincerely  appreciate  the  patience  and  support  of  my  family  and 
especially  Michele,  who  is  soon  to  become  my  family.  This  thesis  is  a 
oroduct  of  their  sacrifice  as  well  as  mine. 

Finally,  I  would  like  to  thank  Mrs.  Joanne  Crane,  who  graciously 
volunteered  many  long  hours  to  help  me  prepare  this  document.  Without 
her  assistance  this  thesis  would  have  been  nearly  impossible  to  complete 
on  time. 


1  Accesion  For  ^ 

NTIS  CRA&I 
DllC  TAB 

□ 

U.. announced 

□ 

Jj'jtit, cation 

By  . 

1  Diot  ibution/ 

. 

Availability  Codes 

t  1  Avail 

and  1  or 

VI.  Rtsponsi  Surfaca  Analysis  . . . . .  65 


Introduction  .  65 

Rasponsa  Surfaca  Intarpratation  .  65 

Coaputing  Multipliars  . . . .  69 

Policy  Siaulation  . . .  70 

Factor  laportanca  .  72 

Optiaization  Applications  . . . . .  74 

VII.  Concluding  Raaarks  .  87 

Rasaarch  Suaaary  .  87 

Racoaaended  Furthar  Rasaarch  .  90 

Appandix  A:  Solving  the  Klein-Goldberger  Modal  .  92 

Appendix  £t:  The  Equivalence  of  Multipliers  and  Least  Squares 

CoeHicients  for  Linear  Systeas..... .  109 

Appendix  C:  Decoding  Second  Order  Coa'f'f iciants  .  110 

Appendix  D;  Decoding  the  Experiaental  Design  .  115 

Appendix  E:  Response  Surface  Coefficients  .  121 

Appendix  F:  Qptinization  Problen  Solution  .  130 

Bibliography  .  132 


List  of  Figurfs 


Exanple  Response  Surface  . 

The  Coaaodities  Market  . 

Incone  Effects  of  Increased  Sovernnent  Spending  . 

The  Evolving  Econony  . 

Dependence  of  dQ/dTc  on  6  in  Five  Years  . 

Relationship  Between  Q  and  Tc  at  Different 
Levels  of  Q  in  Period  Five  . 

A  Graphical  Conparison  of  Multipliers  for  Changes 
in  Q  Due  to  Changes  in  G  in  Period  Five . 


List  of  Tablts 


The  Kltin-Soldbtrgtr  nodtl  . 

Glossary  of  Variablss  for  the  Klein-Goldberger  Model  .... 

The  Klein  Model  I  Linear  Econonetric  Model  . 

Three  Factor  Three  Level  Factorial  Experinental  Design  .. 

Sunnary  Table  of  Stepwise  Regression  Results  for  Nunber 
of  Workers  Enployed  in  Period  Zero  . 

Sunnary  Table  of  Stepwise  Regression  Results  for  the 
Price  Index  in  Period  Zero  . . . . . 

Sunnary  Table  of  Stepwise  Regression  Results  for  Gross 
National  Product  in  Period  Zero  . . . 

Sunnary  Table  of  Stepwise  Regression  Results  for  Nunber 
of  Workers  Enployed  in  Period  Five  . 

Sunnary  Table  of  Stepwise  Regression  Results  for  the 
Price  Index  in  Period  Five  . 

Sunnary  Table  of  Stepwise  Regression  Results  for  Gross 
National  Product  in  Period  Five  . 

Design  Point  Fit  Check  for  Second  Order  Response  Surface. 

Randon  Point  Fit  Check  for  Second  Order  Response  Surface. 

R‘  Values  for  the  First  Order  Response  Surface 
Equation  .  . . . 

Multipliers  for  Unit  Increase  in  Governnent  Spending  .... 

Response  Surface  Coefficients  for  a  Unit  Increase  in 
Governnent  Spending  . . . 

Multiplier  Comparison  for  Changes  in  Q  Due  to  Changes 
G  i n  Per i od  Five  . . . 

Optimal  Fiscal  Policy  for  the  Example  Problem  . 

K1 ei n-Gol dberger  Model  Solutions  in  the  Area  of 
Alleged  Optimality  . 


Data  for  Selecting  the  Best  Variable  to  Alter 
Shadow  Prices  for  Fiscal  Policy  Problem  . 


✓  ✓  /y 


economic  pirf ornenci. 


To  formulate  an  econometric  model,  economiets  first  hypothesize 
equations  to  describe  theoretical  economic  relationships  far  sectors  of 
the  economy.  For  example,  in  one  model,  manufacturing  output  is  a 
function  of  the  amount  of  hours  worked  in  the  manufacturing  sector,  the 
amount  of  money  invested  in  manufacturing,  and  average  productivity. 
Honey  invested  in  manufacturing  is  a  function  of  manufacturing  capacity 
used  in  the  previous  period,  manufacturing  output  in  previous  periods, 
cash  flow  in  the  manufacturing  sector,  and  interest  rates  on  bonds.  The 
number  of  hours  worked  is  a  function  of  previous  period  manufacturing 
output,  wage  rates,  and  percent  of  manufacturing  capacity  used.  Wage 
rates  depend  on  past  wage  rates  and  the  cost  of  living,  and  so  on 
(Evans,  1969:433-442).  Economists  postulate  the  form  of  functions  such 
as  these  and  then  use  historical  data  to  estimate  unknown  coefficients 
in  the  equations  with  statistical  techniques.  Finally,  they  put  equa¬ 
tions  representing  all  sectors  of  the  economy  together  and  solve  them 
simultaneously  to  obtain  predictions. 

Klein  and  Evans  ennumerate  three  major  uses  for  econometric  models 
in  their  work  entitled.  The  Wharton  Econometric  Forecasting  Hgdel  (Klein 
and  Evans,  1968:50).  First,  economists  can  use  them  for  prediction. 
Second,  econometric  models  can  simulate  the  consequences  of  economic 
policies  such  as  tax  increases,  government  spending,  and  Federal  Reserve 
actions  for  periods  in  the  past.  In  many  ways  the  most  importaht  use 
for  econometric  models,  claim  Klein  and  Evans,  is  computing  multipliers 
for  fiscal  and  monetary  policy  alternatives.  Fiscal  policy  is  manipula¬ 
ting  the  economy  by  government  spending  and  taxation,  Monetary  policy 


is  manipulating  tne  economy  through  actions  ot  the  Federal  Reserve  such 
as  selling  government  bonds,  changing  the  required  reserves  that  banks 
must  hold  against  demand  deposits,  etc.  A  multiplier  is  a  constant 
which,  when  multiplied  by  a  change  in  an  econometric  input  variable, 
gives  the  change  in  an  output  variable.  Multipliers  are  especially 
important  because  ot  the  increasing  size  and  complexity  ot  econometric 
models.  For  example,  the  Uarton  Econometric  and  Forecasting  Unit  has  S3 
equations  and  29  identities  (Evans,  1969:442).  A  model  designed  at 
Brookings  has  over  ISO  equations  (Evans,  1969:S03).  The  ettects  ot 
changes  in  certain  input  variables  are  ditticult  to  trace  through  to  the 
tinal  output  ot  such  large  models  unless  multipliers  are  computed. 
Unfortunately,  multiplier  analysis  does  not  account  well  tor  nonlinear 
systems  with  interactions  among  input  variables.  Response  surface 
methodology  might  provide  a  way  ot  overcoming  these  difficulties. 
5§sS9Dse  Surface  Mfthodglogy 

Response  Surface  Methodology  (RSM)  is  an  analytical  tool  tor  model¬ 
ling  a  very  complex  or  unknown  process  with  a  single  mathematical  equa¬ 
tion  and  exploring  the  resulting  relationship  between  the  inputs  and  an 
output  ot  the  process.  A  response  surface  results  from  plotting  output 
values  obtained  from  the  response  surface  equation  against  input 
variables  as  they  vary  over  continuous  ranges.  Chemical  reactions  are 
classic  RSM  applications.  Factors  which  affect  the  yield  of  a  chemical 
reation  are  temperature,  pressure,  amount  of  reactants,  and  reaction 
time.  If  one  were  to  set  the  amount  of  reactants  at  some  set  level  and 
fi:;  the  time  the  reactants  are  allowed  to  react,  a  set  of  experiments 
could  be  conducted  at  various  combinations  of  temperature  and  pressure. 
Bv  recording  the  vield  of  each  experiment  and  plotting  that  yield 


against  tenperature  and  prissura,  tha  rasult  mould  be  soecthing  like  the 
graph  of  the  response  surface  in  Figure  1.1. 


Figure  1.1.  Exaeple  Response  Surface. 


Box  and  Wilson  first  introduced  response  surface  nethodology  in 
their  1951  paper  entitled,  “On  the  Experieental  Attainnent  of  Optieue 
Conditions"  (Box  and  Wilson,  1951).  Since  then,  nany  researchers  have 
profitably  applied  the  technique  to  problens  in  chemistry,  foodstuffs, 
tool  life  testing,  and  other  areas  (Hill  and  Hunter,  1966:576).  Smith 
and  Mellichamp  first  demonstrated  that  RSH  could  provide  valuable  in¬ 
sight  into  complex  deterministic  analysis  models  (Smith  and  Hellichamp, 
1979).  Based  on  Smith  and  Mel  1 i champ  '  s  work,  students  at  the  Air  Force 
Institute  of  Technology  have  applied  RSn  to  several  deterministic  and 
probabilistic  models  (Manacapi 1 1 i ,  1964;  Qraney,  1984;  Heitzler,  1964; 
Sparrow,  1984).  Deterministic  models  are  mathematical  representations 
of  an  underlying  process  for  which  given  inputs  to  the  process  yield  the 
same  output  every  time  the  model  is  run.  The  outputs  of  probabilistic 


aodils  contain  randoa  variation.  Burdick  and  Naylor  fuggiitid  applying 
RSM  to  aconoRCtric  aodels.  Thay  aftoNad  how  to  coabine  a  siapla  six 
equation  econoaatric  aodal  and  a  utility  function  for  aconoaic  policy 
optimization  by  response  surface  techniques  (Burdick  and  Naylor, 
1969:29).  Houever,  they  did  not  actually  estiaate  response  surfaces  for 
this  system.  Their  ideas  aerit  a  more  thorough  development. 

Applying  RSN  to  a  problem  involves  several  steps.  First  the  pro¬ 
blem  is  defined  and  variables  of  interest  are  specified.  Next  a  response 
surface  equation,  usually  a  Ion  order  polynomial,  is  selected  to  aodel 
the  process  under  study.  Based  on  the  response  surface  equation 
selected,  an  appropriate  experimental  design  is  chosen  which  compro¬ 
mises  between  economy  of  design  points  and  orthogonality.  Then  the 
experiment  or  model  is  run  repeatedly  at  factor  levels  specified  by  the 
exoeriaental  design.  With  the  data  collected  from  the  experiment  or 
model  runs,  a  response  surface  equation  is  estimated  using  ordinary 
least  squares.  After  checking  for  adequate  fit,  the  response  surface  is 
ready  for  interpretation  and  analysis.  These  application  steps  are 
covered  in  much  greater  detail  in  Chapter  III. 

Several  analysis  methods  are  available  for  exploring  response  sur¬ 
faces.  nost  of  these  methods  are  devised  to  enable  the  analyst  to 
optimize  the  underlying  process  or  model  that  the  response  surface 
reoresents.  They  include  the  method  of  steepest  ascent,  classical 
optimization  using  calculus,  Lagrangian  techniques,  mathematical  pro¬ 
gramming,  and  others.  However,  optimization  is  not  the  onlv  use  for 
response  surfaces.  Examining  the  response  surface  equation  itself 
reveals  characteristics  about  the  underlying  process  or  model.  In 
addition,  the  impact  of  tradeoffs  between  inputs  is  easily  assessed. 


Thtfc  analysis  ttchniquis  hava  potintial  applications  for  studying 
nacrotconosttric  sodtls.  For  instancs,  constants  or  sxprissions  siailar 
to  aultiplisrs  could  be  estieatsd.  These  constants  or  expressions  would 
sake  the  inpact  of  changes  in  econoeic  policy  variables  explicit. 
Multipliers  ieplicitly  assuee  that  the  changes  in  response  are  at  least 
approxinately  linearly  related  to  the  change  in  input.  Expressions 
derived  to  serve  as  "nultipliers"  using  the  response  surface  technique 
have  no  such  inplicit  assueption.  Thus  alternative  policy  options  for 
nonlinear  nodels  could  be  evaluated  nore  accurately.  In  fact,  given 
specific  econoeic  objectives,  econoeic  policy  could  be  optinized.  The 
ootential  applications  for  RSM  in  eacroecononics  suggest  several  areas 
for  research.  Below  is  the  specific  olan  for  this  thesis  effort. 
Research  Plan 

^roglen  Econoeists  have  developed  large  econoeetric 

nodels  to  predict  the  perforeahce  of  national  econoeies.  Unfortunately, 
because  of  the  conplexity  of  these  nodels,  economists  have  difficulty 
investigating  the  effects  of  changing  key  input  variables  on  economic 
performance.  Response  surface  methodology  may  be  able  to  reduce  key 
relationships  in  the  model  to  a  single  equation. 

BS§S$!Ich  Question.  How  well  can  response  surface  methodology 
capture  the  predictive  power  of  a  large  econometric  model  and  can 
response  surface  methodology  simplify  sensitivity  analysis  of  such  a 
model  ? 

id^ll^iar^  Questions.  Several  issues  related  to  the  research 
deserve  investigation. 

1.  Can  a  response  surface  based  on  a  simple  function  accurately 
capture  the  relationships  in  a  large  macro-economic  model? 
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f'V  '. 

2.  Is  thtra  a  Unit  to  thi  siza  of  aodal  that  rasponsa 
surfaca  aathodology  can  handla? 

V'  * ' 

3.  Hon  can  tha  coafficiants  of  the  rasponsa  surface  V 

equation  be  interpreted?  >: 


i 


4.  Do  the  rasponsa  surface  coefficients  identify  aconoaic  inputs 
Nhich  are  aost  influential  in  driving  a  national  econoay? 

5.  Can  one  use  rasponsa  surface  aethodology  to  dataraine  fiscal 
and  aonatary  policy  Nhich  tha  aodal  predicts  Nill  optiaiza 
particular  aaasuras  of  aconoaic  parforaanca? 

BS§SiC£d  QbjSStiyas.  To  ansNer  the  research  questions,  several 
objectives  aust  be  aat.  They  are: 

1.  Deteraine  hON  Nell  a  response  surface  can  fit  the  response  of 
iaportant  aconoaic  variables  to  changes  in  fiscal  and  aonetary 
variables  for  an  actual  aacroeconoaetric  aodel.  The  aodel 
should  be  aoderatly  sized,  have  soae  nonl inaar ities,  and  have 
characteristics  Nhich  are  Nell  knoNO  froa  previous  analysis. 

2.  Verify  that  response  surface  does  in  fact  reflect  aodel 
characteristics  by  coaparing  response  surface  equation  para- 
aeters  to  aultipliers  coaputed  for  the  aodel. 

3.  Interpret  tha  response  surface  equations. 

4.  Develop  applications  of  practical  value  for  tha  response  sur¬ 
face  equations. 

i&SDS'  study  deaonstrates  feasibile  applications  for  response 

surface  aethodology  techniques  in  aacroeconoaics.  To  aake  the  research 
effort  aanageable,  several  decisions  are  aade.  The  study  uses  the 
Klein-Goldberger  aacroeconoaetric  aodel  for  the  investigation.  Reasons 
for  selecting  the  aodel  are  given  in  Chapter  II.  Based  on  previous 
experience  with  response  surface  aethodology,  a  second  order  polynoaial 
response  surface  equation  is  assumed.  The  model  is  constructed  to 
mirror  the  characteristics  of  the  national  econoay  which  changes  slowly 
in  response  to  changes  in  policy  and  so  a  second  degree  polynomial 
should  adequately  fit  the  model.  Horeover,  Goldberger  argues  that  the 
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Kltin-6oldbirgtr  nodel  is  nearly  linear  (Goldberger,  1959:136-138). 
Response  surfaces  are  built  for  three  output  variables  in  teres  of  five 
input  variables.  Although  response  surfaces  Mith  aore  variables  could 
be  constructed,  no  particular  advantage  is  seen  in  this.  Although  any 
econometric  model  is  stochastic  in  nature,  it  is  assumed  that  the 
model  is  deterministic.  Finally  it  is  assumed  that  the  model  is  a 
reasonably  valid  representation  of  the  economy's  behavior  with  the 
exception  of  deficiencies  in  the  monetary  sector  Mhich  will  be  discussed 
in  Chapter  II.  No  effort  is  made  to  evaluate  the  model's  forecasting 
record . 

QtQSCiL  tl&ttlSdal.gg^.  The  general  plan  of  attack  for  accomplishing 
the  research  objectives  is: 

1.  Develop  a  computer  program  which  solves  the  Klein-Soldberger 
model  for  output  variables  in  terms  of  given  input  variable 
values. 

2.  Solve  the  model  for  values  of  input  variables  reqiured  by  the 
experimental  design  selected. 

3.  Fit  a  second  order  polynomial  response  function  to  the  data  and 
check  fit. 

4.  Fit  a  first  order  polynomial  model  to  data  for  a  direct 
comparison  of  response  surface  coefficients  to  multipliers 
computed  by  Goldberger. 

5.  Interpret  response  surface  equations  and  develop  ways  for 
summarizing  the  information  contained  in  the  eouations. 

6.  Develop  practical  applications  for  derived  response 
surfaces.  Specifically,  develop  optimization  aopl ications. 

ItltSLS  The  chapters  in  this  thesis  follow  the  pattern  in 

this  chapter.  Chapter  II  examines  macroeconomic  models  in  general  and 
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the  Klein-Goldberger  model  in  particular.  It  also  describes  solution 
techniques  and  methods  for  deriving  multipliers.  Chapter  III  describes 
the  steps  in  applying  response  surface  methodology  and  discusses  how  RSN 


can  bi  appliid  to  aacroiconoaic  aodtli.  Chaptir  IV  ditailf  aithodology 
for  thif  rasaarch  effort.  Chaptar  V  addratias  hoM  wall  tha  ratponaa 
surfacai  fit  tha  aodal,  and  coaparat  first  ordar  rasponta  surfaca 
coafficianta  to  Goldbarg's  aultipliars.  Chaptar  VI  intarprats  faaturas 
of  tha  darivad  reaponaa  surfacas  and  davalops  an  optiaazation  problaa 
application  for  the  rasponsa  surfacas.  Finally,  Chaptar  VII  suaaarizes 


findings  and  recomnands  further  rasaarch. 
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Intrgdygtign 

Chaptir  I  daacribtd  what  aacrotconoaetric  aodala  art  and  hoa  thay 
ara  uaad.  Tha  chaptar  axplalnad  that  macroaconoaatrlc  aodala  ara  a  aat 
O'f  siaultaneous  aquations  basad  on  theory  and  historical  data  ahich 
alloNS  aconoaists  to  pradict  tha  perforaance  of  the  national  econoay. 
Uses  for  econoaetric  aodels  include  forecasting,  policy  siaulation  in 
historical  periods,  and  aost  laportantly,  for  coaputing  aultipliers 
which  relate  changes  in  fiscal  and  aonetary  policy  to  changes  in 
econoaic  perforaanca.  This  chaptar  explains  aspects  of  a  staple  aacro- 
econoaic  aodel ,  and  uses  this  aodel  to  analyze  how  changes  in  fiscal 
policy  affect  econoaic  perforaance.  Next  the  Klein*6oldberger  (KG) 
aodel  is  introduced.  After  discussing  aacroaconoai c  aodels,  this  chap¬ 
ter  discusses  aethods  for  solving  aacroaconoaetric  aodels  and  deriving 
aultipl iers. 

Describing  how  aacroeconoaic  aodels  are  built  is  often  the  subject 
of  an  entire  college  course.  The  discussion  of  aacroeconoaics  here  is 
aerely  aeant  to  be  a  quick,  siaplified  review  of  points  relevant  to  the 
research  effort.  The  aaterial  presented  is  condensed  froa  Baird  and 
Cassuto's  introductory  text  entitled  Hacroecgnga j.c5L  hgnetaryj.  Sgarchi 
iQcgae  Ihegries  (Baird  and  Cassuto,  19B4).  The  text  is  thorough  yet 
extreaely  readable  with  plenty  of  helpful  exaaples  and  illustrations. 

The  reader  who  is  unfaailiar  with  aacroeconoaics  is  highly  encouraged  to 
consult  the  Baird  and  Cassuto  text  or  a  siailar  text.  The  discussion 
below  begins  with  econoaic  equalibriua.  Then  the  national  incoae 
identity  is  used  to  develop  the  coaaodities  market  of  a  general 
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Tha  priaary  atauaptlon  bahind  aqualibriua  aacroaconoaic  aodala 
la  that  tha  aconoay  alwayf  aaakt  an  aqualibriua.  A  rigorouf  arguaant 
supporting  this  assartion's  truth  Nill  not  ba  attaaptad;  hoNavar,  tha 
proposition  aakas  intuitiva  sansa.  Natura  is  full  of  axaaplas  of  systaas 
which  saak  aqualibriua.  A  dislodgad  bouldar  rolls  down  tha  aountainsida 
until  is  finds  a  valley  to  rest  in.  Cheaicals  react  until  they  reach 
aqualibriua.  Econoaic  theory  assuaas  that  tha  aconoay  will  also  saak 
aqualibriua  in  tha  absence  of  external  disturbances.  With  this  assuap- 
tion  this  in  aind,  the  aacroaconoaic  aodal  building  discussion  aay 
begin. 

I!lS  (Hitional  iQcgae  Accgunti,ng  Ident^^y.  Perhaps  tha  aost  widely 
used  perforaance  aeasura  of  an  aconoay  is  tha  gross  national  product 
(GNP).  This  nuabar  is  a  aeasura  of  national  incoaa.  It  is  defined  as  tha 
dollar  value  of  all  final  goods  and  services  produced  for  final  consuap- 
tion  during  a  calendar  year,  natheaatically  GNP  can  ba  defined  as: 

Y  »  C  +  Ia  ♦  G  +  Fc  (2.1) 

where 

Y  s  gross  national  product, 

C  -  consuaption  of  goods  and  services, 

1a  s  actual  invastaant, 

G  a  govarnaent  expenditures,  and 
Fc  a  net  exports  to  foreign  nations. 

In  words,  the  equation  says  GNP  is  tha  total  of  a  country's 
expenditures  on  final  goods  and  services  plus  the  value  of  net  exports. 
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Doublt  entry  accounting  procedure!  for  producing  Hrei  require  that  the 
expenditures  far  producing  those  final  goods  and  services  be  equal  to 
the  incoee  froe  the  sales  of  those  goods  and  services,  hatheeatically , 

C  +  Ia  ♦  G  +  Ft  *  Ya  ♦  Sb  ♦  T  (2.2) 

Mhere 

Yd  =  disposable  incoee, 

Sb  3  business  savings,  and 

T  s  taxes. 

The  left  side  of  Eq  (2.2)  is  total  expenditures  and  the  right  side 
is  total  incoee.  Eq  (2.2)  is  knoMn  as  the  national  incoee  accounting 
identity  and  it  fores  the  basis  for  eacroeconoeic  eodel  building.  On  the 
right  side,  disposable  incoee,  Yd,  is  after  tax  after  business  saving 
incoee.  Business  savings,  Sb,  are  the  portion  of  business  net  incoee 
which  fires  do  not  distribute  to  owners.  For  this  discussion,  business 
savings  will  be  assueed  set  a  fixed  level,  Sb.  Finally  taxes,  T,  are 
incoee  appropriated  by  the  governeent.  For  this  sieple  eodel,  it  will 
be  assueed  that  taxes  are  fixed  at  T. 

On  the  left  side  of  Eq  (2.2)  are  the  coeoonents  of  total 
expenditure  including  consueption,  investeent,  governeent  spending,  and 
net  foreign  exports.  If  each  of  the  coeponents  of  GNP  can  be  detereined, 
GNP  can  be  coeputed.  The  coeeodities  market  models  the  relationship 
between  components  of  national  income. 

Jhe  Commodities  Market.  All  of  the  analysis  which  follows  assumes 
fixed  prices.  A  variable  price  model  requires  development  of  other 
markets  in  the  economy.  The  first  component  of  the  expenditures  side  of 


the  national  incofli  accounting  identity  is  consueption.  Consueption  is 
in  large  part  based  on  disposable  incone.  The  more  incone  a  person 
receives,  the  nore  that  person  usually  spends.  Hathenatical ly ,  this  is 
Mritten  C  -  C(Yd).  A  sinple  consueption  function  is  C  -  C  *  bVd  , 
where  C  is  autononous  consunption  that  occurs  without  incone,  b  is  the 
narginal  propensity  to  consume,  and  Yd  is  disposable  incone.  In  general 
b  IS  between  zero  and  one  because  individuals  divide  their  disposable 
incoflie  between  savings  and  consumption. 

The  next  component  of  expenditures  is  investment.  Planned  invest¬ 
ment  IS  also  a  function  of  income  because  the  more  sales  a  firm  receives 
the  more  it  will  want  to  expand  operations  through  capital  investment.  A 
simple  investment  function  is  I  =  I  >  vY  ,  where  I  is  the  autonomous 
part  of  investment  and  v  is  the  marginal  propensity  for  firms  to  invest. 
V  IS  between  zero  and  one  because  firms  expend  income  received  on 
profits,  operating  costs,  etc.,  as  well  as  investment. 

The  next  component  of  expenditures  is  government  spending.  Govern¬ 
ment  spending  is  set  by  the  government  instead  of  market  forces.  Govern¬ 
ment  expenditures  is  assumed  to  be  set  at  a  particular  level,  say 
G. 

The  last  component  of  expenditures  is  net  foreign  exports  which  is 
the  difference  between  exports  and  imports.  For  the  purposes  of  this 
analysis  net  foreign  exports  is  assumed  zero. 

Figure  2.1  shows  the  relationships  between  components  in  the 
national  income  accounting  identity  which  together  comprise  the  commodi¬ 
ties  market.  The  vertical  axis  measures  total  expenditures,  E,  and  the 
horizontal  axis  measures  total  income.  Y. 


Figure  2.1.  The  Connodities  Market 


The  45*^  line  is  the  national  income  accounting  identity  and  nutt 
alMays  hold,  Line  DO  represents  the  sum  of  consumption,  expected 
investment,  and  government  expenditures  known  as  aggregate  demand.  Even 
though  aggregate  demand  is  the  sum  of  consumption,  investment,  and 
government  spending,  it  is  different  from  GNP  because  it  includes 
planned  investment  by  firms  instead  of  actual  investment.  It  is  possible 
that  planned  investment  will  not  equal  actual  investment.  Included  in 
investment  are  inventories  of  goods  produced  for  sale.  If  demand  for 
goods  IS  lower  than  expected,  inventories  will  increase  in  the  short  run 
and  actual  investment  will  be  higher  than  planned.  The  difference 
between  expected  investment  and  actual  investment  is  the  unplanned 
increase  in  inventories.  Expected  investment  equals  actual  investment 
onlv  at  equilibrium  in  the  commodities  market  (ie.,  when  unplanned 
changes  in  inventories  are  zero  and  the  45*^  line  and  the  aggregate 


demand  line  Intersect).  In  Figure  2.1  Y*  is  the  equalibrium  income  and 
E*  is  the  equalibrium  total  spending.  The  system  oi  equations  which 
describes  economic  equalibrium  in  the  hypothetical  economy  is 


Y  =  C  +  Ia  ♦  G  (2.3) 

Y  »  Y.,  +  Sh  +  T  (2.4) 
C  =  C  +  bY  (2.5) 
I  =  I  +  vY  (2.6) 
G  =  6  (2.7) 
T  »  f  (2.8) 
Sh  =  Sa  (2.9) 


G,  T,  and  St.  are  assumed  Hxed.  In  the  system  of  equations  above, 

Y,  C,  I,  and  Ya  are  known  as  endogenous  variables.  Endogenous  variable 
are  variables  which  have  values  determined  within  the  system.  6,  T,  and 
St.  are  exogenous  variables.  Exogenous  variables  are  assumed  to  have  a 
value  determined  outside  the  system.  Although  the  system  above  is  a 
complete  system  with  the  number  of  unknowns  equal  to  the  number  of 
equations,  it  is  complete  only  because  of  the  simplifying  assumptions 
used  to  formulate  the  equations.  For  instance,  prices  are  assumed  to  be 
constant.  Since  Y  =  PQ  where  P  is  the  price  level  and  Q  is  the  real 
output  of  the  economy,  any  increases  in  Y  are  assumed  to  be  increases  in 
real  output,  Q.  Also  it  has  been  assumea  that  consumption  and  investment 
are  functions  only  of  income,  and  that  business  savings  are  constant. 

The  system  of  equations  may  be  solyed  for  any  of  the  endogenous 
variables  in  terms  of  constants  and  exogenous  variables.  Solving  the 
equations  in  this  way  gives  the  explicit  effect  of  exogenous  variables 


on  an  endogenous  variable.  Since  the  analysis  has  centered  on 
deterninants  ot  national  incoee,  Y  is  solved  for.  The  solution  is 

Y  «  1  (E  +  i  ♦  G  -  bf  -  bib)  (2.10) 

(I  -  b  -  V) 

Eq  (2.10)  is  known  as  the  reduced  fore  for  Y  since  Y  is  expressed 
in  teres  of  known,  exogenous  quantities.  The  quantity  l/(l-b-t)  is 
known  as  a  multiplier  because  is  tells  how  many  tines  larger  the  change 
in  Y  will  be  in  response  to  a  change  in  Q,  C,  or  1.  For  instance,  if  b  - 
0.70  and  v  =  O.OS  and  G  changes  by  ten  billion  dollars,  then  Y  increases 
by 

(IIP  billion)  ■  $40  billion 
(1-0.70-0.05) 

under  the  assumptions  set  forth  above.  The  multiplier,  l/(l-b-v)  will 
henceforth  be  denoted  as  m.  Eq  (2.10)  makes  the  effects  of  changing 
exogenous  variables  clear.  For  instance,  if  government  expenditures 
increase  by  AG,  Y  will  increase  by  mAG.  On  the  other  hand,  if  T  is 
decreased  by  AT,  Y  will  increase  by  -bmAT.  Since  b  is  less  than  one,  bm 
is  less  than  m.  An  increase  in  government  spending  changes  total  income 
by  more  that  the  same  decrease  in  taxes.  T  is  a  lump  sum  tax.  A  tax  rate 
function  could  be  introduced  into  the  model,  but  it  is  not  necessary  for 
this  discussion  and  will  be  omitted  for  simplicity.  Changes  in  Y  due  to 
changes  in  the  autonomous  components  of  C  and  I  are  also  easily 
determined.  Figure  2.2  shows  the  increase  in  Y  due  to  an  increase  in  G. 

Increasing  government  spending  shifts  the  aggregate  demand  line 
from  DD  to  DD'which  causes  income  to  increase  from  Y*  to  Y'.  In  fact, 
any  change  which  causes  a  shift  in  C,  I,  or  G  will  shift  the  aggregate 
demand  curve  in  a  similar  manner. 


Figure  2.2.  Income  Effects  of  Increased  Government  Spending 
Changes  in  government  spending  affect  more  in  the  economy  than  just 
the  commodities  market.  Increased  government  spending  increases  demand 
for  output  and  labor  in  the  short  run  and  eventually  raises  prices  and 
interest  rates  in  the  longer  term.  These  effects  are  not  seen  in  the 
simplified  commodities  market  presented  here  because  of  simplyfying 
assumptions.  Host  macroeconomic  models  develop  other  interrelated 
markets  to  model  the  behavior  of  interest  rates,  prices,  money  suply, 
and  labor.  Multipliers  for  larger  systems  are  much  more  complex  than 
the  simple  multiplier  in  Eg  (2.10).  Multipliers  must  capture  the 
effects  of  all  variables  in  the  model. 

Itis  Klsi!ic§aL4fe8C3tc  NafisL 

Econometric  models  like  the  Klcin-Goldberger  model  are  formulated 
to  include  many  variables  in  the  economy.  The  KG  model  consists  of  21 
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sinultaniouf  tquationst  IS  a^  tnhich  art  bthaviaral  and  6  of  which  art 
idtntities  rtlating  variablts.  Bthavioral  tquationt  fuch  at  Eq  (2.3)  and 
(2.b)  havt  parantttrs  which  tutt  bt  aitiaattd  frot  historical  tconotic 
data.  Parattttri  in  tht  KG  lodtl  wtrt  tstitatad  using  tconotic  data 
frot  stvtral  sourcts,  prltarily  tht  Unittd  Statts  Dipartttnt  of 
Conterce.  The  KG  tcdtl  brtaks  the  econoty  into  separate  governetnt, 
corporate,  labor,  and  agricultural  sectors.  The  todtl  characterizes  tht 
dynamic  nature  of  the  econoty  through  a  system  of  lagged  variablts. 

Lagged  variables  give  the  value  of  a  variable  in  previous  years. 

For  instance  investment  lagged  one  year  is  the  value  of  investments  one 
year  prior  to  the  current  year.  Lagged  variables  are  denoted  by  a 
subscripted  negative  number  which  represents  the  number  of  periods  the 
variable  is  lagged.  Investment  lagged  five  years  is  denoted  I-s. 

Lagged  variables  together  with  exogenous  variables  are  know  as  predeter¬ 
mined  variables.  Understanding  lagged  variables  is  essential  to  under¬ 
standing  a  central  feature  of  the  KG  model. 

Table  2.1  lists  the  21  equations  in  the  KG  model  and  Table  2.2 
defines  the  variables  in  the  model.  For  a  detailed  discussion  of  the 
model's  theoretical  development,  one  can  consult  Klein  and  Goldberger's 
original  work  presenting  the  model  (Klein  and  Goldberger,  l?55).  Another 
description  of  the  model  appears  in  Theil's  Econometrics  (Theil,  1971). 

Table  2.1.  The  Klein-Goldberger  Model  (Adapted  from  Goldberger,  19S9:4-7) 

C  =  -22.26+0.55(W,+W3-T„)+0.41(P-Tc-Tn-Sc) 

♦0.34(R.+R2-Tn)+0.26C-,*0.072(L, )-,+0.26Np  (2.1.1) 

I  *  - 1 6 . 7 1 +0 . 70 ( P- Tr- Tn ♦R I ^R 2 -Tr+D ) - I -0 . 07 3K- 1 +0 • 1 4 ( L 2 ) - I  (2.1.2) 
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Sc  »  -3.33+O.72(Pc-Tc)+O.O76(Pc-Tc-Sc)-t-O.O20(S.)-. 

Pc  «  -7.60+0.68P 

D  »  7.25»0.  lQ(K-t-K-t)-»0.044(Q-Ma) 

2 

Nt  »  -1. 40+0. 24(Q-Wa)*0.24(Q-M2)-. +0.291 

(Q-Wa)  »  -26.0B+2. 17th(NM-NQ)+Nc]+0. lA(K+K-» )+2.05t 

2 

W-W-.  »  4.  ll-0.74(Ni.-N«-Nc)+0.32(p-i-p-2)+0.54t 
Fi  *  0 • 32+0 • 0060 ( M"T w “ Tc*Tn "Tr ) ( p /pr ) +0. 8 1 ( F T ) - 1 
Rt(p/pR)  *  “0 i 36+0 • 054 ( W I +W2*Ta+P“Tc“TN "Sc ) (p/pp) 

“0.007t (Wt  +  Ma“TM+P“Tc~T  m~Sc ) (p/pp) ]-»+0.012Fr 
Pp  a  -131. 17+2. 32p 

Li  =  0. 14(M-TM-Tc-Tp-Sc-TR)+76.03(ic-2.0)-‘>-“‘» 

La  »  -0.34+0.26Wt-l,02i.-0.26(p-p-t)+0.6l(L2)-i 
1l  *  2 . 58+0 . 44 ( i ■ ) - 3+0 . 26 ( i ■ 1 -9 

lQ0Ci.-(i.l-i]  =  11.17-0.67Lb 

i. 

K-K-,  a  I-D 
Sb” ( Sb ) -  1  a  Sc 
Wi+Wa+P+Ri+Ra  a  M 
C+I+G“Fi  a  M+Te+D 
h(M/p)N.  a  Wi+Wa 


(2.1.3) 

(2.1.4) 
(2.1.3) 

(2.1.6) 

(2.1.7) 

(2.1.8) 
(2.1.9) 

(2.1.10) 

(2.1.11) 

(2.1.12) 

(2.1.13) 

(2.1.14) 

(2.1.15) 

(2.1.16) 

(2.1.17) 

(2.1.18) 

(2.1.19) 

(2.1.20) 


Q  a  H+Te+D 


(2.1.21) 


► 


r 


k 


Tabli  2.2.  Sloiiary  of  Variabltt  for  thi  Kltin-Goldbtrgir  Model 
(Adapted  froe  Soldberger,  1959:3-6) 

Nee  Brief  Category 

Syebol  Definition 


c 

Coneueption 

Endogenous 

D 

Depreciation 

Endogenous 

F. 

leports 

Endogenous 

Fr 

Fare  exports 

Exogenous 

G 

Governeent  expenditures  and 
exports 

Exogenous 

h 

Hours  of  Mork 

Exogenous 

I 

Investment 

Endogenous 

iu 

Long-term  interest  rate 

Endogenous 

1  a 

Short-term  interest  rate 

Endogenous 

K 

Capital  stock 

Endogenous 

Li 

Household  liquid  assets 

Endogenous 

Lj 

Business  liquid  assets 

Endogenous 

Lb 

Percentage  excess  reserves 

Exogenous 

M 

National  income 

Endogenous 

Nb 

Entrepreneurs 

Exogenous 

No 

Government  employees 

Exogenous 

Ni 

Labor  force 

Exogenous 

Np 

Population 

Exogenous 

Nu 

Employees 

Endogenous 

P 

Nonwage  nonfarm  income 

Endogenous 

Pr. 

Corporate  profits 

Endogenous 

*•  N  , s’  •/  •/  %•  •/ 
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p 

Price  level 

Endoqinous 

Pf 

Import  price  level 

Exoqenous 

P" 

Farm  price  level 

Endoqenous 

Q 

Gross  national  product 

Endoqenous 

Ri 

Faro  income 

Endoqenous 

R:> 

Farm  subsidies 

Exoqenous 

Sb 

Corporate  surplus 

Endoqenous 

Sc 

Corporate  savinqs 

Endoqenous 

t 

Time  trend 

Exoqenous 

Tc 

Corporate  taxes 

Exoqenous 

Te 

Indirect  taxes 

Exoqenous 

Tn 

Nonmaqe  nonfarm  noncorporate 

taxes  (less  transfers) 

Exoqenous 

T, 

Farm  taxes  (less  transfers) 

Exoqenous 

Tu 

Waqe  taxes  (less  transfers) 

Exoqenous 

H 

Waqe  rate 

Endoqenous 

Ui 

Private  Naqe  bill 

Endoqenous 

Ws 

Government  eaqe  bill 

Exoqenous 

The  text 

beloH  briefly  describes  each  equation  in  the  model. 

Eq 

(2.1.1)  is  the  consumption  function.  It  qives 

consumption  as 

function 

of  labor,  corporate,  and  aqricultural  disposable  income.  Thi 

equation 

includes  factors  from  the  money  market  throuqh 

the  household 

liquid  assets  (LJ  term.  Plso  included  is  a  population 

trend. 

Eqs 

(2.1.2)  throuqh  (2.1.S)  model  the  behavior  of 

the  corporate 

sector. 

The  investment  function,  Eq  (2.1.2)  is  similar 

in  form  to  the 

consumption  function.  Thi  investment  function  depends  on  corporate  and 


invenstaant  to  factor!  in  thi  iioncy  aarkot.  Inviitaant  also  dtpandi  on 
capital  Mhich  is  the  accueulation  of  undepreciated  capital  as  stated  in 
Eq  (2.1.16).  Eq  (2.1.3)  links  corporate  savings  to  business  incoee  froe 
the  current  and  previous  year  and  business  suplus  froe  the  previous 
year.  Corporate  savings  are  also  defined  as  the  change  in  business 
suplus  in  Eq  (2.1.17).  Eq  (2.1.4)  relating  corporate  profits  to  nonuage 
nonfare  incoee  is  really  just  an  eepirical  relationship  used  to  close 
the  systee.  Depreciation  naturally  depends  on  the  existing  capital  stock 
in  Eq  (2.1.S).  The  second  tere  in  the  equation  shows  that  depreciation 
increases  when  there  is  a  high  degree  of  capacity  utilization. 

Eq  (2.1.6),  (2.1.7),  and  (2.1.8)  eodel  the  behavior  of  the  labor 
earket.  Eq  (2,1.6)  is  the  labor  denand  function  which  gives  private 
demand  for  labor  as  a  function  of  private  sector  output.  Eq  (2.1.7)  is 
the  production  function  which  shows  how  labor  and  capital  combine  to 
produce  private  output.  Eq  (2.1.8)  is  the  labor  supply  function,  but  is 
expressed  in  terms  of  unemployment.  The  lagged  prices  in  the  equation 
indicate  that  wages  are  slow  to  change  in  response  to  price  changes. 

Foreign  imports  in  Eq  (2.1.9)  increase  when  national  disposable 
income  is  high  and  foreign  prices  are  low. 

Eqs  (2.1.10)  and  (2.1.11)  model  the  agriculture  sector  of  the 
economy.  Eq  (2.1.10)  relates  farm  income  to  domestic  customer  prosperity 
and  foreign  exports.  The  ratio  between  the  general  price  index  and  the 
agricultural  price  index  accounts  for  the  terms  of  trade  of  agriculture. 
Eq  (2.1.11)  relates  farm  prices  to  the  general  price  level. 

Eqs  (2.1.13)  through  (2.1. IS)  comprise  the  money  market.  Eos 
(2.1.12)  and  (2.1.13)  are  the  household  and  business  demand  for  monev 


fquationi.  Tht  dioiand  for  aoney  hat  two  coiponintit  spiculativt  dtaand 


Nhlch  is  rslatid  to  Intsrtit  ratot  and  pricts,  and  transaction  diaand 
which  is  rtlattd  to  incoss.  Both  coapontnts  appear  in  Eqs  (2.1.12)  and 
(2.1.13).  Eqs  (2.1.14)  and  (2.1.15)  show  the  relationship  between  short 
tern  interest  rates,  long  tore  interest  rates,  and  bank  excess  reserves. 
Long  tern  interest  rates  are  nerely  a  weighted  average  of  past  period 
short  tern  interest  rates.  The  percent  change  in  short  tern  interest 
rates  depends  on  bank  excess  reserves  which  are  supposedly  deternined 
outside  the  systen.  Host  equilibriun  nacroecononic  nodels  relate  bank 
excess  reserves  (which  deternine  the  supply  of  noney)  to  other  narkets 
in  the  econony  through  prices.  The  Klein-Goldberger  nodel  does  not. 
Soldberger  suggests  that  this  is  a  deficiency  in  nodeling  the  link 
between  the  noney  narket  and  the  connodities  narket.  When  the  KG  nodel 
was  linked  together  to  nake  extended  period  forecasts  with  excess  bank 
reserves  set  at  a  constant  level,  interest  rates  increased  without 
bound.  This  caused  investnent,  consumption,  and  6NP  to  dininish  to 
zero.  To  remedy  this  deficiency,  Goldbergcr  set  liquid  reserves  (L,  and 
L?)  equal  to  a  constant  for  all  studies  of  the  dynanic  nature  of  the 
nodel.  This  measure  effectively  deletes  the  money  narket  from  the  model 
(Soldberger,  1959:04-95).  Consequently,  interest  rates  oust  be  ignored 
and  the  commodity  and  labor  markets  are  linked  directly  through  output 
and  prices. 

Eqs  (2.1.10)  through  (2.1.21)  are  identities  relating  variables  in 
the  nodel.  Eqs  (2.1.19)  and  (2.1.21)  together  are  KG  model  version  of 
the  national  income  accounting  identity. 

In  order  to  use  the  Klein-Goldberger  model  for  the  current  state  of 
the  economy,  economists  must  solve  the  system  of  simultaneous  eouations 
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in  Tabli  2.1.  Tha  folloMing  sactian  ihONs  hoN  aconoaatric  aodala  ara 


aolvad  and  linkad  togathar  to  provida  aktandad  pariod  foracatts.  In 
addition  tha  saction  discuifai  hoN  tha  foracaata  ara  uaad  to  coaputa 
nultipliara. 

l&8Q8!etCic  (!8dala  an^  £8fl88tl!!8  QultlBlieCf 

Hacroaconoaatric  aodala  ara  claaai'fiad  aa  linaar  or  nonlinaar. 

Linaar  aodala  ara  solved  with  diHarant  aathoda  than  nonlinear  aodala. 

Tha  K6  aodal  is  an  exaaple  o^  a  nonlinaar  aodal.  Nonl inearitiaa  appear 

Mhen  current  endogenous  variables  ara  raised  to  powers  other  than  one  or 

are  aultiplied  together.  Nonl inear itics  appear  in  K6  aodal  Eqs  (2.1.7), 
+ 

(2.1.?),  (2.1.10),  (2.1.12),  (2.1.15),  and  (2.1.20). 

iSSQSSSiEiS  SS^bIs.  An  exaaple  of  a  linaar  aodal  is 
the  Klein  Model  I  shown  in  Table  2.3  below. 

Table  2.3.  Tha  Klein  Model  1  Linear  Econoaetric  Modal 
(Adapted  froa  Theil,  I97l! 432-435) 

C  *  16. 78+0. 02+0. 23P-. +0. 0O(W+W )^a 
I  »  17.79+0.23P+0.55P-.-0. 15K-.+e' 

W  =  1.60  +  0.42X  +  0. 16X-,+0.  13(T-1931)+e' ' 

X  =  c+i+e 

P  =  X-W-T 
K  =  K-i+I 

where 

C  a  consuaption, 

P  =  profits, 

W  3  mge  bill  paid  by  private  industry. 


I  3  net  invastaant. 


K  s  capital  itock, 


,  X  -  total  production  ot  privatt  induttry,  and 

'  And  •*'  >  randoa  trror  tarai. 

h 

t 

I  The  eodel  predict!  current  endogenous  variable  values  given  current 

exogenous  variables  and  lagged  endogenous  values. 

Linear  eodels  can  be  solved  Mith  eatrix  algebra.  Any  linear 
I  econoeetric  eodel  can  be  put  into  the  fora 

]  Gy  ♦  Bz  ■  E 

where 

I  y  *  the  a  elaaant  coluan  vector  of  a  endogenous  variables, 

'  G  a  the  a  X  a  coefficient  aatrix  with  a  coefficient  for  each  of  a 

» 

endogenous  variables  for  each  of  the  a  equations 
I  z  «  the  n  eleaent  coluan  vector  of  n  predeterained  (lagged 

f 

endogenous,  exogenous,  and  lagged  exogenous)  variables, 

B  -  the  a  X  n  coeficient  aatrix  with  one  coefficient  for  each  of  n 
^  predeterained  variables  for  each  of  the  a  equations,  and 

E  »  the  a  eleaent  error  vector. 

I  For  the  Klein  Model  I, 

y"  =  CC,  P,  W,  I,  K,  X] 


'1.00 

o 

o 

I 

1 

o 

CD 

o 

0.00 

0.00  0.00' 

0.00 

>0.23 

0.00 

1.00 

0.00  0.00 

G’’  » 

0.  00 

0.00 

1.00 

0.00 

0.00  -0.42 

l.OO 

0.00 

0.00 

1.00 

0.00  -1.00 

0.00 

1.00 

1.00 

0.00 

0.00  -1.00 

0.00 

0.00 

0.00 

1.00 

-1.00  0.00. 

z’’  = 

Cl,  P 

-  1  •  K  - 

It  X  - 1 

,  W  , 

T,  G,  t] 

I 
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•  .*  V 


-16.78 

-0.23 

0.00 

0.00 

-0.80 

0.00 

0.00 

0.00' 

-17.79 

-0.55 

0. 15 

0.00 

0.00 

0.00 

0.00 

0.00 

249.43 

0.00 

0.00 

-0.16 

0.00 

0.00 

0.00 

-0. 13 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Ce,  e', 

e"  , 

o 

o 

o 

0.00, 

0.00] 

For  the  deterministic  case,  it  can  be  assumed  that  E  -  0.  Using  matrix 
algebra  one  can  solve  For  y. 


y  =  -S-‘Bz 

letting  D  =  Eg  (2.12)  may  be  expressed  as 

y  »  Dz 


(2.12) 


(2. 13) 


Eg  (2.13)  gives  the  unknown  endogenous  variables  in  terms  of  linear 
functions  of  known  predetermined  variables.  The  m  by  n  0  matrix  contains 
constants  which,  when  multiplied  by  a  predetermined  variable,  give  the 
level  of  an  endogenous  variable.  These  constants  are  known  as 
multipliers.  Multipliers  will  be  discussed  more  fully  later. 

The  Klein  Model  I  forecasts  only  the  present  period  from  past 
periods  and  current  policy.  It  would  be  useful  to  forecast  a  future 
period  based  on  the  current  state  of  the  economy  and  the  expected 
external  forces  influencing  the  economy.  Figure  2.3  shows  how  the  state 
of  the  economy  evolves  from  period  to  period  as  influenced  by  exogenous 
forces. 
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Figure  2.3.  The  Evolving  Economy 


The  solution  for  an  extended  forecast  can  be  found  by  decomposing 
Eg  (2.13)  into  its  components.  For  example,  in  the  case  Nhere  there  are 
only  one  period  lags,  Eg  (2.13)  can  be  rewritten, 

y  =  do  +  D,y-i  t  DjXo  D:,x-, 

where 

yo  =  vector  of  current  endogenous  variables  to  be  solved  for, 
y-i  =  vector  of  endogenous  variables  lagged  one  year, 

Xo  =  vector  of  current  exogenous  variables, 

x-i  3  vector  of  exogenous  variables  lagged  one  year,  and 

do.  Di .  D?.  and  Dx  =  coefficient  matricies. 


Ds  contains  nultipliers  which  give  the  iepact  o^  current  exogenous 
variables  on  current  endogenous  variables  and  are  known  as  inpact  multi¬ 
pliers. 

It  is  assumed  that  xo  and  x-i  are  determined  outside  the  system  by 
uncontrollable  events  or  by  policy.  However,  y-,  is  the  result  of 
previous  period  activity.  Writing  this  explicitly, 

yo  -  d o +  0 1  ( do +D 1 y - a +  D aX -  1 +D3X - 2 ) ♦DaX o ♦DsX -  1 

®  ( I ♦ D I ) d o +0 1 ^ y - 2+D2 X 0+ ( D 1 D2+D3 ) X - 1 ♦  D1D3  (2.14) 

By  decomposing  the  y-a  into  its  components,  one  can  express  the  current 
state  of  the  economy  in  terms  of  exogenous  variables  and  endogenous 
variables  lagged  three  periods.  The  process  is  repeated  to  obtain  fore¬ 
casts  for  any  number  of  periods  in  the  future. 

The  coefficient  matricies  for  the  lagged  exogenous  variables 
Ce.g.,x-i,  x-2  in  Eg  (2.14)3  are  especially  important  for  policy  analy¬ 
sis  because  they  give  the  impact  of  exogenous  variables,  including 
fiscal  and  monetary  variables.  They  are  different  for  each  lagged  period 
Ce.g.,  Da,  (DiDa+Ds)  in  Eq  (2.14)3  indicating  the  changing  influence  of 
the  exogenous  variable  over  time.  The  numbers  in  the  Da  and  DiDa^D3 
matricies  are  known  as  interim  multipliers. 

ISliLHa  NoQlLQSiE  i£9QS(S®Tr i.c  Models,  flultipliers  computed  by 
matrix  algebra  completely  characterize  both  dynamic  and  static  aspects 
of  a  linear  system.  They  give  the  impact  of  exogenous  changes  on  the 
state  of  the  economy  at  any  point  in  time.  Unfortunately,  the  national 
economy  cannot  be  accurately  describee  by  a  linear  system.  Economic 
theory  prescribes  inherently  nonlinear  functions.  For  example,  nominal 
endogenous  variables  are  often  divided  by  a  price  index  to  obtain  real 


values.  The  price  index  Is  itself  an  endogenous  variable,  and  so  equa¬ 
tions  containing  noeinal  values  converted  to  real  values  are  nonlinear. 
Another  exanple  of  inherent  nonlinearity  is  production  functions.  Func¬ 
tions  such  as  the  Cobb-Douglas  production  function  and  the  constant 
elasticity  of  substitution  production  function  have  proven  to  charac¬ 
terize  real  world  econoeics  quite  well.  They  are  nonlinear. 

Solving  and  analyzing  nonlinear  systees  is  not  as  simple  as  solving 
linear  systems.  Nonlinear  systems  cannot  be  solved  by  simple  matrix 
algebra.  Usually  some  numerical  technique  must  be  used.  However,  if  the 
system  is  aggr gx i matel y  linear,  a  derivative  technique  may  be  used  to 
linearize  the  system.  Goldberger  used  this  technique  to  linearize  the  K6 
model.  The  derivations  which  follow  are  a  condensation  of  Goldberger's 
work.  (Goldberger,  1959:17-20) 

A  Nonlinear  equation  in  a  single  y  and  a  single  z  can  be  written  in 
the  form, 

f (y,z)  =  0 

The  total  differential  of  the  function  f  is  also  equal  to  zero. 

df  »  l<f/i»y)dy  (of/'zldz 

Solving  for  dy, 

dy  =  -(df /3z)dz  (2.15) 

(df /dy) 

Eq  (2.15)  gives  the  explicit  dependence  of  changes  in  y  on  changes  in  z. 
The  expression 


nay  not  be  a  constant,  but  if  sanple  neans  are  substituted  for 
variables,  the  expression  can  be  evaluated  at  a  point.  If  the  equation 
is  approximately  linear,  this  constant  Mill  be  approxinately  correct  for 
a  large  range  of  z.  Next,  consider  the  nonlinear  system  of  equations 
written  in  the  operator  form, 

F(y,z)=0 

where  F  is  a  matrix  of  functional  operators.  Taking  the  total 
differential  and  solving  for  dy, 

dF  =  (SF/3y)dy  ♦  (^F/^zldz  »  0 
-(f^F/^y)  =  -(^F/^z)dz 

dy  *  -<dFA^y)-‘ (c»F/3z)  (2.16) 

Eq  (2.16)  is  like  the  solution  to  the  linear  system  in  Eq  (2.12) 
except  it  is  expressed  in  terms  of  differentials.  Also,  ^F/^y  and  ^F/(^z 
are  not  always  constant  matricies  but  can  be  function  matricies.  By 
evaluating  these  matricies  at  some  value,  say  at  the  sample  mean  of  each 
predetermined  variable,  these  matricies  can  be  converted  to  constant 
matricies.  The  elements  of  the  constant  matricies  are  multipliers  for 
changes  in  predetermined  variables.  In  general,  they  are  guaranteed  to 
be  valid  only  for  small  changes  about  the  point  at  which  they  are 
evaluated.  However,  if  the  system  is  approximately  linear,  the  multi¬ 
pliers  may  be  approximately  correct  over  a  wide  range  of  values.  If  so, 
they  may  be  used  in  a  manner  similar  to  the  D  matricies  computed  for 
linear  systems.  Extended  period  forecasts  and  interim  multipliers  are 


conputtd  in  a  nanner  analaqous  to  tha  lintar  aconoMatric  aodal.  Untortu- 
nataly,  ona  cannot  alNayi  count  on  tha  aconoaatric  nodal  baing  avan 
approxinatal y  linaar. 

H  the  nodal  cannot  ba  linearized,  than  a  nunerical  technique  for 
solving  tha  nodal  is  usually  used.  Thera  are  several  nunerical  tech¬ 
niques  available  including  Neaton-Raphson,  Gauss-Sidel,  and  others.  The 
nethod  chosen  for  this  study  is  the  Sauss-Sidel.  Klein  reconnandad  tha 
nethod  over  the  Neaton-Raphson  nethod  because  although  the  Neaton- 
Raphson  nethod  usually  converges  in  fewer  iterations,  each  iteration 
requires  significantly  nore  conputation  than  each  iteration  of  the 
Gauss-Sidel  nethod  (Klein,  1974:238-240).  The  Gauss-Sidel  nethod  it  alto 
easy  to  progran  and  debug  and  does  not  require  tha  conputation  of  a 
derivative.  Appendix  A  describes  tha  Gauss-Sidel  nethod  in  detail. 

Solving  the  nonlinear  econonetric  nodel  for  current  endogenous  in 
terns  of  predeternined  variables  produces  a  forecast  of  current  endoge¬ 
nous  variables.  To  be  even  nore  useful,  a  nethod  nust  be  devised  to 
produce  extended  period  forecasts.  Extended  period  forecasts  can  be 
conputed  fron  current  period  forecasts  by  setting  lagged  endogenous 
variables  equal  to  the  current  endogenous  solution,  updating  the  exoge¬ 
nous  variables,  and  resolving  the  systen.  For  exanple,  after  each  solu¬ 
tion  is  conputed,  current  values  of  consunption,  investnent,  etc.,  are 
deternined.  To  extend  the  forecast,  the  consunption,  investnent  etc., 
variables  lagged  one  year  are  set  equal  to  the  current  solution  for 


consunption,  investnent,  etc.  Lagged  exogenous  variables  are  updated  in 
a  sinilar  nanner.  Current  exogenous  variables  are  set  to  whatever  the 


policy  under  investigation  requires.  Then  the  nodel  is  resolved.  In  this 


M«y,  tht  nodtl  can  bi  linkad  taqethtr  to  obtain  foricatti  for  any  nuabtr 
of  periods  in  the  future. 

Because  eultipliers  are  so  valuable  for  policy  analysis,  a  eethod 
of  coeputing  multipliers  for  models  salved  by  numerical  techniques  is 
needed.  Evans  and  Klein  describe  a  eore  general  method  of  detereining 
multipliers  which  can  be  used  with  models  solved  by  numerical  techniques 
(Evans  and  Klein,  1968:48-49).  To  calculate  the  multipliers,  a  control¬ 
led  solution,  yc,  is  computed  with  all  predetermined  variables  (exoge¬ 
nous  and  lagged  endogenous)  at  a  given  level,  and  with  the  input 
variable  of  interest  set  at,  say,  Xe.  Next,  a  new  solution,  yd,  is 
computed  with  x  at  a  disturbed  level,  Xd.  The  multiplier,  e,  is  then: 


a 


(2.17) 


A  generalized  multiplier  such  as  m  can  be  computed  for  a  “package"  of 
changes  in  predetermined  variables.  However,  they  are  valid  only  for  the 
changes  and  variable  levels  used  to  estimate  them.  A  separate  run  for 
each  combination  of  input  variable  changes  must  be  made  to  estimate  each 
multiplier. 

tlGlLtations  of  Econometric  Models  and  M:ilti.ELi8r  AnaLySlS. 

The  limitations  of  econometric  models  and  multiplier  analysis  are 
summarized  below. 

1.  Linear  econometric  models  can  be  constructed  which  are  easily 
salved  and  analyzed:  however,  they  do  not  accurately  reflect 
the  national  economy  in  theory  or  in  practice. 

Near  linear  models  more  accurately  predict  the  performance  of 
the  national  economy,  and  they  may  be  anlayzed  with  minor 
accuracy  degadation,  but  they  may  not  adequately  model  the 
inherent  nonlinearities  of  the  actual  economy. 


2. 


3.  Lets  aggragatid,  ADra  nonlintar  aodels  aay  bt  daviaad  which 
accurately  predict  the  econoeyi  but  euch  eodeli  are  difHcult 
to  analyze.  Multipliers  eay  be  calculated  by  coeputing  control 
and  disturbed  solutions  and  then  dividing  the  diHerence  in 
these  two  solutions  by  the  diHerence  between  input  variables. 
However  these  eultipliers  are  good  only  for  saall  changes  about 
the  specific  disturbed  solution  for  which  they  were  coeputed. 

The  model  must  be  rerun  for  each  policy  alternative  is  examined. 


Response  surface  methodology  is  one  way  of  overcoming  some  of  these 
difficulties.  Response  surface  methodology  accomodates  nonlinearities. 

By  fitting  a  response  surface  to  an  econometric  model,  one  could 
investigate  the  effects  of  varying  one  or  more  key  input  variables 
singly  or  jointly  over  their  entire  ranges.  The  next  chapter  describes 
response  surface  methodology  and  explains  how  it  might  be  applied  to 


econometric  models. 


III.  Agel^ing  Response  Surfggg 


Intro^aStign 

Ch«pt*r  I  dtfcrlbi*  what  rasponit  lurfaci  ■■thodology  ii  and  hoN 
it  can  ba  usad.  Thii  chaptar  diacuatai  tha  stapt  in  applying 
rasponaa  lur^aca  aathodology  to  a  problaa  in  ganaral  to  provida  a 
background  for  tha  aathodology  davalopad  in  Chaptar  IV.  This  chaptar 
also  suggasts  spacific  Mays  to  apply  rasponsa  surfaca  aathodology 
to  tha  analysis  of  aconoaatric  aodals. 

Applying  ASM  can  be  divided  into  eleven  distinct  steps.  The  steps 

are 

1.  Define  the  problea  and  determine  that  RSn  is  an  appropriate 
analysis  technique. 

2.  Determine  tha  input  and  output  variables  of  interest. 

3.  Determine  tha  operating  region  of  interest. 

4.  Select  a  response  surface  equation. 

5.  Select  an  experimental  design. 

6.  Translate  the  coded  design  points  to  actual  factor  levels. 

7.  Run  the  experiment  or  model  to  obtain  responses  for  each  set  of 
factor  levels. 

8.  Regress  the  coded  experimental  design  on  the  responses. 

9.  Check  the  response  equation  fit. 

10.  Decode  the  response  surface  coefficients. 

11.  Perform  analysis  on  the  fitted  response  surface. 

The  discussion  below  amplifies  each  step. 

The  first  step  in  applying  RSM  is  to  define  the  problem  and  to 
decide  that  RSH  is  an  appropriate  method  for  analysis.  Not  all  problems 


lend  theneelvis  to  antlyeis  by  RSn.  RSn  relates  multiple  Inputs  to  a 
single  output.  Meyers  'further  points  out  fundaeental  assumptions  under¬ 
lying  RSM  in  his  text,  Resgonse  Surface  Methodology  (Meyers,  1976:62). 
RSM  is  appropriate  for  problems  in  which  the  relationship  between 
input  variables  and  output  variables  is  either  very  complex  or  unknown, 
but  the  variables  are  quantitative  and  continuous.  Also,  the  'functional 
relationship  between  inputs  and  the  response  must  be  approximated  by  a 
low  order  polynomial  or  other  simple  function  whose  parameters  are 
estimable.  Finally,  the  input  variables  must  be  controllable  and  all 
variables  must  be  measured  with  negligible  error. 

Deterraini.ng  Vari.abl.es  of  Interest 

The  second  step  in  applying  RSM  is  to  determine  the  input  and 
output  variables  of  interest.  The  input  variables  selected  for  the 
analysis  must  include  all  the  important  factors  which  bear  on  the 
problem.  Properly  defining  the  problem  should  make  these  important 
factors  obvious.  However,  the  size  of  the  experimental  design  reguired 
to  estimate  response  function  coefficients  increases  rapidly  with  the 
number  of  factors  (more  on  this  below).  All  important  factors  should  be 
included  in  the  response  surface  equation,  but  the  number  of 
experimental  design  points  required  must  be  considered. 

5§t§C!DiDiD9  tdi  QB§C§tiQ9  BS9i9D 

Once  the  variables  of  interest  have  been  selected,  their  ranges 
must  be  specified.  The  ranges  of  the  input  variables  must  be  feasible 
and  independent  of  one  another.  In  addition,  the  ranges  should  be  narrow 
enough  so  that  the  response  does  not  contain  too  manv  inflection  points. 
Too  many  inflection  paints  in  the  response  require  a  complicated 
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rRsponse  function  with  higher  order  teres  and  a  large  experieental 
design  to  accurately  capture  the  input  output  relationship. 

Selecting  the  Resggnsg  Surfact  IflgsiLSQ 

The  response  surface  equation  selected  in  the  fourth  step  is 
usually  a  first  or  second  degree  polynoeial.  An  exaeple  of  a  first 
degree  polynoeial  response  surface  equation  with  two  input  variables  is 

y  *  Bo  f  BiXi  ♦  B2X2 

and  an  example  of  a  second  degree  polynoeial  response  surface  equation 
Mith  tMO  input  variables  is 

y  *  Bo  f  BiiXi^  f  BtXi  +612X4X2  +82X2  f  822X2^ 

where 

y  s  output  (response)  variable, 
xi,  X2  •  input  variables,  and 

Bi  -  response  surface  coefficients  to  be  deternined 
(i  =  0,  I,  2.). 


There  are  a  number  of  advantages  to  using  a  low  order  polynomial  as 
a  response  surface  equation.  First,  the  coefficients  of  a  polynomial  are 
estimable  by  the  method  of  least  squares,  the  most  commonly  used 
regression  technique.  Also,  experience  has  shown  that  a  first  or  second 
degree  polynomial  worxs  well  as  a  response  surface  function  because  a 
polynomial  is  a  truncated  form  of  the  Taylor  series  (heyers,  1976:62). 
Another  argument  for  the  use  of  low  order  polynomials  is  that  many 
experimental  designs  have  been  developed  for  collecting  data  to  estimate 
the  coefficients  of  polynomials.  However,  functions  other  than 


polynomials  can  serve  as  the  response  surface  equation.  Theoretical 


considirations  may  dlctati  tha  utc  of  a  certain  type  of  eathenatical 
function.  Such  functions  eay  approxieate  the  response  very  closely  and 
should  be  used;  hoMever,  these  functions  should  reeain  sieple  so  that 
aspects  of  the  underlying  process  can  be  easily  explored  and  interpreted 
(Hill  and  Hunter,  1966:573). 

iSlS£^i.Q9  iHBiCLSSQtil  QSSLSQ 

The  next  step  in  applying  RSH  is  selecting  an  experimental  design. 
An  experimental  design  is  a  set  of  specifications  of  input  variable 
levels  for  repeated  experimental  runs  of  the  process  under  study.  Each 
combination  of  input  levels  is  called  a  design  point.  Table  3.1  contains 
an  example  of  a  three  level  three  factor  experimental  design  with  27 
design  points. 

The  design  in  Table  3.1  is  in  coded  form.  Factor  levels  are 
represented  by  1,  0,  and  -1  for  three  level  experimental  designs  and  1 
and  -1  for  a  two  level  design.  For  a  three  level  design,  a  1  represents 
the  factor  high  level,  a  -1  represents  a  low  factor  level,  and  a  0 
represents  the  average  of  the  high  and  low  values.  In  a  two  level 
design  the  1  represents  the  high  level,  and  the  *1  represents  the  low 


Table  3.1.  Three  Factor  Three  Level  Factorial  Experleental 

Design 


Factor  1 

Factor  2 

Factor  3 

-1 

-1 

-1 

-1 

-1 

0 

-1 

-1 

1 

-1 

0 

-1 

-1 

0 

0 

-1 

0 

1 

-1 

1 

-1 

-1 

1 

0 

-1 

1 

1 

0 

-1 

-1 

0 

-1 

0 

0 

-1 

1 

0 

0 

-1 

0 

0 

0 

0 

0 

1 

0 

1 

-1 

0 

1 

0 

0 

1 

1 

1 

-1 

-1 

1 

-1 

0 

1 

-1 

1 

1 

0 

-1 

1 

0 

0 

1 

0 

1 

1 

1 

-1 

1 

1 

0 

1 

1 

1 

By  using  an  appropriate  experimental  design,  one  can  estimate  the 
coef 'f icients  in  the  response  function  with  a  minimum  number  ai  experi¬ 
mental  runs.  Economy  runs  is  an  important  criteria  for  choosing  an 
experimental  design.  However,  the  minimum  number  of  runs  required  to 
accurately  estimate  response  surface  coefficients  depends  on  the  type  of 
response  surface  equation.  Equations  with  higher  powers  of  variables 
and  interaction  terms  (products  of  input  variables)  require  more  runs  to 
estimate  the  unknown  coefficients.  As  noted  in  the  paragraph  on  selec¬ 
ting  variables  of  interest,  the  size  of  the  experimental  design  limits 
the  number  of  variables  in  the  response  surface  equation.  For  example, 


if  k  is  tht  nusbtr  oi  fsctors,  a  three  level  factorial  experieental 
design  requires  3**  experieental  runs  to  estieate  the  coefficients  of 
nain,  inteaction,  and  squared  factor  effects.  A  3^  factorial  design 
requires  729  experieental  runs.  If  it  is  knoen  that  soee  coefficients 
are  insignificant,  soee  experieental  runs  nay  be  elieinated.  In  addi¬ 
tion,  soee  runs  nay  be  elieinated  at  the  expense  of  having  sone  of  the 
variation  attributed  to  the  wrong  tern. 

Attributing  sone  of  the  variation  to  the  wrong  tern  is  caused  by 
nulticolinearity  in  the  experieental  design.  HuI ti col ineari ty  occurs 
because  of  correlation  between  input  variables.  Ideally,  the  experinen- 
tal  design  should  be  orthogonal,  which  neans  that  there  is  zero  correla¬ 
tion  between  input  variables,  host  orthogonal  designs  require  nueerous 
design  points.  Consequently,  selecting  an  experieental  design  involves  a 
tradeoff  between  econoey  of  experieental  runs  and  orthogonality.  Box  and 
Benkhen  devised  soee  three  level  designs  which  eake  a  very  reasonable 
tradeoff  between  orthogonality  and  econoey.  These  designs  can  be  found 
in  the  paper  entitled  "Soee  New  Three  Level  Designs"  (Box  and  Behnken, 
1960:460-463) . 

One  additional  point  worth  eentioning  is  that  eost  experieental 
designs  are  devised  for  experieents  in  which  there  is  randoe  variation 
in  the  response  due  to  uncontrollable  factors.  They  contain  extra 
points,  usually  center  points,  to  estimate  the  size  of  this  randoe 
variation.  If  the  response  has  no  random  variation,  then  these  extra 
points  are  redundant  and  may  be  eliminated.  Work  in  this  study  assumes 
that  there  is  no  random  variation  in  the  Kl ei n-Gol dberger  econometric 
model  (i.e.,  it  is  a  deterministic  model),  and  so  repetitive  center 


points  are  not  needed. 


Decoding  the  Exger ieental  Desi.gn 

Before  the  experieent  can  be  run  to  obtain  data  for  estieating  the 
response  surface,  the  experieental  design  nust  be  decoded  froe  I's,  O's, 
and  -Is  to  actual  factor  levels.  The  range  of  the  variables  of  interest 
determines  the  factor  levels  for  the  cxperiecnt.  For  example,  if  the 
range  of  a  variable  is  10  to  60,  then  a  1  in  the  experimental  design 
represents  a  factor  level  of  60,  a  -1  represents  10  and  a  0  represents 
the  average  of  the  two  factors,  i.e.,  35. 

BmODIDB  |xB§!Ii!!l5Ctal  Design 

To  obtain  the  input  and  output  data  necessary  to  estimate  the 
response  surface  coefficients,  the  experiment  must  be  run  at  the  levels 
specified  in  the  experimental  design  and  the  response  recorded.  This 
step  IS  straight  forward.  When  running  an  actual  experiment  with  random 
variation,  it  is  advisable  to  randomize  the  order  in  which  the  experi¬ 
ments  are  run.  However,  when  obtaining  data  from  a  deterministic 
mathematical  model,  order  is  unimportant. 

EiiiiDS  idf  BfSESDse  Surface 

Once  an  appropriate  experimental  design  has  been  selected  and  the 
experiment  run  to  collect  data,  the  response  surface  must  be  fit  to  the 
data.  The  surface  is  fit  to  the  data  by  computing  the  response  surface 
coefficients.  The  method  of  least  squares  regression  is  the  usual  method 
for  fitting  a  surface  to  data.  This  method  computes  coefficients  which 
minimize  the  distance  from  the  observed  responses  to  the  response  sur¬ 
face.  For  a  more  thorough  discussion  of  least  squares  estimation,  one 
can  consult  a  statistics  textbook.  One  excellent  source  is  Mathematical 
iiiSiSiiSS  SJiTti  BBEiiSItigns  by  Mendenhall,  Scheaffer,  and  Wackerlv 
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(nendenhall,  Scheaffer,  and  Nackerly,  1981:425).  Rtsponfc  surface  coef¬ 


ficients  are  estinated  by  regressing  the  coded  experieental  design 
matrix  on  the  response  variable.  Using  the  coded  design  matrix  preserves 
orthogonal i ty . 

Checking  Fit 

Once  the  response  surface  has  been  constructed,  it  must  be  checked 
for  proper  fit.  There  are  at  least  four  ways  to  check  the  fit  of  the 
response  surface.  They  include  checking  the  values,  checking  the  sum 
of  squares  error  (SSE) ,  checking  the  residuals  of  the  design  points,  and 
checking  the  residuals  of  random  points.  The  R^  value  gives  the  fraction 
of  total  variation  explained  by  the  response  surface  equation.  R^ 
always  increases  with  the  number  of  factors  in  the  response  surface 
equation.  SSE  gives  much  the  same  information  as  R^,  except  it  may 
increase  with  the  number  of  factors  after  a  certain  point.  The  R^  and 
SSE  criteria  are  the  easiest  and  quickest  way  to  check  for  fit. 

Another  way  to  check  fit  is  to  examine  residuals.  By  dividing  the 
residual  by  the  actual  resoonse  value,  a  measure  of  the  error  can  be 
computed.  These  errors  can  be  averaged  for  all  the  design  points  and 
then  subtracted  from  one  to  give  a  value  similar  to  R‘.  Mathematically, 
this  relation  is 


percent  fit  =  100  7.  1 


1 

n 


2 

1 


where 

y,  s  macroeconometr 1 c  model  output  for  input  combination  i, 
v>  =  response  surface  output  for  input  cpmbination  i,  and 


n  =  number  input  combinations. 


It  is  perhaps  even  more  useful  to  look  at  the  largest  errors  to 
determine  where  the  greatest  lack  o^  fit  occurs.  Residual  plots  are  also 
helpful  in  determining  where  lack  of  fit  occurs. 

As  a  final  check  of  the  response  surface  fit,  experiments  with 
random  factor  levels  can  be  run  and  the  responses  compared  with 
responses  predicted  by  the  response  surface.  An  error  measure  similar  to 
the  error  measure  based  on  residuals  can  also  be  computed.  The  advantage 
of  checking  random  points  is  that  it  may  reveal  anomalies  in  the 
response  which  were  missed  by  the  experimental  design. 

The  four  measures  of  response  surface  fit  mentioned  here  are  used 
to  determine  whether  or  not  the  response  surfce  fits  well  enough  for  the 
purpose  intended.  If  the  response  surface  does  not  fit  well  enough,  a 
new  response  surface  equation  is  usually  postulated.  Steps  4  through  9 
are  repeated  until  the  fit  meets  requirements. 

5S£9^LQa  Its  Resagnse  Surface  Iguatign 

Qnce  an  acceptable  response  surface  is  obtained,  the  response 
surface  equation  is  decoded.  The  equation  requires  decoding  because  the 
response  equation  computed  in  step  eight  was  computed  from  coded  input 
variables.  To  make  the  response  surface  equation  i nterpretabl e  it  must 
be  expressed  in  terms  of  the  original  variables.  Let  x>„.«  be  the  high 
factor  level  for  the  ith  factor,  x....r.  be  the  low  factor  level  for  the 
ith  factor  level,  and  ;< .  o  be  the  decoded  factor  level  for  the  ith 
factor,  and  x.c  be  the  coded  factor  level  for  the  ith  factor.  By 
substituting  the  expression 


in  the  coded  response  equation  ior  the  ith  coded  t/ariable  and  collecting 
terms,  the  equation  is  decoded.  Decoding  equations  in  this  manner  by 
hand  is  tedious  and  prone  to  error.  To  make  the  job  easier,  a  second 
order  polynomial  response  surface  equation  can  be  decoded  using  matrix 
algebra.  (See  Appendix  C) 

BD§iySiDg  the  Response  Surface 

The  final  step  in  applying  RSli  is  analyzing  the  response  surface. 
The  methods  available  for  analysis  are  discussed  in  detail  in  the 
section  below.  But  first,  a  few  prelimiary  comments  are  in  order. 
Obtaining  a  good  response  surface  equation  fit  implies  that  the  response 
surface  is  an  accurate  representation  of  the  underlying  process  or 
model.  More  faith  can  be  placed  in  the  validity  of  the  response  surface 
if  the  check  of  random  points  confirms  a  good  fit.  However,  the  response 
surface  is  only  valid  for  input  variable  values  within  the  ranges 
originally  used  to  estimate  the  response  surface.  A  response  surface 
fit  to  a  deterministic  model  may  fit  well,  but  the  response  surface  is 
only  as  good  as  the  underlying  model  used  to  construct  it.  If  the  model 
does  not  capture  the  process  it  is  supposed  to  represent,  the  response 
surface  will  not  either. 

There  are  several  well  developed  uses  for  response  surfaces.  The 
most  commonly  useo  technique  is  optimization.  The  explicit  form  of  the 
response  surface  equation  giving  the  response  variable  as  a  function  of 
known,  controllable  input  variables  lends  itself  to  optimization  prob¬ 
lems.  Moreover,  Granev  showed  how  several  response  surfaces  could  be 
comPined  for  constrained  optimization  problems  (Graney,  1984).  Another 
use  of  response  surfaces  is  for  performing  "what  if"  analysis.  One  can 


determine  the  impact  of  changing  input  variables  individually  or  jointly 
without  running  the  model  or  experiment.  In  addition,  tradeoffs  between 
factors  can  be  displayed  to  decision  makers  graphically  for  valuable 
insights  into  the  problem.  Response  surfaces  can  also  be  used  to  make 
predictions.  One  must  remember  that  the  predictions  made  are  only  valid 
if  independent  variables  not  in  the  response  surface  equations  actually 
assume  the  levels  that  they  were  at  when  the  response  surface  was 
estimated. 

ftfifilicat i.gns  for  tlicrgecgngmetr i c  Mgdels 

These  uses  suggest  several  applications  for  macroeconomic  models. 
Multiplier  analysis  completely  captures  the  input-output  relationship  of 
linear  models  well,  but  it  cannot  handle  nonlinear  models  as  well.  RSM 
should  be  able  to  capture  input-output  relationships  in  nonlinear  models 
for  key  variables  easily.  Response  surfaces  describing  key  economic 
performance  variables  such  as  gross  national  product,  inflation,  and 
unemployment  could  be  used  to  assess  the  impact  of  changing  fiscal  or 
monetary  policy  variables  such  as  government  spending,  taxes,  and  money 
supply.  The  response  surface  could  be  used  to  answer  "what  if"  questions 
in  policy  simulation.  Optimal  policy  for  obtaining  specific  economic 
goals  could  be  determined.  All  of  this  analysis  could  be  done  using  a 
limited  number  of  experimental  runs.  A  response  surface  captures  the 
relationship  between  variables  of  interest  over  the  entire  region  of 
interest.  The  controldisturbed  solution  method  of  computing  multipliers 
only  characterizes  how  specific  changes  in  input  variables  affect  output 
variables.  Every  time  a  new  combination  of  policies  is  considered,  the 
model  must  be  rerun  and  a  new  multiolier  comouted. 

This  chapter  describes  the  general  procedures  for  generating 


response  surfaces  and  suggtsts  Mays  of  applying  RSM  to  the  analysis  of 
nacroecononic  nodels.  Chapter  IV  describes  how  the  general  application 
steps  are  actually  applied  to  the  Klein-Goldbergcr  econonetric  nodel  in 
this  study. 


IV.  Mfthodolog^ 


Introduction 

This  thtsis  ef'fort  proposes  to  bring  together  concepts  fron  two 
areas,  response  surface  methodology  and  macroecononetric  modeling. 
Chapter  I  lists  the  specific  objectives  to  be  met  by  this  research 
effort.  They  include  determining  whether  RSli  can  accurately  fit  a 
macroeconometric  model,  reproducing  multipliers  for  the  macroeconometric 
model  via  RSM,  interpreting  the  response  surface  equations,  and 
developing  RSn  applications  for  the  response  surface  equations.  This 
chapter  develops  the  methodology  by  which  these  objectives  are  achieved. 

To  accomplish  the  research  objectives,  a  scheme  must  be  developed 
for  generating  and  checking  response  surface  equations  from  the  Klein- 
Soldberger  model.  The  scheme  used  here  follows  the  general  steps  for 
applying  RSli  to  any  problem  as  discussed  in  Chapter  III.  This  chapter 
discusses  each  step  in  detail. 

iOd  Qbfcking  B§§Q9D§§  iuCidSSS 

QsfiDina  tbf  Pcsbiss  selecting  Variables  of  Interest. 
Determining  which  variables  to  use  in  the  response  function  depends  on 
the  specific  purpose  for  which  the  response  surface  is  to  be  used. 
Because  applications  developed  are  geared  toward  determining  the  best 
economic  policy  for  the  federal  government,  the  exogenous  variables 
selected  for  manipulation  are  corporate  taxes  (Tc) ,  wage  taxes  (Tu), 
government  nonwage  spending  (G)  ,  government  wage  bill  (Wz) ,  and  number 
of  government  employees  (No).  These  variables  are  instruments  of  federal 
policy  which  are  broken  out  in  the  KG  model.  It  is  also  desirable  to  see 
how  the  economy  can  be  manipulated  through  monetary  policy  as  well. 


Unfortunately,  Soldberger  adnitted  that  Nhile  the  KG  nodel  did  an 
adequate  job  of  forecasting  current  eonetary  variables  (such  at  short 
tern  interest  rates,  long  tern  interest  rates  and  liquid  assets),  it  did 
a  poor  job  of  sinulating  econonic  changes  further  in  the  future 
(Goldberger,  1959:84-03).  Goldberger  clained  that  short  tern  interest 
rates  would  increase  without  bound.  In  addition,  the  nodel  failed  to 
capture  the  relationships  between  the  nodel  s  nonetary  sector  and  the 
other  sectors  adequately.  Consequently,  the  nonetary  sector  of  the  nodel 
IS  onitted  fron  further  consideration. 

There  are  several  possible  endogenous  variables  for  which  response 
surfaces  could  be  built.  Econonic  performance  indicators  conmonly  used 
to  asses  the  health  of  the  national  economy  include  percent  growth  in 
gross  national  product,  percent  unemployment,  and  percent  inflation. 
Other  measures  are  interest  rates  and  federal  deficits.  The  KG  model 
does  not  include  federal  deficits.  Because  of  limitations  mentioned  in 
the  preceding  paragraph,  the  model  cannot  be  used  to  study  the  dynamic 
properties  of  interest  rates.  In  this  study,  response  surfaces  are  built 
for  gross  national  product  (Q) ,  total  number  of  workers  employed  (Ni«) , 
and  the  price  index  <p). 

isLSCtinq  the  Resgonse  Surface  Eguatign.  Once  the  variables  of 
interest  have  been  selected,  the  next  step,  selection  of  the  response 
surface  equation  can  be  accomplished.  Because  of  the  arguments  enumer¬ 
ated  in  Chapter  111,  a  second  order  polynomial  has  been  selected  as  the 
response  surface  equation.  Since  Goldberger  argued  that  the  model  was 
nearly  linear  (Goldberger,  1959:136-138),  a  second  order  polynomial 
should  have  no  difficulty  approximating  model  ouputs.  In  addition,  if 
there  are  any  significant  interactions  between  variables  or  second  order 
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e'f'fects,  the  second  order  equation  Mill  pick  thee  up.  If  the  second 
order  polynoeial  does  not  fit  well  enough,  a  higher  order  polynoeial  can 
be  used. 

Because  one  of  the  objectives  of  the  study  is  to  reproduce 
nultipliers,  the  first  order  polynoaial  is  also  of  interest.  It  can  be 
shown  that  coefficients  coeputed  for  a  first  order  response  surface  for 
linear  econometric  models  are  mathematically  equivalent  to  the 
multipliers  derived  in  Chapter  II  (see  Appendix  6).  Because  the  KG 
model  is  nearly  linear,  the  coefficients  of  the  variables  in  the  first 
order  decoded  response  surface  should  compare  quite  closely  with 
multipliers  computed  by  Goldberger. 

itl^cting  the  Exger |.mentaL  Qesi,gn.  Selecting  the  variables  of 
interest  and  the  form  of  the  response  equation  narrows  the  choices  of 
experimental  designs.  A  three  level  full  factorial  design  with  five 
factors  would  require  3”  »  243  runs  of  the  KG  model  just  to  build 
response  functions  for  period  zero.  For  each  succeeding  period,  another 
set  of  243  runs  would  have  to  be  made.  A  more  economical  design  is  found 
in  Box  and  Behnken's  1960  paper  entitled,  "Some  New  Three  Level  Designs 
for  the  Study  of  Quantitative  Variables"  (Box  and  Behnken,  1960:460). 

Box  and  Behnken's  five  factor  design  has  only  46  design  points.  Of 
these,  five  are  redundant  center  points  which  can  be  eliminated  for  a 
deterministic  model.  The  design  is  highly  orthogonal  with  only  a  slight 
correlation  oetween  the  squared  and  intercept  terms.  The  Box  and 
Behnken  design  works  especially  well  for  variables  with  small  variance. 
Since  the  K-G  model  is  deterministic  without  error  terms,  the  output  has 
no  variance.  Consequently,  this  design  is  especially  appropriate.  The 


Box  and  Bahnkcn  dasign  is  aeant  for  use  with  a  second  degree  polynoaial, 
but  the  design  can  be  used  for  the  first  order  nodel  because  the  nain 
effects  are  uncorrelated.  Using  a  three  level  design  to  estieate  a  first 
order  model  has  the  added  advantage  of  enabling  assesseent  of  second 
order  effects  not  accounted  for  in  the  first  order  nodel.  If  other 
models  with  terms  higher  than  second  degree  are  used  as  response 
equations,  another  design  must  be  selected.  Appendix  D  shows  the  five 
factor  Box  and  Behnken  design  with  extra  center  points  deleted. 

QitSCliQLQa  5iQaS  at  iQaeaeQdeQt  The  next  step  in  RSM 

application  is  to  determine  the  range  of  the  independent  variables.  In 
this  study,  several  different  ranges  are  appropriate.  To  check  whether  a 
second  order  response  function  can  fit  the  model  output,  the  entire 
range  of  data  for  the  years  over  which  the  model  was  estimated  is 
appropriate.  For  those  exogenous  variables  in  the  response  surface 
equation,  the  maximum,  minimum,  and  average  of  maximum  and  minimum 
sample  values  comprise  the  three  levels  used  in  model  runs.  Variables 
not  in  the  response  equation  are  set  at  the  sample  mean.  The  objective 
of  these  runs  is  to  get  the  best  fit  possible. 

Goldberger's  multipliers  were  computed  at  the  sample  mean  for  unit 
changes  in  the  predetermined  variables.  For  runs  reproducing 
multipliers,  the  sample  mean  plus  or  minus  one  unit  are  used  as 
factor  levels. 

For  runs  used  to  evaluate  specific  policies,  response  surfaces 
should  be  constructed  for  exogenous  variable  levels  which  are  considered 
politically  feasible  by  the  decision  makers  using  the  analysis.  All 
predetermined  variables  not  included  in  the  response  surface  equation 
should  be  set  to  current  or  forecast  values  so  that  predicted  values  for 


tndogtnoua  varlablts  will  b*  raaaonablc. 

&S&9di!ia  t!lt  l>i98CLS8Q&iL  Q8iLftQ>  O^ct  thi  factor  level  a  have  been 
apeclfied,  the  coded  experieental  deaign  euat  be  decoded.  The  decoded 
deaign  apecifiea  the  actual  exogenoua  variable  valuea  for  which  a  aolu- 
tion  ia  to  be  coeputed.  To  reduce  drugery,  save  time,  and  decreaae 
arithmetic  errors,  a  FORTRAN  program  was  developed  to  automate  the 
task.  The  program  reads  a  coded  design  file  and  writes  a  decoded  version 
of  the  experimental  design  to  a  new  file.  The  actual  code  for  this 
program  is  in  Appendix  D. 

Q988l8tiQ9  Q988I  RuQI*  Creating  a  file  with  factor  levels 
specified  is  one  prerequisite  for  the  next  step  in  developing  a  response 
surface,  salving  the  model  for  each  design  point.  Also  needed  is  a  file 
specifying  the  values  of  the  other  predetermined  variables.  Once  these 
are  specified,  the  model  becomes  a  system  of  twenty-one  equations  in 
twenty-one  unknown  endogenous  variables.  Since  the  system  is  large 
and  nonlinear,  a  numerical  approach  to  solving  the  system  is  used. 
Because  of  the  arguments  set  forth  in  Chapter  III  regarding  the  best 
method  to  solve  a  maroeconometr ic  model,  the  Qauss-Sidel  method  is 
selected  as  the  algorithm  for  salving  the  K6  model.  Appendix  A  discusses 
the  mathematical  aspects  of  the  method  and  contains  the  FORTRAN 
implementation  of  the  method  as  applied  to  the  KG  model.  The  program 
reads  data  from  five  data  files  including  the  coded  design  file,  the 
decoded  design  file,  a  file  containing  predetermined  variable  values  not 
included  in  the  experimental  design,  an  initial  trial  solution  file,  and 
a  file  containing  control  language  for  the  program.  The  program  has  the 
capability  to  link  solutions  together  to  produce  extended  period 
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forecasts.  The  proqran  produces  a  data  file  Mhich  includes  the  coded 
experinental  design,  a  sequential  case  nunber  for  each  design  point,  the 
period  number  for  each  solution,  and  the  solution  for  all  twenty-one 
endogenous  variables.  This  output  file  eay  be  read  directly  by  the  BMDP 
statistical  package. 

EittiQS  Data.  Once  the  data  sets  are  generated,  response 
functions  are  fit  using  BMDP's  stepwise  regression  routine.  The  BhDP 
package  is  used  for  this  research  project  because  it  has  all  the 
capabilities  needed  for  fitting  response  surfaces,  is  familiar  to  the 
researcher,  and  is  available  to  the  researcher.  Stepwise  regression  is 
used  because  it  brings  in  independent  variables  one  at  a  time  in  order 
of  influence  on  the  dependent  variable. 

The  second  order  response  surface  equations  have  higher  order 
terms.  To  estimate  the  coefficients  for  these  terms,  appropriate 
transformations  are  made  in  BMDP's  control  language.  One  point  to 
reiterate  is  that  regressions  are  made  in  terms  of  the  coded  variables 
to  preserve  orthogonality.  The  coefficients  computed  by  BMDP  must 
therefore  be  decoded  (except  in  the  case  where  multipliers  are 
reproduced).  Before  decoding  the  coefficients,  however,  it  is  convenient 
to  check  the  fit  of  the  response  surface. 

Qtl§£kLQ9  ttlS  ELt-  The  four  methods  of  checking  the  response  func¬ 
tion  fit,  R‘,  S5E,  residuals,  and  random  points  are  all  useful  for  this 
study.  The  and  SSE  values  are  given  automatically  for  each  step  in 
the  stepwise  regression  procedure.  The  R^  value  is  the  primary  criteria 
for  deciding  which  variables  to  keep  in  the  response  surface  equation 
used  in  this  study.  However  the  SSE  is  checked  to  insure  that  it  is  not 


increasing  as  more  variables  are  brought  into  the  response  surface 


equation.  Variables  which  reduce  the  negligibly  (less  than  0.0001) 
can  be  omitted  from  the  response  surface  equation. 

While  the  value  gives  a  good  overall  average  measure  of  fit,  it 
is  also  useful  to  examine  the  percent  deviation  of  the  response  surface 
equation  from  the  K6  model  solution  for  each  design  point.  This  quantity 
is  computed  from  the  residuals  by  dividing  the  residual  by  the  actual 

response  value  and  multiplying  by  100  percent.  Of  interest  are  the 

largest  percent  deviations  and  where  those  deviations  occur. 

Finally,  the  most  stringent  test  of  response  surface  fit  is  the 
percent  deviation  between  the  model  responses  and  response  surface 
responses  for  random  points.  If  a  coded  random  experimental  design  with 
random  values  on  the  interval  (>1,1)  for  design  points  is  created,  the 
random  design  can  be  treated  just  like  a  regular  experimental  design. 
This  random  design  can  be  used  to  compute  new  data  points.  The  resulting 

data  file  can  then  be  appended  to  the  data  used  to  fit  the  response 

surface.  If  the  random  points  are  given  a  weight  of  zero  in  BUDP, 
residuals  are  computed  for  the  random  points,  but  the  points  are  not 
used  to  compute  response  surface  coefficients.  The  percent  error  is 
computed  from  the  random  point  residuals  to  provide  another  assessment 
of  response  function  fit. 

If  the  response  surface  fits  well,  (within  96  percent)  the 
coefficients  may  be  decoded  and  analysis  can  begin.  If  not,  a  new 
response  surface  function  and  experimental  design  must  be  selected  and 
the  fitting  procedure  repeated  until  a  satisfactory  fit  is  obtained. 
Decoding  the  response  surface  coefficients  is  simply  a  straight 
application  of  the  procedure  discusssed  in  Chapter  III. 


Once  response  sur'faces  have  been  iit  satisfactorily  and  decoded, 
the  response  surface  equations  aust  be  interpreted.  Sone  important 
questions  about  the  response  surface  equations  follow. 

1.  What  does  the  response  function  imply  about  the 
relative  contribution  of  each  input  and  the  relation¬ 
ship  between  inputs?  Do  these  implications  make  economic 
sense? 

2.  If  there  are  significant  higher  order  or  interaction 
terms,  why  do  these  occur?  Can  economic  theory  explain? 

3.  How  do  the  response  surface  coefficients  compare  with 
multipliers? 

4.  The  K-6  model  is  composed  of  mostly  linear  equations 
with  some  products  of  input  variables.  How  does 

this  affect  the  response  function? 

If  the  response  surface  equations  appear  valid,  applications  for 
the  equations  may  be  developed.  Chapter  III  suggests  several  uses  tor 
response  surfaces  fit  to  a  macroeconomic  model.  They  include  policy 
simulation,  trade  off  analysis,  and  optimization.  Describing  the  details 
of  these  applications  is  deferred  until  Chapter  VI. 

The  methodology  outlined  in  this  chapter  describes  what  steps  must 
be  taken  to  meet  the  research  objectives  set  forth  in  Chapter  I.  The 
Gauss-Sidel  numerical  method  of  solving  simultaneous  nonlinear  equations 
is  implemented  in  a  computer  program  to  solve  the  KG  model.  Second  order 
polynomial  response  surfaces  are  built  for  important  economic  indicators 
to  asses  how  well  response  surface  can  fit  the  macroeconomic  model.  A 
first  order  response  surface  is  estimated  and  the  coefficients  compared 
to  Goldberger's  multipliers.  Finally,  the  response  surfaces  are 


interpreted  and  applied. 


V.  Results 


Introduction 

Chapter  IV  described  what  data  are  required  to  achieve  the  research 
objectives  and  outlined  analysis  to  be  performed  with  the  data.  Computer 
runs  were  made  on  the  Air  Force  Institute  of  Technology's  VAX  11/780 
computer  to  collect  required  data.  Appendix  A  contains  the  actual 
FORTRAN  code  used  to  obtain  the  data.  This  chapter  summarizes  results  of 
the  regression  and  comments  on  significant  aspects  of  the  results.  It 
also  discusses  research  objectives  one  and  two  in  light  of  the  results, 
ifcond  Order  Model  Fit 

The  first  research  abjective  is  to  see  how  well  a  second  degree 
polvnomial  can  approximate  the  output  of  the  KG  model  when  five  factors 
are  changed  jointly.  To  satisfy  this  objective,  the  KG  model  is  solved 
for  period  :ero  and  period  five  at  factor  levels  required  by  the  Box  and 
Bennken  experimental  design.  (For  a  discussion  of  the  Gauss-Sidel 
numerical  technique  used  to  solve  the  Klein-Goldberger  model,  see 
Appendix  A.)  Second  order  polynomial  coefficients  are  estimated  for 
number  of  workers  employed  (Nw),  price  index  (p),  and  gross  national 
product  IQ)  in  terms  of  wage  taxes  (Tw),  corporate  taxes  (Tc),  govern¬ 
ment  nonwage  spending  (G)  ,  government  wage  bill  (Wz) ,  and  number  of 
government  workers  (Nnl  using  the  BMDP  2R  program,  stepwise  regression, 
'''he  general  form  of  each  equations  is 


Q  ®  Aoo  +  AoiTc+AosT  w'^dosG'^  Ao  ^  AosNb  laTcTw+disTcG  +  di^TcWa 

+  dl  ■5TcNs+d33Tw6  +  d24T  wW^  ^d2sTuNa^d34  GWz^d  SsGNG+d45W2NB+dl iTc^ 

+  d22Tw  ^  + A33  e=+a  44W2^+dssNB^ 

Nw  ®  boo  +boiTc^h  o2T  m  ‘^b  osG^b  o  4W2^bosNa^bi2TcTw^bi3TcG‘^bt4TcM2 

+bx5TcNB+b23TwG+b24  T^Wa+b  23  IvNo+b  34  GUa-t-b  35  GNa'^b  45  WaNQ+b»iTc= 

+  b22Tu‘‘  +  b33G^  +  b44W2^+b55No* 

p  *  C0O+C01TC+C02TU  +  C03G  +  C04W2  +  C05NS  +  C  12TcTw  +  Ci3TcG  +  Ci4T  cW2 

^-CisTcNg+C  23  T  gmG'^c  34  T  uWa  +  CasTwNG  +  C  3  4  GWa  +  C  3sGNo^C4sM2NQ^Ci iTc^ 

^C22Tw'‘  +  C33G'‘  +  C4  4  +  C  ssNs  “■ 


where  dij,  b,j,  and  Cij  are  the  coeffecients  to  be  determined.  The 
tallowing  conditions  are  applied  in  astimating  the  coefficients. 

1.  The  three  factor  levels  for  corporate  taxes  (Tc).  wage  taxes 
(Tu),  government  nonwage  spending  (G>,  government  wage  bill 
(Wa) ,  and  number  of  government  employees  (Nn)  are  the  maximum 
sample  value,  the  minimum  sample  value,  and  the  average  of  the 
maximum  and  minimum  samole  values. 

2.  The  Box  and  Behnken  five  factor  three  level  design  discussed  in 
the  methodology  chapter  with  redundant  center  points  deleted  is 
used. 

3.  All  other  predetermined  variables  are  set  at  sample  mean  values. 
In  computing  sample  means  for  lagged  variables,  the  appropriate 
data  values  from  the  periods  1923-1951  are  used,  (e.g.,  the 
sample  mean  for  the  price  index  lagged  one  year  includes  the 
price  indicies  for  1928  and  1944,  but  excludes  the  price 
indicies  from  1940  and  1952.) 

4.  For  each  design  point,  all  current  (nonlagged)  exogenous 
/enables  are  held  fixed  for  extended  period  forecasts  (beyond 
period  zero).  Lagged  variables  are  updated  with  new  values 
after  each  period's  forecasts  are  computed. 

5.  The  monetary  sector  is  suppressed  bv  excluding  the  liquidity 
forecasting  eouations,  Eqs  (2.2.12)  and  (2.2.13).  This  step  is 
necessary  to  match  Goldberger  s  analysis. 

6.  The  time  trend  variable  is  updated  bv  one  each  year. 

Tables  5.2a-f  summarize  the  results  of  stepwise  regression  for  each 


response  function.  They  include 

step  number,  entering  variable,  multiple 

R  and  R^,  and  change  in  R^. 

Table  5.1a. 

Sumnary  Table 

of  Stepwise  Regression  Results  for 

Nunber  of  Workers  Enployed  in  Period  Zero 

Step 

Variable 

Multiple 

Change 

No. 

Entered 

R 

R  = 

in  R^ 

1 

3  a 

.9091 

.8264 

.8264 

2 

5  No 

.  9681 

.9372 

.  1108 

3 

I  T  w 

.9843 

.9689 

.0317 

4 

4  W2 

.9993 

.9985 

.0296 

5 

2  TC 

1.0000 

1.0000 

.0015 

6 

41  6W2 

1. 0000 

1.0000 

.0000 

7 

31  G  = 

1.0000 

1.0000 

.0000 

a 

42  GNo 

1.0000 

1.0000 

.0000 

9 

35  TwG 

1.0000 

1.0000 

.  0000 

10 

36  TuM2 

1.0000 

1.0000 

.0000 

11 

38  TcG 

1.0000 

1.0000 

.0000 

Table  5.  It 

,  Summary  Table 

of  Stepwise  Regression  Results  for 

the  Price 

Index  in  Period  Zero 

•  w"* 

Step 

Variable 

Multiple 

Change 

w  * 

No. 

Entered 

R 

R* 

in  R* 

1 

3  G 

.6919 

.4787 

.4787 

2 

4  W? 

.9107 

.8294 

.3507 

3 

5  No 

.9887 

.9775 

.  1481 

v’!>' 

4 

1  T  w 

.  9979 

.9959 

.0184 

im 

5 

32  W,= 

.9991 

.9981 

.0023 

6 

2  Tc 

.9995 

.9990 

.0008 

7 

43  W,No 

.9997 

.9994 

.0004 

a 

31  Q  = 

.9999 

.9998 

.0004 

'-'.**** 

9 

41  GW, 

1.0000 

.9999 

.0002 

10 

35  TwG 

1.0000 

1 . 0000 

.  0000 

11 

33  No= 

1.0000 

1.0000 

.0000 

12 

42  GNo 

1.0000 

1.0000 

.0000 

13 

36  TwW, 

1.0000 

1.0000 

.0000 

14 

38  TcG 

1 . 0000 

1.0000 

.  0000 

Table  S.lc.  Suoifflary  Table  oi  Stepwise  Regression  Results  for 


Step 

No. 

Gross  National 

Variable 

Entered 

Product  in  Period  Zero 

Multiple 

R  R= 

Change 
in  R  = 

1 

3 

6 

.  9738 

.9484 

.9484 

2 

1 

Tu 

.9923 

.9847 

.0363 

3 

4 

Us 

.9992 

.9983 

.0136 

4 

2 

Tc 

1.0000 

1.0000 

.0017 

5 

5 

No 

1.0000 

1 . 0000 

.0000 

6 

41 

GW, 

1.0000 

1.0000 

.0000 

7 

31 

G  = 

1.0000 

1.0000 

.0000 

8 

42 

GNo 

1.0000 

1 . 0000 

.0000 

9 

35 

T„G 

1.0000 

1.0000 

.0000 

10 

36 

T  uWs 

1.0000 

1.0000 

.0000 

1 1 

38 

TcG 

1.0000 

1.0000 

.0000 

12 

37 

T  wNo 

1.0000 

1 . 0000 

.  0000 

13 

30 

Tc  = 

l.OOOO 

1 . 0000 

.0000 

5. Id.  Summary  Table 
Number  o-f  Workers 

Step  Variable 

No.  Entered 

of  Stepwise  Regression  Results  fi 
Employed  in  Period  Five 

Multiole  Change 

R  in  R^ 

1 

3 

G 

.  9286 

.8622 

.8622 

2 

2 

Tc 

.  9605 

.9226 

.0604 

3 

1 

T„ 

.  9903 

.9807 

.0501 

4 

5 

Ng 

.  9951 

.  9903 

.  0096 

5 

4 

w. 

.  9990 

.  9901 

.0078 

6 

31 

G  = 

.  9996 

.  9993 

.0012 

7 

35 

TwG 

.9997 

.9995 

.  0002 

8 

38 

TcG 

.  9998 

.  9997 

.  0002 

9 

41 

GW, 

.  9999 

.9998 

.0001 

10 

42 

GNc 

1 . 0000 

.  9999 

.0001 

11 

36 

T  wWs 

1.0000 

.  9999 

.  0000 

12 

39 

TcW^ 

1.0000 

.  9999 

.0000 

13 

34 

T  uiT  c 

1 . 0000 

l.OOOO 

.0000 

14 

40 

T  cNg 

1 . 0000 

1 . 0000 

.  0000 

'.V.*. 


'  ^ 
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Table  S.le.  Suninary  Table  of  Stepwise  Regression  Results  for 
the  Price  Index  in  Period  Five 


Step  Variable  Multiple  Change 


No. 

Entered 

R 

R= 

in  R= 

1 

3 

6 

.8852 

.7836 

.7836 

2 

4 

W:: 

.9268 

.8590 

.0754 

3 

S 

Ng 

.9554 

.9128 

.0538 

4 

1 

Tm 

.9806 

.9617 

.0489 

5 

2 

Tc 

.9991 

.9982 

.0365 

6 

32 

W3  = 

.9993 

.9986 

.0004 

7 

43 

UzNq 

.  9995 

.9990 

.0004 

8 

41 

GW:, 

.  9997 

.9994 

.0004 

9 

42 

GNo 

.9998 

.  9996 

.0002 

10 

33 

Nq  = 

.  9999 

.9997 

.0001 

11 

38 

TcG 

.9999 

.  9998 

.0001 

Table  S.lf.  Summary  Table  of  Stepwise  Regression  Results  for 
Gross  National  Product  in  Period  Five. 

Step  Variable  Multiple  Change 


Na. 

Entered 

R 

R  = 

in  R  = 

1 

3 

G 

.  9305 

.8658 

.8658 

2 

2 

Tc 

.  9701 

.9411 

.0753 

3 

1 

Tw 

.  9990 

.9980 

.0569 

4 

31 

G= 

.  9995 

.9990 

.0010 

5 

5 

Ng 

.  9997 

.9993 

.0003 

6 

35 

TwG 

.9997 

.9995 

.0002 

7 

38 

T  cG 

.9998 

.  9996 

.0001 

a 

41 

GW,. 

.  9999 

.  9998 

.0002 

9 

4 

Wa 

.  9999 

.9999 

.  0001 

10 

42 

GNo 

1.0000 

.  9999 

.0000 

1 1 

36 

T  wWa 

1.0000 

.  9999 

.0000 

12 

39 

T  cWa 

1.0000 

1 . 0000 

.0000 

13 

34 

T  wT  c 

1.0000 

1.0000 

.  0000 

14 

40 

T  cNo 

l.OOOO 

1.0000 

.  0000 

15 

37 

T  wNg 

1.0000 

1 . 0000 

.0000 

16 

43 

W^No 

1 , 0000 

1.0000 

.  0000 

The  data  in  the  Tables  5.la-f  yield  two  important  conclusions. 
First,  for  both  period  zero  and  period  five,  the  response  surfaces  tit 
the  data  well  as  reflected  by  the  R*  column.  In  no  case  does  the 
exceed  0.9983.  A  second  check  of  fit  is  the  percent  error  for  the  design 
points.  Table  S.2  lists  the  design  points  with  the  largest  percent  error 


ior  each  response  sur'face  as  computed  from  the  residuals 


Table  3.2.  Design  Point  Pit  Check  for 
Second  Order  Response  Surface. 


Resid¬ 

ual 

Predict 

Value 

Pet 

Error 

Tw 

Factor  Levels 

Tc  S 

M:! 

Ng 

Period 

Zero 

Nw 

-0.0041 

42.69 

0.0097 

1.0000 

-l.OOOO  0.0000 

0.0000 

0.0000 

P 

0.  1631 

95.00 

0.1717 

0.0000 

0.0000  -1.0000 

1.0000 

0.0000 

Q 

-0.0059 

88.30 

0.0067 

1.0000 

0.0000  -1.0000 

0.0000 

0.0000 

Period 

Five 

Mw 

-0.  1745 

14.06 

1.2567 

1.0000 

0.0000  -1.0000 

0.0000 

0.0000 

P 

3.  1580 

35.53 

8. 1627 

0.0000 

0.0000  -l.OOOO 

l.OOOO 

0.0000 

3 

-0.4030 

46.91 

0.8140 

l.OOOO 

0.0000  -l.OOOO 

0.0000 

0.0000 

The  data  in  Table  3.2  imply  good  fit.  Mith  the  exception  of  price 
index  in  period  five,  design  point  error  is  less  than  1.3  percent.  Price 
index  in  period  five  has  a  larger  percent  error  of  8.16  percent.  Fifty 
random  points  in  the  operating  region  were  also  run  to  check  fit.  The 
points  with  the  largest  percent  error  for  each  surface  are  shown  in 
Table  5.3.  The  lack  of  fit  is  extremely  pronounced  for  period  five  price 
index  (26.22  percent ) . 


59 


Table  S.3.  Random  Point  Fit  Check  for 
Second  Order  Response  Surface 


Resid¬ 

ual 

Predict 

Value 

Pet 

Error 

T„ 

Factor  Levels 

Tc  G 

U: 

Ng 

Period 

Zero 

Nu 

-0.0053 

19.38 

0.0116 

0.9154 

-0.8486  -0.2112 

-0.8846 

0.8269 

p 

-0.2690 

111.7 

0.2291 

0. 1516 

0.9346  -0.8707 

-0.8653 

-0.8905 

S 

0. 1082 

91.54 

0.0118 

0. 1516 

0.9346  -0.8707 

-0.8653 

-0.8905 

Period 

F  i  ve 

Nw 

0.8136 

19.38 

4.0290 

0.8024 

-0.7953  -0.9008 

0.9026 

-0.7324 

P 

6. 2620 

17.62 

26.2206 

0.2778 

0.1031  -0.9972 

0.8670 

-0.4157 

0 

2.0090 

74.96 

2.6101 

0.8024 

-0.7953  -0.9008 

0.9026 

-0.7324 

The  large  percent  errors  are  evidence  that  the  resoonse  surface 
does  not  fit  the  price  index  response  for  the  entire  operating  region 
aoecified.  Of  all  the  response  surfaces,  one  might  expect  the  response 
surface  for  p  to  be  the  most  difficult  to  fit.  Of  the  six  nonlinear 
equations  in  the  KG  model,  p  appears  multiplied  with  other  endogenous 
variables  in  Egs  (2.1.9),  (2.1.10),  and  (2.1.20).  If  one  were  to  solve 
for  0  in  terms  of  the  other  endogenous  variables,  an  endogenous  variable 
would  be  in  the  denominator.  Perhaps  a  higher  order  polvnomial  or 
logrithmic  function  can  provide  a  closer  approximation  for  the  p 
resoonse  surface.  Although  the  resoonse  surface  theoretically  should  fit 
the  response  throughout  the  whole  region,  it  is  interesting  to  note  that 
large  errors  for  p  occured  at  small  values  of  p.  The  design  points  and 
the  random  points  with  the  largest  percent  error  also  hao  the  smallest 
/alues  for  p.  In  fact,  for  the  point  with  the  largest  percent  error,  the 


value  for  p  was  23.88,  but  the  seallest  saople  value  for  p  fron  the 
1929-1952  data  was  90.7.  In  addition,  other  values  of  exogenous 
variables  were  extreecly  far  reeoved  fron  the  sample  data  for  this  case. 
It  seees  unlikely  that  real  world  analysis  would  be  conducted  in  this 
region  of  the  response  surface.  The  largest  percent  error  for  any  design 
or  random  point  with  a  p  value  over  90.0  was  0.83  percent  indicating  a 
good  fit  in  the  range  of  real  world  response. 

Apparently,  the  output  of  the  KG  model  can  indeed  be  approximated 
by  a  low  order  polynomial.  To  be  absolutely  certain  on  this  point, 
response  surfaces  would  have  to  be  built  for  all  endogenous  variables 
which  included  all  predetermined  variables  for  all  periods.  For  practi¬ 
cal  applications,  however,  if  closely  fitting  response  surfaces  can  be 
built  for  the  endogenous  variables  of  interest  which  include  the  prede¬ 
termined  variables  of  interest  and  which  cover  the  time  frame  of  inter¬ 
est,  this  is  all  that  is  necessary  to  proceed  with  analysis.  Further¬ 
more,  there  is  no  reason  to  believe  that  other  closely  fitting  response 
functions  cannot  be  developed  for  any  endogenous  variables  in  terms  of 
any  predetermined  variables. 

The  second  major  conclusion  to  be  drawn  from  Tables  5.la-f  is  that 
first  order  terms  account  for  most  of  the  variation  in  the  data.  Table 
5.4  lists  the  percent  of  variation  explained  by  first  order  terms  for 
each  resoonse  surface  equation.  These  values  were  obtained  by  fitting  a 
first  order  model  to  a  second  order  Box  and  Behnken  experimental  design 
(Box  and  Behnken,  1960:460) 


Table  3.4.  Values  for  The  First  Order  Response 
Surface  Epuation. 

Period  Period 

Zero _ Five 

Nw  1.0000  0.99B1 

p  0.9967  0.9982 

Q  1.0000  0.9985 

The  high  R-^  values  indicate  that  the  nodel  is  very  nearly  linear. 
One  eight  expect  price  index  to  show  significant  nonlinearity.  Of  the 
six  equations  in  the  K6  eodel  which  contain  nonlinearities,  price  is 
involved  in  three.  Although  the  response  function  for  price  index  in 
periods  zero  and  five  are  eore  nonlinear  than  either  number  of  employed 
workers  or  gross  national  product,  99.7  percent  of  the  variation  is 
explained  by  linear  terms.  These  observations  are  consistent  with  Gold- 
berger  s  argument  that  the  model  is  very  nearly  linear  and  that 
multipliers  that  he  computed  at  the  sample  mean  are  accurate  for  a 
large  range  of  time  series  data. 

MyLlLBlitCS  Response  Surface  Coefficients  CgmBared 

The  second  research  objective  is  to  reproduce  Goldberger’s  multi¬ 
pliers  using  a  first  order  response  surface  equation.  To  do  this,  runs 
are  made  with  the  following  conditions: 

1.  The  three  factor  levels  of  corporate  taxes,  wage  taxes, 
government  nonwage  spending,  government  wage  bill,  and  number 
of  government  employees  are  the  sample  means  plus  or  minus  one 
unit. 

2.  The  experimental  design  used  is  the  Box  and  Behnken  five 
factor  three  level  design  discussed  in  the  methodology  chapter 
with  redundant  center  ooints  deleted. 


3.  All  other  predetermined  variables  are  set  at  sample  mean 
values.  For  lagged  variables,  the  appropriate  data  values  from 
before  the  periods  1929-1940,  and  1945-1952  are  used,  (e.g.,  t 
sample  mean  for  the  price  index  lagged  one  year  includes  the 
price  indicies  for  1928  and  1944  but  excludes  the  price 
indicies  from  1940  and  1952.) 

4.  For  each  design  point,  all  current  (nonlagged)  variables  are 
held  fixed  for  extended  period  forecasts  (beyond  period  :ero). 
Lagged  variables  are  updated  with  neu  values  after  each 
period's  forecasts  are  computed. 

5.  The  monetary  sector  is  suppressed  by  excluding  the  liquidity 
forecasting  equations.  This  step  was  necessary  to  match 
Soldberger's  analysis. 

6.  The  time  trend  variable  is  not  updated  since  Goldberger 
computes  a  separate  multiplier  to  account  for  the  time  trend. 

The  conditions  were  applied  to  correspond  to  the  assumptions  made 
by  Goldberger  in  developing  his  multipliers.  Tables  5.5a  and  5.5b  sum¬ 
marize  the  results,  They  show  multipliers  computed  by  Goldberger  for  a 
unit  increase  in  government  spending  and  the  corresponding  response 
surface  equation  coefficients.  Multipliers  are  taken  from  Table  5.2  of 
Imgact  MuLtLalitrs  and  Dynamic  Progerties  of  the  Kl_ein;Ggl,dberg|r  Model 
(Goldberger,  1959:80). 


Table  5.5a.  Multipliers  for  Unit  Increase  in  Government  Spending, 


Period 
2  3 


1.628  1.842 
4.631  5.911 
3.884  4.565 
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VI.  BSSBSDSS  iyCi^ce  Analysis 


Introduction 

Chapter  V  demonstrates  that  a  response  surface  can  indeed  fit  the 
output  of  the  KG  model  with  a  Iom  order  polynomial  with  a  high  degree  of 
accuracy.  Further,  RSM  verifies  that  the  model  is  very  nearly  linear 
with  first  order  terms  accounting  for  over  99  percent  of  the  total 
variation  for  all  response  surfaces  generated.  Since  these  surfaces  fit 
so  well,  one  may  conclude  that  they  are  accurate  representations  of  the 
model's  characteristics  and  may  be  used  as  an  approximation  to  the  model 
for  analysis.  This  chapter  examines  what  the  response  surfaces  mean  and 
sxolores  some  analysis  possibilities  emphasising  practical  applications. 
B§SG9D§§  §yCi§Sf 

The  response  surfaces  generated  for  the  Klein-6ol dberger  Model 
summarice  relationships  in  the  model  presenting  the  impacts  of  predeter¬ 
mined  variable  changes  on  current  endogenous  variables  explicitly.  An 
example  serves  to  illustrate.  For  period  five,  the  decoded  response 
surface  equation  for  gross  national  product  (Q)  in  terms  of  wage  taxes 
(Tu),  corporate  taxes  (Tc),  government  nonwage  spending  (G),  government 
wage  bill  (W^),  and  number  of  government  employees  (No)  is 

Qs  =  54.4254  -  4.0951T«  -  4.3058Tn  ♦  5.0849G  -  0.5864W=  -  0.4145Ng 

-  0.01346=  ♦  0.0244T„G  +  0.0226TcG  ♦  O.OlTSGWj  (6.1) 

The  coefficients  in  Eq  (6.1)  are  computed  by  stepwise  regression. 

It  is  assumed  that  the  accuracy  afforded  by  Eq  (6.1)  is  sufficient  for 
purposes  of  discussion.  The  R=  value  for  this  response  surface  is 
0.9998. 


L-; 
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Uhat  is  the  signi'f icance  Eq  (6.1)?  First,  the  equation  gives 
Mhat  the  KG  inadel  Qredic^on  will  be  far  any  coebination  of  T..,  Tc,  6, 
M2,  and  No  in  the  operating  range  of  interest.  The  equation  does  not 
clain  to  make  any  predictions  about  the  econoey  that  the  original  eodel 
could  not  nake.  The  response  surface  is  only  as  good  as  the  underlying 
model.  The  response  surface  also  only  purports  to  characterize  the  gross 
national  product  in  terms  of  the  predetermined  variables  in  Eq  (6.1) 
with  other  predetermined  variables  at  their  sample  means.  There  very 
well  could  be  interactions  between  variables  not  included  in  Eq  (6.1) 
(e.g.,  gross  national  product,  investment,  and  prices  from  the  previous 
period)  and  the  variables  appearing  in  Eq  (6.1)  (i.e.,  T..,  Tc,  6,  N2, 
and  No),  For  analysis  using  a  response  surface,  the  predetermined 
variables  not  included  in  the  response  surface  equation  should  be  set  at 
values  close  to  what  they  would  be  for  the  particular  economic  simula¬ 
tion  under  study.  For  instance,  if  a  study  is  to  be  made  of  the  effects 
of  government  spending  and  taxes  on  gross  national  product  two  years  in 
the  future,  then  the  macroeconomic  model  used  to  generate  the  response 
surface  should  have  lagged  endogenous  and  exogenous  variables  set  at 
their  appropriate  current  levels  or  what  they  are  expected  to  be.  Mhat 
Eq  (6.1)  does  give  is  the  relationship  between  Tu,  Tc,  6,  M2,  and  Ng, 
and  Q  in  the  Klein-Goldberger  model  in  the  operating  region  of  interest 
for  a  period  five  years  in  the  future  with  all  other  oredetermi ned 
variables  at  samole  means. 

Eq  (6.1)  contains  interaction  and  squared  terms.  These  terms 
suggest  that  the  change  in  Q  due  to  a  change  in  a  particular  predeter¬ 
mined  variable  is  dependent  on  its  own  or  another  variable's  level.  For 
example,  one  may  want  to  know  the  effect  of  increasing  corporate  taxes 


on  grois  national  product  Hve  years  in  the  future.  Taking  the  first 
partial  derivative  of  Eq  (6.1)  with  respect  to  Tc  yields  an  expression 
relating  the  change  in  Q  to  a  change  in  Tc. 

ao  =  -4.3058  ♦  0.02266  (6.2 

STc 

The  right  hand  side  of  Eq  (6.2)  is  a  nonconstant  "multiplier".  Eq 
(6.2)  suggests  that  the  change  in  Q  due  to  a  change  in  Tc  is  dependent 
on  the  level  of  6  as  shown  in  Figure  6.1.  It  is  important  to  note  that 
Eq  (6.2)  is  valid  only  for  the  ranges  of  Te  and  6  used  to  build  the 
response  surface  (Tc  ranges  from  $  0.40  to  $  11.88  billion  and  6  ranges 
from  t  11.5  billion  to  $  41.7  billion). 


Figure  6.1.  Dependence  of  dQ/aTc  on  G  in  Five  Years 
It  IS  of  interest  to  determine  why  the  change  in  Q  due  to  a  change 


in  Tr.  should  depend  on  G.  Corporate  taxes  directly  effect  gross 


national  product  through  confunption  and  invcstnint  in  Eqs  (2.1.1), 
(2.1.2),  (2.1.19),  and  (2.1.21).  However,  corporate  taxes  also  af'fect 
another  elenent  o^  gross  national  product  in  Eq  (2.1.9),  foreign  in¬ 
ports.  It  is  inportant  to  note  that  it  is  not  corporate  taxes  alone 
that  affect  foreign  investnent,  but  the  product  of  the  price  index  and 
corporate  taxes.  Economic  theory  asserts  that  government  spending  has  a 
strong  effect  on  prices.  Consequently,  there  is  an  interaction  between 
corporate  taxes  and  government  spending  in  determining  foreign  imports 
and  hence  gross  national  product.  Eq  (2.1.10),  which  models  the  deter¬ 
minants  of  farm  income,  also  has  a  similar  interaction  between  corporate 
taxes  and  prices.  One  way  to  visualize  the  magnitude  of  the  Tc6  interac¬ 
tion  term  is  to  plot  Q  versus  Tc  at  different  levels  of  G.  Figure  6.2 
shows  Q,  at  the  five  year  point,  as  a  function  of  Tc  for  three  levels  of 
G.  The  TcG  term  causes  a  change  in  slope  at  different  G  levels.  The 
change  in  slope  is  barely  discernable. 


Figure  6.2.  Relationshio  Between  Qr.  and  Tc  at 
Different  Levels  of  G  in  Period  Five. 


To  generate  Figure  6.2,  all  predeternined  variables  except  Tc  and  6 
are  set  at  sample  eeans. 
j  QoSBUtiQS  t^yltiBlisrs 

Chapter  Five  shows  the  close  correspondence  between  multipliers 
and  decoded  first  order  response  surface  coefficients.  In  fact,  Appendix 
I  B  shows  that  they  are  equivalent  for  linear  systems.  Since  both  RSM 

and  the  derivative  method  yield  essentially  the  same  numerical  values 
for  multipliers,  Goldberger's  extensive  analysis  applies  to  RSn  derived 
multipliers  as  well. 

I 

If  the  K6  model  were  more  nonlinear,  interaction  and  squared  terms 
would  became  more  significant.  It  is  here  that  response  surface  methodo¬ 
logy  provides  an  advantage  over  multiplier  analysis.  By  using  response 
surfaces,  one  can  detect  interactions  between  predetermined  variables  as 
noted  in  the  last  section.  To  generate  "multipliers"  from  response 
surface  equations,  one  computes  the  partial  derivative  of  the  response 

i 

surface  equation  with  respect  to  the  variable  of  interest.  The  last 
section  comouted  a  multiplier  for  changes  in  Q  due  to  changes  in  Tc. 

This  multiplier  together  with  other  multipliers  computed  from  Eq  (6.1) 
are  listed  below. 

^Q_  =  -4.3058  +  0.0226G  (6.2) 

Sfe 

ao  =  -4.0951  ♦  0.0224G  16.3) 

StL 

9  Q  =  5.0894  +  0.0244Tw  ♦  0.0226Tr.  -  0.0268G  +  0.0175Wr  (6.4) 

d  G 


SQ  =  -0.5847  ♦  0.0175G 


(6.5) 


so  =  -0.4145  (6.6) 

SNg 

Eqs  (6.2)  through  (6.6)  give  "function'*  nultipliers  which  capture  the 
relationship  between  0  and  T..,  Te,  &i  Us,  and  No  more  accurately  than 
traditional  multipliers. 

Eolith  iL™uLiti.Qn 

Using  the  multipliers  computed  above,  economists  can  answer  "what 
if"  questions  easily.  For  Instance,  if  an  economist  wants  to  know  the 
impact  on  gross  national  product  in  five  years  of  increasing  government 
spending  by  five  billion  dollars  and  paying  for  it  with  a  five  billion 
dollar  increase  in  wage  taxes,  he  can  use  the  multipliers  to  forecast 
the  answer.  Assuming,  for  illustration  purposes,  that 

Tu  =  $8  billion 
Tc  =  #10  billion 
G  =  $40  billion 
Wt  =  $16  billion 

and  all  other  predetermined  variables  are  at  sample  means,  then  the 
multiolier  relating  changes  in  Tw  to  changes  in  Q  is 

-4.0951  +  0.0244(40)  =  -3.1191 

from  Ea  (6.3).  The  multiplier  relating  changes  in  G  to  changes  in  Q  is 

5.0894  ♦  0.0244(8)  +  0.0226(10)  -  0.0268(40)  +  0.0175(16)  =  4.7186 

from  Eq  (q.4).  The  assumed  values  for  Tu,  Tc,  G,  and  W?  are  close  to  1952  samole 
values  from  the  data  used  to  estimate  the  model  (Klein  and  Goldberger, 


1955:131-132)  and  all  other  predetermined  values  are  at  sample  means.  A 
Hve  billion  dollar  increase  in  wage  taxes  changes  gross  national 
product  by  (15  bi 1 1 i on)  (-3. 1 191 )  »  -115.60  billion.  A  five  billion 
dollar  increase  in  government  spending  increases  gross  national  product 
by  ($5  bi 1 1  ion) (4. 7186)  s  $23.59  billion.  The  net  change  is  $23.59 
billion  -  $15.60  billion  =  $7.99  billion. 

It  IS  interesting  to  compare  the  multiplier  computed  above  to 
Goldberger's  multipliers  and  the  corresoonding  first  order  response 
surface  coefficients  (see  Table  5.5).  As  an  example,  Table  6.1  compares 
the  three  types  of  multipliers  for  changes  Q  due  to  changes  in  G  in 
oeriod  5. 

Table  6.1.  Multiplier  Comparisons  for  Changes  in  Q 
Due  to  Changes  in  Q  in  Period  5. 

Goldberger's  Multiplier . ...4.922 

First  Order  Response  Surface  Multiplier . 4,766 

Second  Order  Response  Surface  Mul t i pi i er . , . . 4, 716 

The  difference  shown  between  Goldberger's  multiplier  and  the  first  order 
response  surface  coefficient  is  the  accumulated  error  from  the  way  in 
wnich  Goldberger  linearized  the  model.  The  difference  between  the  first 
order  response  surface  coefficient  and  the  second  order  multiplier  is 
that  the  first  order  multiplier  is  computed  at  the  sample  mean  but,  the 
second  order  multiplier  is  computed  at  values  given  in  the  example 
above.  If  the  other  variables  in  Ea  (6.1)  had  been  at  different  levels, 
the  second  order  response  surface  multiplier  would  also  have  been 
different.  If  the  KG  model  were  more  nonlinear  the  difference  would  have 
oeen  more  pronounced.  Figure  6.3  shows  the  difference  between  the 
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multipliers  graphically.  Depicted  are  multipliers  for  the  full  range  of 
6  with  other  variables  fixed  at  the  levels  specified  above. 


Figure  6.3.  A  Graphical  Comparison  of  Multipliers  for  Changes  in  Q 
Due  to  Changes  in  G  in  Period  Five. 

By  using  response  surfaces  and  multipliers  generated  from  multi¬ 
pliers,  economists  can  answer  many  questions  without  repeated  runs  of 
the  macroeconomic  model.  Furthermore,  interaction  and  squared  terms  are 
identified  with  response  surfaces,  but  not  with  traditional  multiplier 
analysis.  There  are  still  other  valuable  uses  for  response  surfaces. 
l!DB9C5§D£f 

Response  surface  equations  can  be  used  to  evaluate  factor  imqor- 
tance  in  determining  the  response  variable  value.  P.W.  Smith  and  J.  M. 
heiiicnamp  show  how  to  evaluate  factor  importance  for  a  nuclear  exchange 
model  in  their  paper  entitled  "A  Methodology  for  Multidimensional  Impact 
Analysis  for  Military  Problems"  (Smith  and  Mellichamp,  1979).  In  this 
paoer  the  authors  ooint  out  that  the  size  of  the  factor  coefficient 


gives  the  relative  inpact  on  response  per  unit  of  factor.  The  factor 
with  the  largest  coefficient  has  the  nost  influence  per  unit  of  factor. 

The  relative  nagnitude  of  coefficients  in  the  response  surfaces 
generated  for  the  K6  model  do  give  some  neasure  of  the  influence  of 
factors  which  are  measured  in  the  sane  units.  For  instance,  the  first 
order  coefficients  presented  in  Appendix  E  indicate  that  for  each  dollar 
increase  in  government  spending,  gross  national  product  increases  by 
4.5200,  but  for  each  dollar  decrease  in  wage  taxes,  gross  national 
product  increases  by  3.1234  in  period  five.  Qoldberger  pointed  out  that 
determining  which  factors  are  most  influential  in  causing  endogenous 
variable  changes  from  an  historical  point  of  view  also  involves  the 
amount  by  which  the  factor  changes  from  period  to  period.  Two  factors 
with  the  sane  response  surface  coefficient  or  multiplier  do  not  have  the 
same  influence  on  an  endogenous  variable  if  one  changes  by  only  a  small 
increment  and  the  other  changes  by  a  large  increment.  To  measure  the 
relative  importance  of  predetermined  variables  in  determining  endogenous 
variable  values  Qoldberger  formulated  an  index  which  was  equal  to  the 
appropriate  multiplier  multiplied  by  the  sum  of  the  absolute  values  of 
the  changes  from  one  period  to  the  next  during  the  sample  period  and 
divided  by  the  number  of  periods  (Qoldberger,  1959:72-73).  Computing  an 
eauivaient  index  with  response  surface  coefficients  is  certainly 
possible.  Such  an  influence  index  is  useful  in  quantifying  the  histori¬ 
cal  impact  of  predetermined  variables  on  endogenous  variables.  However, 
from  a  policy  simulation  point  of  view,  another  measure  might  provide 
more  useful  information. 

If  a  policy  maker  has  influence  to  change  economic  oolicv  variables 
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over  a  limited,  politically  feasible  range,  then  the  policy  maker  would 
be  interested  in  which  policy  variable  at  his  disposal  would  be  most 
influential  in  bringing  about  desired  objectives.  If  the  policy  maker 
built  a  response  surface  using  the  maximum  and  minimum  politically 
feasible  values  as  factor  levels  in  the  experimental  design,  then  the 
coefficient  gives  another  measure  of  influence  of  that  policy 
variable.  For  example,  if  a  policy  maker  feels  the  maximum  government 
expenditures  that  Congress  will  approve  is  $200  billion,  while  the 
minimum  is  $170  billion,  then  he  could  build  a  response  surface  using  an 
experimental  design  with  $200  billion  and  $170  billion  as  factor  levels 
for  running  design  points  through  his  econometric  model.  Factor  levels 
for  other  policy  variables  of  interest  would  be  formulated  in  the  same 
wav.  The  resulting  coded  coefficients  give  the  amount  of  change  that 
could  be  brought  about  by  varying  the  policy  variable  over  its 
politically  feasible  range. 

QESi!!!i?iii9D  ^EEii£iii9D5 

The  explicit  form  of  the  response  surface  equation  with  the  unknown 
endogenous  variable  on  one  side  of  the  equation  and  known  predetermined 
variables  on  the  other  side  of  the  equation  suggests  further  aqplica- 
tions.  Because  response  surface  equations  have  the  form  that  they  do  and 
are  expressed  in  terms  of  actual  levels  instead  of  changes  (as  in  Gold- 
berger's  linearized  K(3  model)  economic  optimization  problems  can  be 
easily  formulated  and  solved.  An  example  serves  to  illustrate. 

Suppose  the  year  is  1952.  The  klein-Goldberger  model  has  just  been 
estimated  and  an  elected  oolicy  maker  wishes  know  what  combination  of 
fiscal  policies  will  maximize  economic  growth  (GNP)  ,  while  holding 
inflation  and  employment  at  or  below  acceqtable  levels.  The  official 
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Mould  likt  these  conditions  to  be  realized  in  about  three  years.  The  KG 
RSH  model  can  provide  some  guidance. 

To  solve  this  problee,  several  response  surfaces  nust  be  generated. 
As  shown  in  Chapter  II,  the  KG  nodel  can  be  linked  together  to  obtain 
forecasts  for  several  periods  in  the  future  by  solving  the  eodel , 
setting  lagged  variables  equal  to  the  current  variable  values,  and  then 
resolving  the  model.  It  is  assumed  that  changes  in  fiscal  policy 
variables  made  in  period  zero  are  sustained  until  period  three.  Salving 
the  problem  requires  construction  of  three  response  surfaces,  one  for 
gross  national  product  (Q),  one  for  price  level  (p),  and  one  for  number 
of  workers  employed  (Nw)  for  a  time  period  three  years  in  the  future. 

The  fiscal  policy  variables  available  for  manipulation  are  government 
nonwage  spending  (6),  wage  taxes  (Tu),  corporate  taxes  (Tc),  and 
government  wage  bill  (W:). 

To  generate  the  response  surfaces  needed,  all  predetermined 
variables  except  the  four  fiscal  policy  variables  are  set  at  expected 
constant  levels,  then  the  policy  variables  are  set  at  the  levels 
required  by  an  experimental  design,  and  the  forecasts  are  computed.  From 
the  resulting  data,  stepwise  linear  regression  is  used  to  estimate  the 
coefficients  of  first  order  response  equations  for  each  economic  per¬ 
formance  indicator.  Shown  below  are  the  equations  generated  from  the  KG 
model  with  the  predetermined  variables  set  at  selected  levels,  based  on 
1952  data  which  were  the  most  current  data  used  to  estimate  the  model 
(Klein  and  Goldberger,  1955:131-133). 

Q  =  51.4040  -  3.1087Tw  -  3.3323Tc  +  4.5025G  -  0.2854W,  (6.7) 

NW  =  13.1882  -  1.2541T„  -  1.2567Tr.  +  1.8021G  -  0.50208W:.  16.8) 


1 


p  =  143.6616  -  4.6051T«  -  4.1155Tc  +  6.6816G  -  5.6552W2  (6.9) 

Since  the  quantitative  relationships  in  the  KG  model  are  approxi¬ 
mated  quite  adequately  by  linear  functions  (See  Table  S.4.),  a  two 
level  four  factor  factorial  design  was  used  to  estimate  the  coefficients 
in  Eq  (6.7),  (6.8),  and  (6.9).  Design  variables  were  varied  over  a 
limited  politically  feasible  range  (IS. 63-11. 63  billion  for  Tu,  17.14- 
13.14  billion  for  Tc,  and  137.7-61.7  billion  for  G,  and  *13.82  to  121.82 
billion  for  Ws.)  These  ranges  were  set  by  looking  at  the  historical 
record  of  change  over  the  sample  period  and  then  making  a  reasonable 
guess  as  to  possible  ranges  of  change. 

From  the  feasible  ranges  and  response  surface  equations  above,  one 


can  formulate  a  linear  programming  problem  as  follows, 
hax i mi ze 

a  =  51.4040  -  3.1087T«  -  3.3323Tc  +  4.5025G  -  0.2B54W2  (6.10) 

Subject  to 

NW  =  13.1882  -  1.2541Tw  -  1.2567Tc 

+  1.8021G  -  0.50198W2  =  58.71  (6.11) 

p  =  143.6616  -  4.6051Tu  -  4.1155Tc 

♦  6.6016Q  -  5.6552W2  i  207.714(1.05)=^  (6.12) 

Tw  <  11.63  (6. 13) 

T„  5.63  (6.  14) 

Tr.  <  13.  14  (6.  15) 

Tc  >  7. 14  (6.16) 

G  <  61.7  (6.17) 


G  ;■  37.7 


(6. 18) 


Us  <  21.82 


(6. 19) 


Us  >  13.82 


(8.20) 


The  abjective  function,  Eq  (6.10),  is  sinply  the  response  surface 


for  Q.  The  first  constraint,  Eq  (6.11),  is  derived  as  follows.  1952  data 


indicate  the  number  of  workers  in  the  labor  force  (N,.)  is  66.6  million, 


the  number  of  workers  employed  (Nw)  is  56.0, the  number  of  self  employed 


workers  (Ne)  is  6.3  million,  and  the  number  of  farm  workers  (Nf)  is 


4.0  million.  Klein  and  Goldberqer  define  the  number  of  unemployed 


persons  (Nu)  to  be  (Klein  and  Goldberger,  1955:19): 


Nu  *  Ni_  *  (.Nw  +  Ne  ^  Nf  1 


For  the  1952  data  Nu  is  0.3  million  workers.  This  translates  to  an 


unemployment  rate  of  0.45  percent  (this  figure  is  clearly  unrealistic). 


The  number  of  self  employed  and  farm  workers  together  have  been 


decreasing  by  about  one  percent  per  year,  and  the  total  labor  force  has 


been  growing  by  about  one  and  one-half  percent  per  year.  Projecting 


these  trends  forward  three  years, 


(Ne  +  Nf)i  =  (6.3  +  4.0)(.99)^  =  9.99 


(Nl)3  =  (66.6) (1.015)-  =  69.64 


where  the  subscript  3  denotes  three  years  in  the  future.  If  the  accept¬ 


able  rate  of  unemployment  is  set  (arbitrarily  for  this  example)  at  one 


percent,  an  expression  for  the  unemployment  rate  in  three  years  can  be 


written 


(N,.)-(Nf+Nf)3-(Nw)3  =  0.01 


Solving  for  (Nw)3  yields 


(Nw)3  =  (1-0.01)  (Nl.)3  -  (Ne  +  Nf)3 


(1-0.01) (i9. 64)  -  (9.99) 


»  58.71 


Setting  the  response  surface  equation  for  Nu  equal  to  this  value  yields 
Eq  (6.8). 

The  left  side  of  the  second  constraint,  Eq  (6.12),  is  the  response 
surface  for  the  price  index.  The  right  side  of  the  inequality  is  the 
currently  forecast  price  index  nultiplied  by  a  five  percent  per  year 
increase  for  each  of  three  years.  This  constraint  keeps  inflation  below 
an  average  of  five  percent  per  year.  The  remaining  constraints,  Eqs 
(6.13)  through  (6.20),  are  the  political  constraints  on  fiscal  policy 
variables.  The  right  hand  side  values  of  the  inequalities  are  1952 
levels  of  the  exogenous  variables  plus  or  minus  the  amount  by  which  the 
variables  can  be  feasibly  changed. 

Solving  this  linear  programming  problem  gives  the  optimal  fiscal 
policy  to  be  followed  by  the  policy  maker.  Table  (6.2)  shows  the  solu- 


Table  6.2.  Optimal  Fiscal  Policy  for  the  Example  Problem 


haximum  Attainable  Q: 


*165.8  billion 


Fiscal  Policy  Variable  Values 


3 


The  (nodel  forecasts  gross  national  product  to  be  >172.0  billion  at 
the  end  of  the  current  year.  The  eaxiAum  attainable  gross  national 
product,  185.8,  translates  to  an  average  growth  rate  of  2.6  percent  over 
three  years.  The  average  inflation  rate  is  4.2  percent,  which  eeans  that 
there  is  "slack"  in  the  inflation  constraint.  Appendix  F  contains  the 
output  from  the  linear  prograsAing  conputer  routine  for  this  problem. 

This  solution  suggests  that  the  best  fiscal  policy  is  to  cut 
corporate  taxes,  raise  wage  taxes,  hire  more  government  employees  (or 
just  pay  them  more)  and  increase  government  expenditures  slightly.  This 
solution  sounds  fairly  plausible,  but  one  might  wonder  why  this  particu¬ 
lar  solution  IS  optimal.  Furthermore,  one  might  wonder  if  the  optimal 
solution  for  the  linear  programming  problem  is  in  fact  the  optimal 
solution  for  the  actual  Klein-Goldberger  model. 

The  answer  to  the  first  question  requires  an  examination  of  the 
coefficients  in  Eqs  (6.7),  (6.8),  and  (6.9).  The  employment  constraint 
is  always  binding  because  it  is  met  with  equality.  Economic  theory  and 
Eqs  (6.7)  through  (6.9)  indicate  that  either  decreasing  taxes  or 
increasing  government  spending  raises  gross  national  product,  employ¬ 
ment,  and  prices.  Increasing  government  wage  bill  decreases  gross 
national  product,  employment,  and  prices  according  to  Eqs  (6.7)  through 
(6.9).  This  IS  counter  intuitive,  but  Soldberger  explained  somewhat 
unconvincingly  that  increases  in  Ws  with  G  constant  represented  "a  shift 
in  the  composition  of  government  expenditures  from  business  produced 


goods  to  purchases  of  labor  services."  (Goldberger,  1959:30)  Because 


the  inflation  constraint  is  not  the  Uniting  factor,  it  is  ignored  for 
the  moment.  The  objective,  then,  is  to  find  the  feasible  combination  of 
fiscal  policy  vaiables  which  maximizes  GNP  for  a  given  level  of  employ¬ 
ment.  To  do  this  one  would  want  to  change  the  fiscal  policy  variable 
with  the  largest  increase  in  gross  national  product  per  unit  increase  in 
employment,  the  fiscal  policy  variable  with  the  next  largest  increase, 
and  so  on  until  the  required  employment  level  is  reached.  For  example, 
the  change  in  Q  per  change  in  Nu  brought  about  by  changes  in  government 
spending  is 


a  Q  =  ^6 

a  Nu  ^  G  ^  Nu 


(6.21) 


Earlier  in  this  chapter  it  was  shown  that  the  quantity  Q/  G  is  simply 
the  coefficient  of  G  in  the  Q  linear  response  function,  Eq  (6.7).  In 
addition,  Q/  Nu  is  the  reciprocal  of  the  6  coefficient  in  the  Nu 
response  surface  equation,  Eq  (6.8).  For  example  the  ratio  between  the 
G  coefficient  in  Eq.(6.7)  and  the  G  coefficient  in  Eq  (6.8)  is 


do  dG  =  4.5025  =  2.4985 

"jG  .3Nu  1.8021 


(6.22) 


The  number  2.4985  gives  the  increase  in  Q  which  occurs  when  G  increase 
enough  to  raise  Nu  by  one  unit.  Similar  ratios  can  be  computed  for  the 
other  factors. 

>Q_  ^T,,  =  -3.  1087  =  2.  4788 

dfu  ,»Nw  -1.2541 

’"nc  =  -3^3323_^  =  2.6516 
.jTc  "^Nw  “1*2567 

=  -0*  2854  =  0.5633 


■ 

>Q 

^T,,  =  -3.1087 

dTu 

.*Nw  -1.2541 

jQ 

''nc  =  -3.3323 

A 

w 

■  3  T  c 

'-*Nu  -1.2567 

'*Wr  =  -0.  2854 

'-.W, 

3Nu  -0.5019 
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Gross  national  product  increases  nost  for  a  given  level  of  employ¬ 
ment  by  a  cut  in  corporate  taxes,  then  by  an  increase  in  government 
spending,  then  by  a  cut  in  wage  taxes,  and  finally  by  a  cut  in  the 
government  wage  bill.  The  optimal  solution  sets  corporate  taxes  at  the 
lower  limit,  government  spending  at  an  intermediate  level,  and  wage 
taxes  and  government  wage  bill  at  the  high  limits.  Thus  the  given  opti¬ 
mal  solution  for  the  linear  programming  problem  does  seem  reasonable. 
However,  a  question  still  remains  as  to  whether  the  optimal  solution 
for  the  linear  programming  problem  is  optimal  for  the  actual  Klein- 
Goldberger  model. 

Chapter  V  shows  that  the  response  surfaces  do  in  fact  closely 
approximate  what  is  going  on  in  the  model  over  the  entire  range  of  data. 
Furthermore,  higher  order  terms  are  not  necessary  to  obtain  a  good  fit. 
Therefore,  what  is  optimal  for  the  response  surface  model  of  the  economy 
should  be  optimal  for  the  KG  model.  Verifying  this  assertion  requires 
searching  the  area  around  the  alleged  optimal  solution  to  see  if  further 
gains  might  be  made  with  an  alternate  policy.  This  search  is  to  be  done 
with  the  original  model.  If  the  solution  given  for  the  linear  orogram- 
ming  problem  is  not  the  optimal  then  one  should  be  able  to  increase 
gross  national  product  and  satisfy  the  constraints  by  adjusting  Tw,  Tc, 
G,  or  W?.  The  table  below  shows  the  results  of  running  the  KG  model 
with  the  fiscal  policy  variables  set  at  values  slightly  different  than 
the  optimal  policy  determined  by  the  response  surface  derived  linear 
programming  problem. 


*  JfcV  A.  L-V  ^ 
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Table  6.3.  Klein-Galdberger  Model  Solutions  in  the  Area  of 
the  Alleged  Qptiital  Solution 


58.71  240.5 
nax  max 


T« 

Tc 

G 

Us 

Nw 

P 

Q 

11.63 

7.14 

44.54 

21.82 

56.5 

222.5 

186.2 

10.63 

7. 14 

44.54 

21.82 

57.7 

226.6 

189.3* 

11.63 

8. 14 

44.54 

21.82 

55.2 

217.6 

182.8 

11.63 

7. 14 

43.54 

21.82 

57.3 

215.0 

181.6 

11.63 

7. 14 

45.54 

21.82 

58.3 

228.7 

190.8* 

11.63 

7. 14 

44.54 

20.82 

57.0 

222.5 

186.6* 

10.63 

8. 14 

44.54 

21.82 

56.5 

222.5 

186.0 

10.63 

7. 14 

43.54 

21.82 

55.9 

219.9 

184.8 

11.63 

8. 14 

45.54 

21.82 

57.0 

224.5 

187.4* 

11.63 

8. 14 

45.54 

20.82 

57.6 

229.5 

187.7 

5.63 

13.14 

44.54 

21.82 

56.5 

224.6 

185.0 

11.63 

13. 14 

61.70 

16.82 

82.4 

337.3 

245.0 

The  first  line  in  Table  6.3  is  the  alleged  optimal  solution. 
However,  the  starred  solutions  yield  greater  gross  national  product 
values  than  the  supposed  optimal  solution.  The  solution  computed  by  the 
linear  programming  algorithm  is  not  optimal  because  the  value  of  Nm 
forecast  by  the  response  surface  function  was  four  percent  too  low.  The 
four  percent  error  is  not  unreasonable  because  although  the  linear 
approximation  to  the  KG  model  is  good,  it  is  not  perfect. 

All  IS  not  lost,  however,  because  the  response  surface  coefficients 
can  be  used  to  "tweek"  the  solution  to  optimality.  As  pointed  out 
above,  economic  theory  and  Eqs  (6.7),  (6.8),  and  (6.9)  indicate  that 
when  corporate  or  wage  taxes  decrease,  gross  national  product  increases, 
prices  increase,  and  employment  increases.  The  effect  of  increased 
government  spending  is  the  same.  Increased  wage  bill  has  a  small  down¬ 
ward  effect  on  gross  national  product  and  employment,  and  a  substantial 
downwara  effect  on  prices  according  to  Eqs  (6.7),  (6.8),  and  i6.9). 


To  increase  Q,  one  could  decrease  taxes  or  governnent  wage  bill  or 
increase  government  spending.  However,  the  same  measures  which  raise  Q 
also  raise  Nw  and  p.  Nu  must  increase  by  58.71  -  56.50  »  2.21  to 
satisfy  the  employment  constraint  with  equality  and  p  may  increase  by 
240.5  -  222.5  -  18.  The  best  variable  to  alter  is  the  variable  which 
increases  Q  the  most  for  the  required  change  in  Nw  without  violating  the 
price  index  constraint.  The  discussion  above  shows  how  to  compute  the 
ratios  for  tradeoffs  between  employment  and  gross  national  product.  For 
example,  if  Nw  must  increase  by  2.21  to  meet  the  employment  constraint, 
then  the  total  change  in  Q  due  to  a  change  in  G  is  (2. 4985) (2. 21 )  » 
5.5217  using  the  ratio  computed  in  Eq  (6.22).  To  find  how  much  6  must 
increase  to  raise  Nw  to  the  required  level,  one  can  divide  the  required 
change  by  the  G  coefficient  in  Eq  (6.8)  (i.e.,  2.21/1.8021  =  $1,226 
billion).  By  Eq  (6.9)  this  increase  in  6  induces  a  price  index  increase 
of  ( 1 . 226) (6. 6816)  =  8.  194.  Since  this  increase  in  G  would  only 
increase  the  price  index  to  222,5  +-8.194  =  230.694,  this  solution  is 
feasible  because  this  price  index  is  below  the  240.5  value  allowed  by  Eq 
(6.  12) . 

To  find  the  best  factor  to  change,  tradeoff  ratios  for  Tw  and 
must  be  compared  to  the  6  tradeoff  ratio.  The  factor  with  the  largest 
tradeoff  ratio  which  does  not  cause  the  price  index  to  exceed  its  maxi¬ 
mum  is  the  best.  Because  Tc  is  already  at  the  lower  limit,  it  need  not 
be  investigated  for  alteration.  Table  6.4  summarizes  the  data  required 
to  select  the  best  factor  to  adjust.  The  first  column  is  the  tradeoff 
ratio.  The  next  column  is  the  change  in  Nw  required  to  satisfy  the 
emolovment  constraint.  The  Q  column  is  simply  the  product  between  the 
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Urst  and  second  columns.  It  is  the  change  in  Q  resulting  when  the 
factor  is  altered  enough  to  bring  about  the  required  change  in  Nw.  The 
F  column  is  the  change  in  the  factor  required  to  increase  Nu  by  2.21. 

It  is  equal  to  the  Nw  column  divided  by  the  factor's  coefficient  in  Eq 
(6.8).  The  p  column  is  the  increase  in  p  caused  by  the  increase  in  the 
factor.  It  is  equal  to  the  F  column  times  the  factor's  coefficient  in 
Eq  (6.10).  Finally  the  p  column  is  the  new  price  index  brought  about  by 


changing 

the  factor 

to  its 

new  level. 

If  the 

figure  in  the 

p  column 

exceeds 

240.5,  the 

solution 

IS  infeasible. 

Table  6.4. 

Data  for 

Selecting 

the  Best 

Variable  to 

Alter 

Q/  Nw 

Nw 

Q 

F 

P 

P 

Tw 

2.4788 

2.21 

5.4781 

-1.7622 

a. 1151 

230.6 

G 

2.4985 

2.21 

5.5217 

1.2263 

a. 1936 

230.7 

M2 

0.5685 

2.21 

1.2563 

-4.4026 

24.897 

247.4 

A  quick  scan  of  Table  6.4  reveals  that  increasing  G  by  $1.2263 
billion  increases  Q  by  $5.5217  billion  while  a  decrease  in  Tw  of  $1.7622 
billion  increases  Q  by  only  $5.4781  billion.  Decreasing  Ws  by  the 
amount  required  to  increase  Nw  by  2.21  million  workers  causes  the  price 
index  to  exceed  the  maximum.  Thus,  the  adjusted  optimal  solution  is 

Tw  =  11.63 
Tr.  =  7.14 
G  =  45.78 
Wt  =  21.02 

A  single  solution  for  a  linear  programming  problem  is  rarely  very 
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useful  without  sensitivity  analysis.  Fortunately,  sensitivity  analysis 
for  linear  progr aoiini ng  problems  is  very  well  developed.  For  instance  the 
shadow  prices  tell  how  much  the  objective  function  will  change  if  the 
right  hand  side  of  a  constraint  is  changed.  Table  6.5  shows  shadow 
prices  for  each  binding  constraint. 


Table  6.5.  Shadow  Prices  for  the  Fiscal  Policy  Problem. 


Constraint 


Maximum  Tm 


Objective  Function  Change 


0.0254 


Maximum  W: 


0.9606 


Minimum  Tm 


0. 1925 


Minimum  Mz 


2.4905 


Another  option  for  conducting  sensitivity  analysis  for  this  example 
linear  programming  problem  is  (believe  it  or  not)  response  surface 
methodology.  A  new  response  surface  can  be  built  for  Q  in  terras  of  Tm, 
Tc.,  3,  nil,  inflation  and  unemployment  by  varying  right  hand  sides  of  the 
constraints  in  accoroance  with  an  experimental  design  and  solving  the 
linear  programming  problem  (Smith  and  Mellichamp,  1979). 

Because  the  KG  model  is  nearly  linear,  first  order  response  surface 
equations  fit  the  model  fairly  well.  The  linear  objective  function  and 
constraints  make  it  oossible  to  formulate  an  optimization  problem  as  a 
linear  program.  If  the  model  were  not  so  linear  and  the  response  sur¬ 
faces  nad  higher  order  terms,  an  optimization  problem  could  still  be 
formulated  and  solved  using  nonlinear  ootimization  techniques  available. 
One  computer  i mo  1 ementat i on  of  nonlinear  techniques  is  the  Sequential 
Unconstrained  Minimization  Technique  iSUMT)  package.  The  program  handles 
nonlinear  objective  functions  and  constraints  with  inequalities.  The 
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This  study  has  uncovered  such  about  Mhat  can  and  cannot  be  done  in 
applying  response  surface  nethodology  to  a  eacroecononic  nodel. 

Response  surface  aethodology  is  a  useful  tool  Mhich  can  be  used  to 
investigate  the  properties  of  a  eacroecononic  nodel  as  long  as  the 
Imitations  of  the  nethod  are  kept  in  nind. 

The  research  denonstrates  hOM  to  generate  response  surfaces  fron  a 
nacroecononic  nodel.  First,  the  problen  to  be  addressed  is  defined  and 
a  deternination  is  nade  that  response  surface  nethodology  is  the  appro¬ 
priate  tool  for  solving  the  problen.  Next,  variables  of  interest  are 
selected  and  their  operating  ranges  specified.  The  forn  of  the  response 
surface  equations  is  decided  upon  and  an  appropriate  expermental  design 
selected.  After  translating  the  coded  expermental  design  to  actual 
predetermined  variable  values,  the  nodei  is  solved  for  each  conbination 
of  predeternined  variable  values  specified  by  the  experinental  design. 
The  data  generated  are  used  to  estinate  response  surface  coefficients  in 
terns  of  the  coded  experinental  design  to  preserve  orthogonality.  The 
response  surface  fit  is  checked  and  the  response  surface  equation  is 
decoded  so  that  it  is  expressed  in  terms  of  the  original  variables. 

After  generating  the  resoonse  surface,  analyses  may  begin. 

This  study  shows  that  a  low  order  polvnonial  can  indeed  fit  the 
responses  of  the  Klein-Goldberger  econometric  model.  The  near  linearity 
of  the  model  is  confirmed.  The  coefficients  from  a  decoded  first  order 
response  surface  fit  to  the  Klein-Goldberger  model  are  compared  to 
multipliers  computed  by  Goldberger  and  found  to  correspond  quite 


clofsly. 

EconoAists  use  Aultipliers  extensively  to  characterize  the  static 

* 

;  and  dynamic  properties  of  econometric  models  and  to  conduct  policy 

I  simulations.  Techniques  for  generating  aultipliers  include  lineariza- 

I  tion  of  the  econometric  model  through  the  derivative  method  and  comput- 

I  ing  control  and  disturbed  solutions.  Response  surface  methodology 

j  offers  another  alternative  for  computing  multipliers. 

,  Multipliers  are  derived  from  response  surface  equations  by  taking 

r 

I  the  partial  derivative  of  the  response  surface  equation  Nith  respect  to 

f 

j  the  variable  of  interest.  The  result  may  or  may  not  be  a  constant. 

Higher  order  terms  in  the  response  surface  equation  cause  multipliers  to 
be  deoendent  on  the  level  of  one  or  more  variables.  Consequently, 

^  multipliers  obtained  from  response  surfaces  are  most  useful  for  investi¬ 

gating  nonlinear  econometric  models. 

The  advantage  of  response  surface  derived  multipliers  over  multi¬ 
pliers  derived  by  model  linearization  is  that  the  model  does  not  have  to 
be  linear  or  near  linear  for  multipliers  to  be  valid  over  a  wide  range 
of  variable  values.  In  addition,  significant  interactions  and  higher 
order  effects  can  be  identified.  The  advantage  of  response  surface 
multipliers  over  control -di sturbed  multipliers  is  that  thev  more 
completely  characterize  the  relationships  in  the  model  and  reduce  the 
number  of  runs  reauired  to  estimate  multipliers.  Also,  multipliers 
'  computed  Pv  the  control -d i sturbed  method  do  not  identify  interactions 

ana  nigher  order  effects.  These  multipliers  are  onlv  valid  for  small 
ranges  of  the  predetermined  variables. 


Response  surface  multipliers  can  be  used  in  the  same  wavs  that 


Aultiplitrs  dirivtd  by  othir  ttchntqutt  an  ustd.  Usis  includi  policy 
simulation  and  dcttrmining  moat  influmntial  factors  in  the  economy. 

Response  surfaces  can  be  used  for  more  than  just  computing  multi¬ 
pliers.  They  can  also  be  used  to  formulate  optimization  problems.  The 
explicit  nature  of  the  response  surface  equation  giving  endogenous 
variables  in  terms  of  predetermined  variables  facilitates  optimization 
problem  formulation.  Chapter  VI  gives  an  example  problem  in  which  gross 
national  product  is  maximized  while  holding  unemployment  and  inflation 
at  or  below  acceptable  levels.  The  problem  is  formulated  as  a  linear 
programming  problem  and  solved.  Optimization  problems  with  nonlinear 
response  surface  derived  constraints  and  objective  functions  can  also  be 
formulated  and  solved  using  optimization  packages  such  as  SUMT. 

Applying  response  surface  methodology  to  macroeconomic  models  is 
not  without  limitations.  Computing  multipliers  using  response  surface 
methodology  is  more  cumbersome  than  existing  methods  for  deriving  multi¬ 
pliers  from  linear  or  near  linear  models.  Separate  response  surfaces 
must  be  computed  for  each  response  variable  for  each  time  period  for 
each  subset  of  predetermined  variables.  For  nonlinear  models,  response 
surface  multipliers  which  are  functions  better  characterize  input-output 
relationships  than  traditional  multipliers.  However,  special  care  must 
be  taken  to  insure  response  surface  fit  before  drawing  inferences  about 
the  iiiodel  based  on  the  response  surface  generated. 

A  limitation  which  detracts  from  using  response  surfaces  for 
prediction  is  the  few  number  of  factors  which  can  be  inducted  in  the 
response  surface  function.  The  number  of  variables  which  can  be 
included  in  the  response  surface  equation  is  limited  by  the  size  of  the 
experimental  design.  Of  particular  concern  are  interactions  between 
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exogenous  variables  included  in  the  response  surface  and  lagged  endoge¬ 
nous  variables  not  included  in  the  response  surface.  It  is  quite  possi¬ 
ble  that  the  effect  of  exogenous  changes  depends  on  the  current  state  of 
the  economy.  If  variables  are  omitted  from  the  response  surface  equa¬ 
tion,  then  the  response  surfaces  only  capture  model  relationships  at  the 
specific  levels  assumed  in  generating  the  response  surface.  The  other 
methods  of  computing  multipliers  suffer  from  this  deficiency  too. 

On  balance,  the  limitations  of  response  surface  methodology  do  not 
preclude  it  from  being  a  valuable  tool  for  analyzing  certain  aspects  of 
macroeconomic  models. 

Recommended  Further  Research 

There  are  several  areas  available  for  further  research.  First,  it 
nas  oeen  assumed  that  a  low  order  polynomial  could  adequately  fit  an 
econometric  model  more  nonlinear  than  the  Klein-Soldberger  model.  This 
assumption  needs  testing.  Second,  including  more  variables  in  the  model 
bv  using  larger  experimental  designs  has  yet  to  be  explored.  Large 
exoerimental  designs  could  be  develooed  by  computer  algorithm.  A  final 
improvement  of  the  research  presented  in  this  thesis  would  be  to  include 
time  as  an  independent  variable  in  the  response  surface  equation. 
Including  time  in  the  response  surface  equation  would  eliminate  the  need 
to  generate  response  surfaces  for  each  period.  Preliminary  attempts  to 
fit  response  functions  with  time  as  an  independent  variable  to  the 
Kl ei n-Gol dberger  model  yielded  R*  values  of  0.9d10  for  number  of  workers 
empioved.  0.9858  for  price  index,  and  0.9630  for  gross  national  product. 
To  generate  these  response  surfaces,  a  second  order  polvnomial  including 
wage  taxes,  corporate  taxes,  government  nonwage  expenditures,  government 


wage  bill,  number  of  government  employees  and  time  was  fit  to  data  from 
periods  zero  through  five.  Developing  response  surfaces  with  time  as  an 
independent  variable  would  reveal  time  delay  aspects  of  the  econometric 
model  which  could  prove  quite  valuable. 
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Appendix  A:  SQl.ving  the  KLein^goiggtrger  Madel 


Introduction 

Because  the  Klein-6oldberger  eodel  is  not  linear  and  has  no  simple 
analytical  solution,  solving  the  model  requires  an  solution  approxima¬ 
tion  technique.  Goldberger  used  a  derivative  method  to  obtain  a  linear 
approximation  to  the  model  formulated  in  terms  of  variable  changes.  He 
then  derived  impact  and  interim  multipliers  from  the  linear  approxima¬ 
tion.  This  thesis  requires  a  method  for  solving  the  Klein-Goldfaerger 
model  Mithout  linearizing  it.  Economists  commonly  use  some  sort  of 
numerical  technique  for  solving  nonlinear  econometric  models.  Klein 
recommended  the  Gauss-Sidel  numerical  method  for  solving  econometric 
models  (Klein,  1974:238).  The  method  is  a  simple  iterative  procedure 
Mhicn  does  not  require  derivative  computations.  This  appendix  describes 
the  method,  illustrates  it  with  an  example,  and  shows  how  the  method  was 
aoplied  to  solve  the  Klein-Goldberger  model  for  this  research  effort. 
P®§£ciBii9D 

Klein  describes  the  Gauss-Sidel  method  in  ms  text,  A  lextbgqk  of 
i£9D9!!!5£Ei£S  (Klein,  1974:238-239).  The  material  below  restates  Klein’s 
description.  An  econometric  model  with  n  current  endogenous  variables, 
n  endogenous  variables  lagged  up  to  p  periods,  and  m  exogenous  variables 
lagged  up  to  p  periods  can  be  written  in  the  form 


y . .  f. 


g.  ( V 


y  r»  ,  ♦  -  «  f 
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,t  (Aal) 


1  -  (1.2,. ...n) 


where 


vi.t  «  one  0^  n  current  endoenous  variables, 


Vi.t-n  -  one  of  n  endogenous  variables  lagged  p  periods, 

3  one  of  i»  current  exogenous  variables, 

3  one  of  ID  exogenous  variables  lagged  p  periods,  and 
Bi.t  s  a  randoe  error  tere. 

In  Mords,  Eg  (A.l)  says  that  each  equation  in  the  model  gives  a  single 
current  endogenous  variable  (which  must  be  solved  for)  as  a  function  of 
the  other  current  endogenous  variables,  lagged  endogenous  variables,  and 
exogenous  variables.  Lagged  endogenous  variables  together  with  current 
and  lagged  exogenous  variables  are  known  as  predetermined  variables. 

Mith  a  few  exceptions,  the  Klein-Goldberger  model  in  Table  2.1  has  the 
form  of  £q  (A.l).  It  is  usually  possible  to  rewrite  Eq  (A.l)  in  the 
f  orm 

Vi.t  "  g.(vt,t,...yr.,t.yi,t  —  t,...yn.t— 0,...,Xt,t,...Xii..t  —  p)^ei,t  (A. 2) 

1  3  (1,2.. ..,n) 

Eq  '.A. 2)  IS  the  same  as  Eo  (A.l)  except  in  Eq  (A. 2)  y,.t  appears  on  both 
sides  of  the  equation.  Omitting  the  error  term,  e, ,  and  inserting 
superscripts  in  accordance  with  the  Gauss-Sidel  method  converts  Eq  (A. 2) 
to  an  algorithm. 


V ,  . ,  ’  •  ‘  '  =  g .  (  V I  ,  - 


V,,t  —  l....yr.,,-p.  Xl,t....Xm,,— p)  (A.<j) 


I  =  '1.2 . n ) 

where 


*  the  value  of  the  ith  current  endogenous  variable  from 
the  ir*l)th  iteration  of  the  method,  and 


St 


*  th«  valus  ot  tha  ith  current  endogenous  variable  tron  the 
rth  Iteration. 


Iterations  are  per^orned  until 


'yi 


<  tolerance 


A  simple  example  illustrates  the  confusing  notation  in  Eqs  (A.l),  (A. 2), 
and  (A. 3).  The  system  of  nonlinear  equations, 


X  =  -4z  +  2w  ♦  6 


y  .  4^  1/3  +  9 

2  *  x/y  -  2m 


(A. 4) 


where 


X,  V,  2  «  variables  to  be  solved  for,  and 
w  3  j  variable  whose  value  is  known. 


IS  in  the  form  of  Eq  (A.l).  Eqs  (A. 4)  through  (A. 6)  can  be  rewritten  in 
the  form  of  Eq  (A. 2)  by  multiplying  both  sides  of  the  equations  by  a 
constant,  sav  0.5,  and  then  adding  (1-0.5)  times  the  left-hand  side 
variables  to  both  sides  of  the  equations. 


X  =  0.5(-42  +  2w  ♦  6)  ♦  0.5> 
V  =  0. 5  (4x  ' '  '  ♦  8)  *  0.  5v 
2  =  0.5(x/y  -  2w/  *  0.52 


(A. 7) 


(A. 8) 


Bv  arranging  terms  and  inserting  superscripts  denoting  iterations,  Eqs 
(A. 7;  through  iA.9)  become  algorithms  for  computing  a  solution. 


>  O-Sx'*-*  -  2:“"  +  Iw  ♦  3  (A. 10) 

y'"**'  »  2x-‘'=*'-*‘*  +  0.5y‘^‘  ♦  4  (A. 11) 

z‘-*“  s  0.5(x  )  -  w  ♦  0.5z‘"’  (A. 12) 

The  Qauss-Sidel  nethod  requires  an  initial  solution,  a  specification  for 
M,  and  a  specification  of  the  error  tolerance.  If  x‘“’  «  0,  y‘°'  «  0, 
z‘o>  =  0,  w=5,  and  tolerance®. 01  then  the  first  iteration  of  algorithm 
Eqs  (A. 10)  through  (A. 12)  is 


X ' ‘ '  =  (0.5)  (0)  -  2(0)  ♦  1 (5)  +3  =  8 

v'  “  =  2(8)  +  (0.5)  (0)  +  4  =  4 

z '  “  =  0. 5(8/4)  -  5  ♦  (0.5) (0)  =  -4 

y<t>=4  j<l'a-4 

The  second  iteration  is 

X '  =  >  =  (0.5)  (8)  -  2(4)  +5  +  3  =  4 
y <3*  s  2(-8) ‘^^  +  (0.5) (4)  +  4  =  2 

»  0.5(-8/-l)  -  5  +  (0.5)(-4)  =  -3 
X‘l>s4  y<i:'a2  z'  =  ’=-3 
Iterations  continue  until 

-  x""*!  <  0.01  ly..*.!!!! _ "_yJ.M.L  ^ 

I  X  • "  '  I  I  V  ‘  '  I 

and 

I  \ 


The  algorithm  is  not  guaranteed  to  converge  for  all  forms  of  Eg  (A. 2) 
and  for  ail  trial  solutions.  There  are  convergence  conditions  for  the 
aauss-Sidel  method,  but  often  the  convergence  conditions  are  extremelv 


difficult  to  comoute.  In  practice,  trial  and  error  usually  reveals 


simple  forms  that  converge. 

The  method  described  above  computes  a  solution  for  current 
endogenous  variables  given  lagged  endogenous  variables  and  exogenous 
variables.  The  method  easily  adapts  to  compute  extended  period 
forecasts.  Once  the  method  yields  a  solution,  lagged  variables  are 
updated  with  current  variable  values  and  exogenous  variables  are  set  to 
values  dictated  by  the  policy  under  investigation.  Then  the  model  is 
resolved.  In  the  notation  of  Eg  (A.l), 


for 

all 

gfl 

and 

k  =  1 , 2 , . . 

•  »P 

for 

all 

i*lg2g*«e  gO 

and 

k=l,2,.. 

.  .P 

The  Klein-Goldberger  model  was  solved  using  the  Gauss-Sidel  method.  The 
K1 ei n-Gol dberger  model  equations  in  Table  2.1  were  put  in  the  form  of  Ed 
(A.l)  with  one  current  endogenous  variable  expressed  as  a  function  of 
the  other  current  endogenous  variables,  lagged  endogenous  variables,  and 
exogenous  variables.  For  the  sake  of  simplified  discussions,  Eq  (A.l) 
can  oe  abbreviated, 


y.  =  g.  (A. 15) 

where  g,  is  a  function  of  the  other  current  endogenous  variables,  lagged 
endogenous  variables,  and  exogenous  variables.  Performing  some  simple 
algebraic  manipulation  converts  Eq  (A. 13)  to  the  form  of  Eq  iA,2).  If  a 
IS  a  constant, 

ay.  =  a  g.  (A. 14) 


av.  (i-a)y,  =  ag,  +  (l-a)y, 


(A.  15) 


k  . 

t . 

k- 


*•_ 


control  language,  an  initial  trial  solution,  values  tor  predetermined 
v'ariaoles  to  be  included  in  the  response  surtace  equation,  values  tor 
oredetermined  variables  not  to  be  included  in  the  response  surface 
equation,  and  coded  values  (-1,  0,  or  i)  tor  the  variables  to  be  in¬ 
cluded  in  the  response  surtace  equation.  The  program  outputs  a  tile 
containing  the  coded  variable  values,  a  case  number,  the  torecast  period 
number,  and  the  solution  tor  endogenous  variables.  The  output  tile  may 
be  read  directly  by  the  BMDP  statistical  package  tor  response  surtace 
coetticient  estimation. 

The  program  tirst  reads  and  echoes  the  control  language  contained 
in  the  tile  "kg.ctl."  Atter  initializing  arrays,  the  program  enters  a 
loop  which  solves  the  rl ei n-Gol dberger  model  tor  each  set  ot  variaole 
levels  soecitied  case  bv  case  tor  the  number  ot  periods  specitied.  The 
program  tirst  reads  one  set  ot  coded  values  (-1,  0,  or  1)  tor  variables 
to  be  included  in  the  response  surtace  eauation  (variables  which  will 


.'Iv! 


m 


i 
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hence'forth  be  referred  to  as  desiQn  variables).  Subroutine  setbas  reads 
predetermined  variable  values  into  the  arrays  P,  IS,  and  PI.  Arrav  P 
contains  exogenous  variable  values  and  endogenous  variable  values  lagged 
one  period.  IS  and  PI  contain  values  of  short  term  interest  rates  and 
prices  lagged  more  than  one  period.  Setbas  also  reads  a  trial  solution 
into  XO.  Next,  subroutine  setdes  resets  the  design  variables  to  the 
levels  specified  in  the  file  "design. cod. "  Subroutine  solve  calls 
subroutine  iterate  in  a  loop  to  compute  Gauss-Sidel  iterations  of  the 
Klein-Soldberger  model  until  the  solution  converges  within  tolerance. 

If  the  control  language  specifies  that  intermediate  period  solutions  are 
to  be  printed,  the  program  writes  the  coded  design  variable  levels,  case 
number,  period  number,  and  the  intermediate  solution  to  the  file 
"kg. out."  The  last  period's  coded  design  variable  levels,  case  number, 
period  number,  and  solution  are  always  printed.  If  the  control  language 
specifies  that  extended  period  forecasts  are  to  be  made,  subroutine 
update  updates  lagged  variables  and  a  new  solution  is  computed.  The 
program  stops  when  solutions  are  computed  for  the  number  of  periods 
specified  for  each  set  of  design  variable  level  specifications. 


*******«*«******#««*«**»<*»«•#««*«*«*«*««•#»«*«*********»*«««*#*«««*««« 

program  kgsolv 

*  Solves  the  K1 ei n-6ol dber ger  model  using  the  Gauss-Sidel  numerical 

♦  technique 

double  precision  XO (21 ) , X 1(21 ) ,P (44 ) , IS (5) , PI (2) 
real  tol ,  C0D0ES(5) 

integer  dl,  d2,  i,  lease,  ioriod,  liter,  ncase,  npriod,  itmax, 

Si  prtall,  guessl,  bascas.  design,  kgout,  kgctl,  codes 

opend, files  'guessl.dat  ) 
rewind ( 1 ) 

ooeni2, files  'bascas.dat') 
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r«win(l(2) 
optn (3,f ilt> 
rewind (3) 
open (4,f ile> 
open (8, file* 
rewind (8) 
open(9,file« 
rewind (9) 
guessl*! 
bascas*2 
design*3 
kgout*4 
l<gctl*8 
codes*9 

*  Initialization 


'defign.dat ' ) 

‘  kg. out ‘ , status* 'new ' ) 
‘ kg. ctl ' ) 

design. cod ' ) 


*  Read  and  echo  control  language 

write(*,»)  '  CONTROL  DATA’ 
read (kgctl ,*)prtall 

write(*,»)'  PRINT  EACH  PERIOD  DATA?  (leYES)  ,ortall 

read (kgctl ,*)ncase 

Mrite(*,»)'  tt  CASES*  ',ncase 

read (kgctl ,»)npriod 

write(*,#)'  4  PERIODS*  ',npriQd 

read (kgctl ,*) tol 

write(*,*l'  TOLERANCE*  ',tol 

read (kgctl , *) i tnax 

write(*,*)'  MAX  ITERATIONS*  ',itiitax 

icase*l 

ipr iod=0 

iiter=l 

*  Initialize  arrays 

do  100  i=l,21 
X0(i )»0 
X 1 ( i ) *0 

100  continue 

do  200  i=l,44 
P(i )=0 

200  continue 

*  Strip  off  data  dimensions  with  dummy  variables. 

read ( codes ,*) d 1 .  d2 

*  Main  Program 

*  While  the  current  case  is  less  than  the  last  case 


400 


If  ucase. le.ncase)  then 


rewind (gucssl ) 
rewind (bascae) 

»  Read  coded  design  variable  levels 

read (codes,*) (CODDES (i ) fi^l  ,S) 

*  Read  predetereined  variable  levels 

call  set bas (guess  1 ,bascas,XO,P, IS,PI > 

*  Read  design  variable  levels 

call  setdes(design,P, IS,PI) 
ipriQd»0 

*  while  current  period  is  less  than  or  equal  to  the  forecast  period 

300  if  (ipriod.le.npriod)  then 

call  solve(X0,Xl , P, i case, ipriod, i t max , tol , i i ter > 
if  (prtall.eq. 1)  then 

write(kgaut,1020) (CODDES (i ) , i =1 ,5) 
wri te (kgout , »)  icase,  ipriod 
write(kgout, 10 20 )(Xl(i), 1=1,21) 
else  if (ipriod. eq.npriod)  then 

wri te (kgout ,  1020) (CODDES (i ) , i  =  1 ,5) 
write(kgout,*)  icase,  ipriod 
write (kgout ,1020)  (Xl(i),i  =  l,21) 
endi  f 

1020  format(lx,5F12.6) 


call  update(X0,Xl ,P, IS,PI) 

ipriod=ipriod+l 
goto  300 
endif 

*  end  while  (ipriod) 
icase=icase*l 
goto  400 

endif 

*  end  while  (icase) 
write(*,*) 

write(*,*)  PROGRAM  COMPLETE,  RESULTS  IN  KG. OUT 

stop 

end 

*»**«««******»***«**«*#**«•»#«**««*##**«*««•«**«***«*«****«*««**« 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦*»*#sijbroutines************************************ 

**«»***«»****«***t******t*«»**«««**tt««*«**«**»**4***«**4«**»***« 
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subroutini  sstbat (guisil ,baicas,XO,P, 1S,PI) 

*  Reads  Initial  guess  solution  and  base  case  predetereined  data. 

double  precision  X0(21),  P(44),  1S(5),  PI(2> 
integer  i,  guessl,  bascas 

»  Read  initial  guess  at  solution. 

read  (guessl ,*) (X0(i),isl,21) 

*  Read  base  case  predetereined  data 

read (bascas,*) (P (i ) ,  i^i ,44) 

*  Set  up  variables  with  more  than  one  year  lag. 

IS(1)=P(29) 

IS(2)=(P(29)+P(30) )/2 
IS(3)=P(30) 

IS(4)  =  (P(30)t'P(31)  )/2 
IS(5)=P(3l) 

PI (1)=P(35) 

PI (2)=P(36) 

return 

end 

t*************************************************************** 

subroutine  setdes (design, P, IS, PI) 

*  Reads  changed  predetereined  data  for  a  new  case 

double  precision  P(44),  IS(5),  PI(2) 
integer  design 

*  Read  in  design  variable  values 

read(design,#)P(2),P(3),P(6),P(7),P(l3) 

*  P(2)=TW,  P(3)=TC,  P(6)=G,  P(7)=W2.  F(13)=NG 

*  Set  up  variables  with  eore  than  one  year  lag 

IS  ( I )  =P (29) 

IS(2)=iP(29)+P(30) ) /2 
IS(3)=P(30) 

IS(4)*(P(30)+P(31) )/2 
IS(5)=P(3l) 

PI  (  n»P(35) 

PI (2)=P!36) 


return 

end 
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subroutine  sol ve(X0, XI ,P, ic*se,ipriod ,ite«K , tol , i iter ) 

*  Coeputes  a  nueerical  solution  to  the  Klein-Soldberger  oodel 

double  precision  X0(21)  ,X1(21)  ,  P(44) 

real  tol,  ERROR (21 ) ,errorO 

integer  i,  lease,  ipriod,  itaax,  liter 

*  Conpute  nueerical  solution 

i i ter »0 

•  Repeat  until  error  <  tolerance 

400  continue 

i^  ( i i ter . 1 1 . i teax )  then 

•  call  iterate  subroutine 

call  iterat(X0,  P,  XI) 
liter*! iter ♦I 

else 

write(*,*)  CASE  ', lease,'  PERIOD  ', ipriod, 
i  ■  FAILED  TO  CONVERGE  AFTER  ' 

write(*,*)iiter,  '  ITERATIONS.  PROGRAM  STOPPED, 
stop 
endi'f 

♦  Check  current  iteration  for  tolerance 

error0“0 
do  200  1*1,21 

ERROR (i )=abs ( (X I (i ) -XO (i ) ) /XO ( i ) ) 
if  (ERRORd  )  .gt.errorO)  then 
err  or 0*ERR0R ( i ) 
end!  f 

200  continue 

if  (errorO. gt. tol )  then 
do  300  i=l,21 
X0(i)»Xl(i) 

300  continue 

goto  400 

♦  iDo  another  iteration  until  error  below  tolerance! 
else 

writel*,*)  CASE  ', lease,'  PERIOD  ', ipriod, 

S,  '  -  filter,'  ITERATIONS' 

endi  f 

return 

end 


subrautini  itiraKXO,  P,  XI) 

*  Per'forns  Gauss-Sidel  iterations 

integer  j 

double  precision  C,  I  ,SC ,PC,D,W1 ,NM,MR,FI ,R1 ,PR,L1  ,L2,  IL,  IS 
&  K,SB,PY,n,PI,Q 

*  double  precision  Cl ,  1 1 ,SC1 ,PC1 ,D1 , W1 1 ,NN1 , WRl ,PI 1 ,R1 1 ,PR1 , 

*  k  L11,L21,IL1,IS1,K1,S61,PY1,I11,P11,Q1 

double  precision  TE,TM,TC,TN,TR,6,W2,R2,T,H,NP,NL,NG,NE,FR 
i  PF,LB 

double  precision  CLl ,SCL1 ,PCL1 ,DL1 , Will ,WRL1 , Fill ,R1L1 , PRLl 
L1L1,L2L1,ISL1,ISL3,1SLS,KL1,SBL1,PL1,PIL1,PIL2,QL1 

double  precision  TELl ,TWL1 ,TCL1 ,TNL1 , TRLl , U2L1 ,R2L1 

double  precisian  X0(21) ,P(44) ,X1 (21) 

*  Initialize  variables  in  the  Klein-Goldberger  model 

*  ENDOGENOUS  VARIABLES 

C  =X0(n 

I  =X0(2) 

SCaX0(3) 

PC=X0(4) 

D  >X0(5) 

Wl>X0(6) 

NW»X0(7) 

WRsxOlS) 

FI=X0(9) 

RI»X0(10) 

PR=X0(ll) 

L1=X0(12) 

L2=X0(13) 

IL»X0(14) 

IS=X0(15) 

K  =X0(16) 

SB=X0(17) 

PY=X0(18) 

II  =X0(19) 

PI=X0(20) 

Q  »X0(21) 

*  EXOGENOUS  VARIABLES 

TE*P(l) 

TWap (2) 


TCsP(3) 

TNsP(4) 

TR=P(5) 

G  =P(6) 

W2»P(7) 

R2«P(8) 

T  »P(9) 

H  »P(10) 

NP*P(lll 

NL«P(12) 

NG»P(13) 

NE«P(14) 

FR»P(15) 

PF»P(16) 

LB»P(17) 

«  LAGGED  ENDOGENOUS  VARIABLES 

CLl  »P(18) 

SCL1=P(19) 

PCL1=P(20) 

DLl  =P(21) 

W1LI=P(22) 

WRL1»P(23) 

FIL1»P(24) 

R1L1»P(25) 

PRL1*P(26) 

L1L1«P(27) 

L2Ll=P(28) 

ISLlsP(29) 

ISL3=P(30) 

ISL5=P(3l) 

KLl  =P(32) 

SBL1=P(33) 

PLl  =P(34) 

PIL1=P(35) 

PIL2*P(36) 

QLl  =P(37) 

♦  LAGGED  EXOGENOUS  VARIABLES 

TEL1*P(38) 

TWL1=P(39) 

TCL1=P(40) 

TNLl=P(41) 

TRL1=P(42) 

W2Ll*P(43) 

R2L1=P(44) 

♦  THE  f,LEIN-GOLDBERQER  MODEL 

♦  1 

C=0.5#(-22.  26  +  0.55*(W:+)il2-TW)-^Ci.4l»(Py-TC-TN-SC) 

+0. 34*  (Rl+R2-TR)^0.26*CLI*0.072*LIL1+0. 26»NP)  +0. 5*C 


A 'a  aV  J  m.'i 


-.**  •  •  /•*  **• 


I»0.5»(-16.71+0.78*(PL1-TCLI-TNLHR1L1+R2L1-TRL1  +  DLI) 

&  -0.073*KL1>0. 14#L2Ll)+0.5*I 

♦  3 

SC»0.5«(-3.S3+0.72*(PC-TC)+0.07A#(PCL1-TCL1-SCL1)-0.028*SBL1) 

Sc  +0.5*SC 

*4 

PC»0.5*(-7.60+0.68#PY)+0.3#PC 

♦5 

D«0.5*(7.25+0. 10*(K+KL1) /2+0.044*(Q-W2) )+O.S#D 

#6 

W1»0.5#(-1.40+0.24#(Q-W2)+0.24*(QL1-W2L1)+0.29#T)+0.5*WI 

♦  7 

NW=0.5«( ( ( {Q-W2) +26.08-0. 16* (K+KLl ) /2. -2. 05*T ) / 2. 17-NE) /H+NG) 

Sc  +0.S*NM 

>8 

WR=0.5#(4. ll-0.74*(NL-NW-NE)+0.52*(PILl-PIL2)+0.54»T+WRLl) 

Sc  +0.5*WR 

♦  9 

Fl=0.5*(0.32+0.006*(M-TW-TC-TN-TR)+PI/PF+0.81»FILl)+0.5tFI 

♦  10 

R1»0. 5* (PR/PI) *(-0.36+0. 054* (W1+W2-TW+PY-TC-TN-SC)*PI /PR 
Sc  -0.007*(  (lillLl+W2Ll-TMLl+PLl-TCLl-TNLl-SCLl)*PILl/PRLl) 

Sc  +0,012*FR)+0.5*R1 

♦  11 

PR»0.5*(-131. 17+2.32*PI)+0.5«PR 

+♦♦**♦*♦♦«♦♦«*»«««**»###««###***«******•*«**************««*»»*»*•****** 

*  The  monetary  sector  is  onitted  (See  Chapter  IV) 

*12 

C  Ll=0.5*(0. 14*(M-TW-TC-TN-SC-TR)*76.  03*( (IL-2.0)**(-0.84) ) ) 

C  Sc  +0.5*Ll 
*13 

C  L2=0. 5* (-0. 34+0. 26*W 1-1.02*15-0. 26* (PI-PILl)+0. 6 1*L2L1) +0.5*12 

♦*♦#♦#♦*♦*♦♦♦♦♦♦♦♦♦#*#♦♦*♦♦♦#*♦**♦♦♦*#»♦«♦♦*♦♦♦♦♦*♦♦♦**♦♦♦♦♦♦♦♦*♦#♦*♦♦♦ 


*14 
*15 
♦  16 

♦  17 

♦  18 

♦  19 

♦  20 


IL=0.5*(2.58+0.44*ISL3+0.26*ISL5)+0.5*IL 

IS  =  0.5*(100*ISL1/ (100-11.  17+0. 67*LB) ) +0.5*  IS 

K  =  0. 5* ( I-D  +  KLl ) +0. 5*K 

SB=0.5* (SC+SBLl ) +0. 5*SB 

PYs0.5*(M-Wl-W2-Rl-R2)+0.5*PY 

M*0.5*(C+I+G-FI-TE-D)+0.5*t1 

PI=0.5*(  (H*NW*WR)  /  (lill+W2)  )+0.5*PI 


♦21 


Q»0. 5* (M+TE+D) +0. 5*Q 


C  SET  XI  s  NEW  VALUES 


XI (l)»C 
XI (2)>I 
XI (3)»SC 
XI <4)»PC 
X1(5)>D 
XI (6)>W1 
XI (7)«NW 
XI (a)»WR 
XI (9)«FI 
XI (10)*Rl 
XI  (n)=PR 
XI (12)=LI 
X1(13)=L2 
Xl(14)sIL 
XI (15)-IS 
X1(16)»K 
XI (17)=SB 
XI (ia)=PY 
XI (19)*M 
X1(20)=PI 
XI (21)=Q 


RETURN 

END 

**«4«*«****«««««****««tt****«««**«*«*«t«*«*tt***«***t #*#***«***# 

subroutine  update(X0,  XI,  P,  IS,  PI) 

*  Updates  values  lor  linking  -forecasts  together 

double  precision  X0(2n,  Xl(2l),  P(44),  IS(5),  PI(2) 


*  Update  lagged  endogenous  variables 


IS(5)=IS(4) 

IS(4)=IS(3) 

ISt3)=IS(2) 

1S(2)=IS(1) 

IS(l)=Xl(15) 

PI  CjsPI  (1) 
PH  I )  =  Xl  (20) 


♦  CLl 

P(18)*Xl  (1) 


♦  SCLl 


P(19)=X1 (3) 


J'VT’ 


♦  PCLl 
#DL1 
tMlLl 
«NRL1 
#FIL1 

♦  RlLl 
*PRLl 

♦  LlLl 

♦  L2L1 

♦  ISLl 

♦  ISL3 

♦  ISL5 

♦  KLl 

♦  3BL1 

♦  PLl 

♦  PILI 

♦  PIL2 
*QLi 


P(20)«X1 (4) 
P(21)»Xl (5) 
P(22)>X1 (6) 
P(23)>X1  (8) 
P(24)«X1 (9) 
P(25)=X1 (10) 
P(26)-X1(11) 
P(27)=X1 (12) 
P(2a)=Xl  (13) 
P(29)=IS(1) 
P(30)*IS(3) 
P(3l)=IS(5) 
P(32)=X1 (16) 
P(33)»X1 (17) 
P(34)sXl  (18) 
Pi35)*PI(l) 

P (36) »PI (2) 
P(37)*X1 (21) 


*  uoaati  laggtd  axagenaus  variables 

♦  TELl 


♦  ThILI 

♦  TCLl 

♦  'NLl 

♦  TRlI 

♦  w2lI 

♦  ftiLl 


P  i30)  »P  U  ) 
P i39) »P  12) 
P (40; *P (3) 
P ;41 ; =P (4) 
P(42;=P(5) 
P (43) =P (7) 
P  i44) =P (8) 
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*TII1E  TREND  (Suppressed  for  sone  runs) 
C  P(9)=P(9)+1 

*  Use  XI  as  new  starting  guess 

*  do  100  isi,21 

*  X0(i)-Xl(i) 

tlOO  continue 

return 

end 


Appendix  B.  Ea\iL^4LgQ£8  9i  duLtlBlitCi  tS4St  Squares  Cgetiicients 

iac  tLQSiC  l^stemSs. 


Chapter  II  showed  how  the  coefficients  of  a  first  order  response 
surface  equation  could  estimate  multipliers  for  a  linear  system.  Here 
least  squares  coefficient  will  be  shown  to  be  equivalent  to  the 
multipliers. 

Let 

Y  =  the  n  by  m  response  matrix  containing  n  observations  on  a 
endogenous  variables 

X  =  the  n  by  k  predetermined  variable  data  matrix  with  k 
predetermined  variables  and  n  observations. 

D  =  the  k  by  ffl  multiplier  matrix 

Then  the  linear  system  can  be  written 

i  -  XD 

Let  B  be  the  k  by  m  least  squares  coefficient  matrix.  B  is  defined  as 

B  =  (X'X)-‘X'Y 
=  (X  •  X) - ' (X  ■  X)  D 
=  ID 
B  =  0 

Therefore,  the  least  squares  coefficients  computed  for  a  first  order 
response  surface  fit  to  a  linear  model  are  identical  to  the  multipliers 


for  'he  same  oerioo. 


Appendix  C.  Decoding  Second  Order  BSSBSDSS  Coefficients 

After  estimating  response  surface  coefficients  in  terms  of  the 
coded  experimental  design,  the  response  surface  equation  must  be  re¬ 
expressed  in  terms  of  the  original  variables.  Decoding  the  coded 
response  surface  coefficients  for  a  second  order  response  surface  equa 
tion  IS  time  consuming,  tedious,  and  prone  to  errors.  This  appendix 
outlines  a  method  for  decoding  coded  coefficients  using  matrix  algebra 
which  simplifies  the  decoding  process.  If  a  computer  with  routines 
caoable  of  matrix  inversion  and  multiplication  is  available,  decoding 
can  be  made  much  easier. 

Chapter  II  gave  the  formula  for  translating  the  ith  original 
decodeo  variable  to  coded  form. 

X  c  »  ®  Xoi  —  (  3  .  • 

wner  9 

xc>  =  the  coded  x  value, 

y.a,  =  the  original,  decoded  x  value, 

=  the  maximum  factor  level,  and 
=  the  minimum  factor  level. 

Let 

X.  =  (X,m.«+’X,m,r.)/2 

^Xi  *  (X.mai.— X.minl/^ 

Eo  (3.2)  can  be  rewritten 

Xc.  =  5.1,.  ■  iC.l) 

i  X  . 

The  coded  second  order  response  surface  quation  is  a  auadratric 
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farin  Mhich  can  be  written 


y  = 


bc>  jXci  Xcj 


(C.2) 


where 


y  -  the  response  variable, 

xc>  and  xc,  =  the  ith  and  jth  coded  independent  variables, 
bcij  =  the  coeHicient  ai  the  product  of  the  ith  and  jth  coded 
var  tables, 

k  =  the  number  of  factors, 

;<co  =  1,  and 

bcoo  =  the  intercept  term. 

In  regression  program  outputs  bc>  ,  and  bo,  are  summed  because  X( 
xr,.  Consequently  bc.j  in  Eq  (C.2)  is  half  the  value  appearing  as  a 
regression  coefficient  in  a  regression  package  output. 

Bv  defining  appropriate  vectors  and  matricies,  Eq  (C.2)  can  be 
written  in  matrix  form.  If  k  is  the  number  of  factors  (independent 
variables)  in  the  response  surface  equation,  then  let  xo  be  the  k+l 
element  column  vector  whose  first  element  is  one  ano  the  remaining 
elements  are  the  k  decoded  independent  variables.  Similarly,  let  xr:  b 
the  k+l  element  column  vector  whose  first  element  is  one  and  the 
remaining  elements  are  tne  k  coded  independent  variables.  For  example 


Let  Be  be  the  matrix  whose  elements  are  bcu  and  Bo  be  the  matrix  whose 
elements  are  boij.  Then  Eq  (C.2)  can  be  rewritten 


y  -  xcBcXc 


Similarly, 


y  «  Xo’’BcXd 


It  follows  that 


Xc^'BcXc  =  Xo’^BdXo 


(C.3) 


Bo  contains  the  coefficients  of  the  decoded  independent  variables  which 
are  desired.  It  is  convenient  to  solve  for  Bo  in  terms  of  Be,  x,,  and 
Xi,  Let  A  be  the  matrix  which  transforms  xo  to  xc. 


AXo  *  Xc 

tnen  Eq  (C.3)  can  be  rewritten 


Xc’’BcXc  *  (AXd)^BdAXo 
=  XoMA^BcA)xo 


Then  it  follows  that 


Xo’^BoXo 


Bn  =  A’^BcA 


It  will  be  demonstrated  but  not  proven  that 

A  =  C  -‘  d-xu’’)  I.C.4) 

wner  e 

A  *  the  transformation  matrix, 

C  =  a  x-t’i  bv  kti  diagonal  matrix  whose  ith  diagonal  element  is  x, 
(define  xo  =  1), 


I  =  a  k+l  by  k*l  identity  matrix 


H  -  a  k-t'l  slenint  colu«in  vector  Mhose  ^irst  elenent  is  zero  and  the 
remaining  elements  are  x. ,  and 

u  -  a  k'M  element  column  vector  whose  first  element  is  one  and  the 
remaining  elements  arc  zero. 

An  example  demonstrates  the  validity  of  Eq  (C.4).  If  k32  then 


1 

0 

0 

’0“ 

c  = 

0 

^X  1 

0 

X  a 

X  1 

0 

0 

AX  I 

-Xa 

_ 1 

[1  0  0] 

0 

0 

0" 

X  1 

X  t 

0 

0 

I  z 

XT' 

0 

0 

I 


0 


0" 

0 


Lo 

0 

U 

0 

Oj 

L-xj  0 

ij 

1  /AXo 

0 

0 

1  /AX  1 

0 

0 

0 

1  /A  X 

'1 

0 

0 

'1 

0 

o' 

-  xu^ 

)  = 

0 

1  /Ax  1 

0 

• 

-X  , 

1 

0 

0 

0  1/axt-_ 

-XT' 

0 

1 

U3 


by  Eq  (C. 1 ) .  Then 

Bd  =  CC  ‘(I-xaMJ^BcCC  MI-KuM]  (C.5) 

which  IS  the  desired  result.  Eq  (C.5)  is  used  to  decode  coded  resoonse 
surtace  coeHicients  in  this  thesis.  It  saves  time  and  eHort. 
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Appendix  D.  Decoding  the  |xger i.Bgntal^  Design 

Coded  experimental  designs  are  used  to  collect  data  ^or  Htting 
response  surfaces.  In  order  to  determine  what  factor  levels  to  run  an 
experiment  or  model  at,  the  experimental  design  must  be  decoded  from 
Is,  O's,  and  -Is  to  actual  factor  levels.  Below  is  a  FORTRAN  program 
which  accomplishes  this  task.  The  program  reads  the  coded  experimental 
design  and  factor  levels  from  user  specified  files,  and  writes  the  coded 
experimental  design,  factor  levels,  and  decoded  experimental  design  to  a 
user  specified  file.  At  the  end  of  this  section  is  sample  output.  The 
output  contains  the  Box  and  Behnken  three  level  five  factor  coded 
experimental  design  used  so  extensively  in  this  research  effort.  The 
output  also  contains  the  factor  levels  which  are  the  design  variable 
nigh  and  low  sample  values.  After  the  factor  levels  is  the  decoded 
experimental  design. 

PROGRAM  DECODE 

*  Translates  a  coded  two  level  or  three  level  experimental 

*  design  matrix  to  a  design  matrix  with  actual  factor  levels. 

*  The  factor  level  file  must  have  MAXIMUM  values  followed  by 

*  MINIMUM  values. IMPORTANT '  The  data  in  each  of  the  input 

*  design  files  must  have  the  array  dimensions  (m  n)  as  the 
»  first  line. 


♦DECLARATIONS* 

integer  mD,  nD,  mF,  nF,  maxm.  maxn, 
run  1 1 1 ,  run  1 1 2 ,  wuni t 
parameter imaxm=200,maxn=50) 
double  precision  D t maxm , maxn )  ,  F(maxm,maxn) 
double  orecision  Dc (maxm, maxn) 
character*24  expdes.  faciev,  outfil 
run  1 1 1  =  I 
runit2=2 


'  ^  if  •  •‘  V  ,«.^ 


Munlt»3 


♦USER  INPUT# 


call  gatdat (axpclcs,  faclev,  outfill 
open  (runitl,  Hie  -  expdes,  status 
reMind  (runitl) 

open  (runit2,  Hie  *  faclev,  status 
rewind  (runit2) 

open  (wunit,  Hie  -  outfil,  status  » 


=  'old'  ) 


=  'old ' ) 


♦END  OF  USER  INPUT# 

call  nread (runitl ,  D,  eO,  nD,  Aaxn,  naxn) 

Mr i te (Muni t , 1 ) 

I  formate  EXPERIMENTAL  DESIGN  MATRIX') 

call  fflwritelMunit,  D,  mD,  nD,  maxm,  maxn) 

call  mread (runi t2,  F,  mF,  nF,  maxm,  maxn) 

write(wunit,#) '  FACTOR  LEVEL  MATRIX' 

call  mwr i te (wuni t ,  F,  mF,  nF,  maxm,  maxn) 

call  switch(D,  mD,  nD,  F,  maxm,  maxn,  Oc) 

write(Munit,#) '  DECODED  MATRIX' 

call  mMrite(Munit,  Dc,  mD,  nD,  maxm,  maxn) 

write!#,#)  '  All  done!  ' 

stop 

end 

»*##*#«*####*«#«*«*SU6R0UTINES##########*«##################«**#######* 


subroutine  getdat (expdes ,  faclev,  outfil) 


DESCRIPTION:  Requests  experimental  design  file  name, 
factor  level  file  name,  and  output  file 
name  from  user  at  terminal. 

INPUT:  Input  and  output  file  names  supplied  by  the  user. 
OUTPUT:  Variable  values  for  expdes,  faclev,  and  outfil. 

cnaracter#24  expdes,  faclev,  outfil 

write(#,»)  Enter  coded  experimental  design  file  name, 
read  (#,10)  expaes 
format (a24) 

write!#, ♦)'  Enter  factor  level  file  name. 
read(#,10)  faclev 

write(#,#)'  Enter  output  file  name. 
read(#,10)  outfil 
format (II) 
return 


f*tf *«#*#*##*«*«*«#««««*«••*««*«••««««*«•*«««**««*«**«§**#*#«•#**«**#«* 


subroutine  ereadCrunit,  A,  a,  n,  aaxa,  aaxn) 

DESCRIPTION:  Reads  diaensions  and  values  ior  aatrix 
A. 

INPUT:  Values  'for  diaensions  and  values  for  eleaents 
of  aatrix  A  froa  a  file.  Max  sise  froa  calling 
progr aa. 

OUTPUT:  Matrix  A,  diaensions  a,  n  to  calling  orograa. 

integer  runit,  a,  n,  aaxa,  aaxn,  i,  j 
double  precision  A(aaxa,  aaxn) 

read  (runit,  *)  a,  n,  ((A(i,  j),  j  «  1,  n),  i  »  I,  a) 
return 
end 


i^************************************* **«**#**•**««**«*«*«#«««**««***«* 


subroutine  awriteiwunit.  A,  a,  n,  aaxa,  aaxn) 

DESCRIPTION:  Writes  a,  n,  and  then  the  aatrix  A  in 
rows. 

INPUT;  Array  A  with  diaensions  a,  n  and  aax  size  aaxa,  aaxn 
froa  calling  prograa. 

OUTPUT:  Values  for  a,  n,  and  eleaents  of  A  to 
standard  output. 

integer  a,  n,  aaxa,  aaxn,  i,  j,  wunit 
double  precision  A(aaxa,  aaxn) 

writeiwunit,  1000)  a,  n 
f oraat (  '  ,  12, '  BY  ,  12) 

do  1  ISO  1  =  1 ,  a 

write(wunit,  1100)  (A(i,  j),  j  =  1,  n) 
foraat ( '  '  ,5E14.5) 
continue 
wr 1 te (wuni t , ♦) 
retur n 


"’TJ 


subroutine  switchiDc 


7 


V 


♦ 

* 

♦ 

» 

» 


’  I  »  •-  .-«  k  •  ».-•  k-*  k»  i.'Si'V  C~VV> 


DESCRIPTION:  Creates  a  new  array  containing  factor 
place  of  coded  values  fron  the  design 
INPUT:  Coded  design  natrix,  Dc,  with  dieensions  nD 
naximua  array  dieensions  naxn,  and  naxn 
OUTPUT:  Decoded  design  matrix,  Dd 


levels  in 
natrix. 
by  nD, 


integer  mD,  nD,  oiaxoi,  maxn,  i,  j 

double  precision  Dc (aaxffl,aaxn) ,  F(maxn,naxn)  ,  Dd(naxn, 
*  fflaxn) 


I 


do  40  i  -  1,  inD 

do  30  j  =  1,  nD 

Dd (i , j )  =  (F(l , j)+F (2, j) ) /2  +  Dc (1 , j ) ♦ (F  ( 1 , j ) -F  (2,  j )  )  /  Z 
30  continue 

40  continue 

return 
end 


♦  End  of  FORTRAN  code. 

4****«**«***««*******«**««*«««««*****«**«*****««****t*«**«*«««****«*«** 


Sample  Output  Of  DECODE  (Box  and  Behnken  Three  Level 
Five  Factor  Oetign. 


EXPERIMENTAL  DESIGN  MATRIX 
41  BY  3 

0.  e+00  0.  e+00  0.  e+00  0.  e+00  0.  e+00 

O.lOOOOe+Ol  O.lOOOOe+Ol  0.  e+00  0.  e+00  0.  e+00 

O.lOOOOe+01  -O.lOOOOe+Ol  0.  e+00  0.  e+00  0.  e+00 

-O.lOOOOe+01  O.lOOOOe+Ol  0.  e+00  0.  e+00  0.  e+00 

-O.lOOOOe+01  -O.lOOOOe+01  0.  e+00  0.  e+00  0.  e+00 

0.  e+00  0.  e+00  O.lOOOOe+Ol  O.lOOOOt+01  0.  e+00 

0.  e+OO  0.  e+OO  O.lOOOOe+01  -O.lOOOOe+Ol  0.  e+00 

0.  e+00  0.  e+OO  -O.lOOOOe+Ol  O.lOOOOe+Ol  0.  e+00 

0.  e+00  0.  e+00  -O.lOOOOe+Ol  -O.lOOOOe+01  0.  e+00 

0.  e+00  O.lOOOOe+01  0.  e+00  0.  e+00  O.lOOOOe+01 

0.  e+00  O.lOOOOe+01  0.  e+00  0.  e+00  -O.lOOOOe+01 

0.  e+00  -O.lOOOOe+Ol  0.  e+OO  0.  e+00  O.lOOOOe+01 

0.  e+00  -O.lOOOOe+01  0.  e+00  0.  e+00  -O.lOOOOe+Ol 

O.lOOOOe+Ol  0.  e+00  O.lOOOOe+Ol  0.  e+OO  0.  e+00 

O.lOOOOe+Ol  0.  e+00  -O.lOOOOe+01  0.  e+00  0.  e+00 

-O.lOOOOe+01  0.  e+00  O.lOOOOe+Ol  0.  e+00  0.  e+00 

-O.lOOOOe+Ol  0.  e+00  -O.lOOOOe+Ol  0.  e+00  0.  e+00 

0.  e+OO  0.  e+00  0.  e+00  O.lOOOOe+Ol  O.lOOOOe+Ol 

0.  e+00  0.  e+OO  0.  e+00  O.lOOOOe+Ol  -O.lOOOOe+Ol 

0.  e+OO  0.  e+00  0.  e+OO  -O.lOOOOe+Ol  O.lOOOOe+Ol 

0.  e+00  0.  e+00  0.  e+00  -O.lOOOOe+Ol  -O.lOOOOe+01 

0.  e+00  O.lOOOOe+01  O.lOOOOe+Ol  0.  e+00  0.  e+00 

0.  e+00  O.lOOOOe+01  -O.lOOOOe+01  0.  e+00  0.  e+00 

0.  e+00  -O.lOOOOe+Ol  O.lOOOOe+Ol  0.  e+00  0.  e+00 

0.  e+00  -O.lOOOOe+Ol  -O.lOOOOe+01  0.  e+00  0.  e+00 

0.  lOOOOe+01  0.  e  +  00  0.  e+00  O.lOOOOe+Ol  0.  e+00 

O.lOOOOe+01  0.  e+00  0.  e+00  -O.lOOOOe+01  0.  e+00 

-O.lOOOOe+Ol  0.  e+00  0.  e+00  O.lOOOOe+Ol  0.  e+00 

-O.lOOOOe+01  0.  e+00  0.  e+00  -O.lOOOOe+Ol  0.  e+00 

0.  e+00  0.  e+00  O.lOOOOe+Ol  0.  e+00  O.lOOOOe+Ol 

0.  e+00  0.  e+00  O.lOOOOe+Ol  0.  e+00  -O.lOOOOe+01 

0.  e+00  0.  e+OO  -O.lOOOOe+Ol  0.  e+OO  O.lOOOOe+Ol 

0.  e+00  0.  e+OO  -O.lOOOOe+Ol  0.  e+00  -O.lOOOOe+01 

O.lOOOOe+01  0.  e+00  0.  e+00  0.  e+00  0. lOOOOe+01 

O.lOOOOe  +  01  0.  e  +  00  0.  e+00  0.  e  +  00  -0.  lOOOOe  +  01 

-O.lOOOOe+01  0.  e+00  0.  e+00  0.  e+00  O.lOOOOe+01 

-O.lOOOOe+Ol  0.  e+00  0.  e+00  0.  e+00  -O.lOOOOe+01 

0,  e+00  O.lOOOOe+Ol  0,  e+00  O.lOOOOe+Ol  0.  e+00 

0.  e+00  O.lOOOOe+Ol  0.  e+00  -0. lOOOOe+Ol  0.  e+00 

0.  e+00  -O.lOOOOe+Ol  0.  e+00  O.lOOOOe+Ol  0.  e+00 

■j.  e+00  -O.lOOOOe+01  0.  e+Oo  -O.lOOOOe  +  01  0.  e  +  00 

FACTOR  level  MATRIX 
-  BY  5 

0.11a30e  +  02  0.l3140e  +  02  0.6l700e+02  0.21820e  +  02  0.  10400e  +  o2 

0.5a300e  +  0l  0.7  1400e  +  0l  0.3"’700e*02  0.l3820e  +  02  0.  10400e  +  02 


V’  ‘C* 

'j 

'w 

DECODED  MATRIX 

r. 

1 

41  BY  5 

0.36300e-*'01 

0.101408+02 

0.497008+02 

0.  170208+02 

0.  104008  +  02 

0.116308+02 

0.131408+02 

0.497008+02 

0.  170208+02 

0. 104008+02 

0.116308+02 

0.714008+01 

0.497008+02 

0. 170208+02 

0.  104008  +  02 

0.563008+01 

0.  131408+02 

0.497008+02 

0. 170208+02 

0.  104008+02 

0.563008+01 

0.714008+01 

0.497008+02 

0.  170208+02 

0.  104008+02 

1 

0.063008+01 

0. 101408+02 

0.617008+02 

0.210208+02 

0. 104008  +  02 

IK. 

0.863008+01 

0. 101408+02 

0.617008+02 

0.  130208+02 

0.  104008+02 

» , 

0.063008+01 

0. 101408+02 

0.377008+02 

0.210208+02 

0.104008+02 

0.363008+01 

0.  101408  +  02 

0.377008+02 

0.  130208+02 

0.  104008  +  02 

-V 

1.;. 

0.063008+01 

0. 131408+02 

0.497008+02 

0.  170208  +  02 

0.  104008  +  02 

'•j 

0.063008+01 

0. 131408+02 

0.497008+02 

0.  170208  +  02 

0.  104008  +  02 

i 

0.863008+01 

0.714008+01 

0.497008+02 

0.  170208+02 

0.  104008  +  02 

p 

0.063008+01 

0.714008+01 

0.497008+02 

0.  170208  +  02 

0.  104008  +  02 

-• 

0.  116308  +  02 

0.  101408  +  02 

0.617008+02 

0.  170208+02 

0.  104008+02 

0. 116308+02 

0. 101408+02 

0.377008+02 

0.  170208  +  02 

0.  104008  +  02 

0.563008+01 

0. 101408+02 

0.617008+02 

0.  170208+02 

0.  104008+02 

0.563008+01 

0. 101408+02 

0.377008+02 

0.  170208  +  02 

0.  104008  +  02 

•  ♦,*  « 

1 

0.363008+01 

0.  101408  +  02 

0.497008+02 

0.210208+02 

0. 104008+02 

0.063008+01 

0.  101408+02 

0.497008+02 

0.210208+02 

0.  104008+02 

•*. 

0.063008+01 

0.  101408  +  02 

0.497008+02 

0.  130208+02 

0. 104008+02 

•  »> . 

0.063008+01 

0.  101408+02 

0.497008+02 

0.  130208+02 

0.  104008+02 

•'V: 

0.363008+01 

0.131408+02 

0.617008+02 

0.170208+02 

0. 104008+02 

• 

• 

0. 063008+01 

0.  131408+02 

0.377008+02 

0.  170208  +  02 

0.  104008  +  02 

■ 

0.363008+01 

0.714008+01 

0.617008+02 

0.  170208+02 

0. 104008+02 

WD  ' 

1 

0.863008+01 

0. 714008  +  01 

0.377008+02 

0.  170208+02 

0.  104008+02 

0. 116308+02 

0.  l0140e+02 

0.497008+02 

0.210208+02 

0. 104008+02 

0. 116308+02 

0.  101408+02 

0.497008+02 

0.  130208+02 

0.  104008  +  02 

0.563008+01 

0.  10140e+02 

0.497008+02 

0.210208+02 

0. 104008+02 

0.563008+01 

0.  101408  +  02 

0.497008+02 

0.  130208+02 

0. 104008+02 

h 

0.863008+01 

0. 101408+02 

0.617008+02 

0.  170208  +  02 

0. 104008+02 

& 

1 

0. 363008+01 

0.  101408+02 

0.617008+02 

0.  170208  +  02 

0. 104008+02 

0.063008+01 

0.  101408  +  02 

0.377008+02 

0.170208+02 

0. 104008+02 

V 

0. 063008+01 

0.  101408  +  02 

0.377008+02 

0.  170208  +  02 

0. 104008+02 

0. 1 16308+02 

0. l0l40e+02 

0.497008+02 

0.  170208  +  02 

0. 104008+02 

0. 116308+02 

0.  I0140e  +  02 

0.497008+02 

0.  170208  +  02 

0. 104008+02 

»-• 

0. 563008+01 

0. 101408+02 

0. 497008+02 

0.  170208  +  02 

0. 104008+02 

i 

0.563008+01 

0.  101408  +  02 

0.497008+02 

0.  170208  +  02 

0. 104008+02 

►  ■ 

0. 363008+01 

0. 131408+02 

0. 497008  +  02 

0.210208+02 

0. 104008+02 

• . 

0. 063008+01 

0. l3l40e+02 

0.  497008+02 

0.  136208  +  02 

0. 104008+02 

■ 

0. 363008+01 

0. 714008+01 

0.497008+02 

210208  +  02 

0. 104008+02 

363008+01 

0. 714008+01 

0.497008+02 

0.  136208  +  02 

0. 104008+02 

i 

>* 

.*'• 

/• 

i 

K\ 

1 

f: 

i  ■* . 

L’- 

•  •*•  •*•  •  -  **.  •  .  •  .  • 

120 

..  1.' 

• 

C*  J 

•  w  +  • 

*'v‘ 

t  • 

‘/-'J 

Appendix  E.  Response  Surface  Coefficients 


Tables  E. la-k  below  contain  both  coded  and  decoded  coefficient 
aatricies  for  the  second  degree  polynoaial  response  surface  equations 
described  in  Chapter  V.  The  Matrix  algebra  Method  in  Appendix  0 
was  used  to  decode  the  coded  coefficients.  Consequently,  off  diagonal 
elements  in  Tables  E.la-k  are  half  of  the  value  normally  given  as  a 
coefficient  for  a  polynomial  (See  Appendix  D.).  The  column  and  row 
marked  with  a  1  contain  first  degree  term  coefficients.  The  upper  right 
hand  corner  element  in  the  tables  is  the  intercept  term. 

Table  E.la.  Coded  Second  Order  Polynomial  Response  Surface  Coefficients 
for  Number  of  Workers  Employed  in  Period  Zero 


(R  = 

=  1.0000 

,  Adjusted 

R=:  8 

1.0000) 

1 

Tw 

___Ic _ 

_ G__ 

_ W? _ 

_ No _ 

1 

44. 1080 

-0.9029 

-0.1937 

4.6116 

-0.8731 

1.6890 

T« 

-0,9029 

0.0000 

0.0000 

0.0039 

-0.0033 

0.0000 

Tc 

-0. 1937 

0.0000 

0.0000 

0.0026 

0.0000 

0.0000 

G 

4.6116 

0.0039 

0.0026  - 

0. 0088 

0.0069 

-0.0042 

W, 

-0.8731 

-0. 0033 

0.0000 

0. 0069 

0.0000 

0.0000 

Nr, 


1.6390  0.0000  0.0000  -0.0042 


0.0000  0.0000 


Table  E.lb.  Decoded  Second  Order  Polynoeial  Response  Surface  Coefficients 
for  Nunber  of  Workers  Employed  in  Period  Zero 


•C'* 


(R  = 

»  1.0000, 

Adjusted  R^  -  1 

.0000) 

1 

T« 

Tc 

6 

No 

1 

23.2089 

-0.1686  - 

■0.0345 

0.3059 

-0. 1366 

0.4989 

Tu 

-0.1606 

0.0000 

0.0000 

0.0000 

-0.0001 

0.0000 

Tc 

-0.0345 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

G 

0.3059 

0.0000 

0.0000 

0.0000 

0.0001 

-0.0001 

W, 

-0. 1366 

-0.0001 

0.0000 

0.0001 

0.0000 

0.0000 

Ng 

0.4989 

0.0000 

0.0000 

-0.0001 

0.0000 

0.0000 

E.  ic. 

Coded  Second  Order 
for  Price 

Polynomial  Response  Surface 
Index  in  Period  Zero 

Coef f  ic 

(R  = 

=  1.0000, 

1  Adjusted  R=  *  1 

.0000) 

1 

Tw 

Tc 

G 

Ws 

Ng 

1 

134.5839 

-2.1629  • 

■0.4619 

11.0416 

-9.4508 

6. 1408 

Tw 

-2. 1629 

0.0000 

0.0000 

0. 1592 

0. 0776 

0.0000 

Tc 

-0.4619 

0.0000 

0.0000 

0.0398 

0.0000 

0.0000 

G 

11.0416 

0. 1592 

0.0398 

-0.7819 

-0.4222  - 

•0.  1191 

-9.4508 

0.0776 

0.0000 

-0.4222 

1,8014  - 

0.6453 

Ng 

a. 1408 

0. 0000 

0. 0000 

-0. 1191 

-0. 6453 

0. 1668 

Table  E.ld.  Decoded  Second  Order  Polynomial  Response  Surface  Coefficients 

for  Price  Index  in  Period  Zero 


(R  = 

1 

-  l.OOOO,  Adjusted  R^  -  1. 

Tw  Tc  G 

0000) 

Ng 

1 

100.3344 

-0.4786 

-0.0927 

0.8‘’42 

-1.5949 

2. 

0708 

Tw 

-0.4786 

0.0000 

0.0000 

0.0020 

0.0022 

0. 

0000 

Tc 

-0.0927 

0.0000 

0.0000 

0.0005 

0.0000 

0. 

0000 

G 

0.3742 

0.0020 

0.0005 

-0.0034 

-0. 0043 

-0. 

0023 

-1.5949 

0.0022 

0.0000 

-0.0043 

0.0431 

-0. 

0294 

Nr, 

2.0700 

0.0000 

0.0000 

-0.0023 

-0. 0294 

0. 

0144 

Table  E.le.  Coded  Second  Order  Polvnomial  Response  Surface  Coefficients 
for  Sross  National  Product  in  Period  Zero 


(R- 

1 

=  1.0000 

T« 

,  Adjusted  R=  =  1. 

Tc  G 

0000) 

w? 

Nr, 

1 

113.3169 

-2.0454 

-0.4398 

10.4517 

1.2531 

-0.0251 

1  M 

-2.0454 

0.0000 

0.0000 

0.0088 

-0. 0073 

.  0044 

T  C 

-0. 4398 

0. 0000 

-0.0008 

0.0051 

0.0000 

0. 0000 

G 

10.4517 

0.0088 

0.0051 

-0.0186 

0.0151 

-0.0094 

1^7 

1.2531 

-0. 0073 

0.0000 

0.0151 

0. 0000 

0. 0000 

Nr. 

-0. 0251 

0 . 0044 

0.0000 

-0.0094 

0.0000 

0 . 0000 

Table  E.lf.  Decoded  Second  Order  Polynoaial  Response  Surface  Coetficients 
for  Gross  National  Product  in  Period  Zero 


(R= 

=  1.0000 

,  Adjusted  R=  »  1 

.0000) 

1 

Tw 

T  c 

6 

Ws 

Ng 

70.6007 

-0.3836 

-0.0780 

0.6933 

0. 1904 

-0.0033 

-0.3836 

0.0000 

0.0000 

0.0001 

-0.0002 

0.0002 

-0.0780 

0.0000 

0.0000 

0.0001 

0. 0000 

0.0000 

0.6933 

0.0001 

0.0001 

-0.0001 

0.0002 

-0.0002 

0. 1904 

-0.0002 

0.0000 

0.0002 

0.0000 

0.0000 

-0.0033 

0.0002 

0.0000 

-0.0002 

0.0000 

0.0000 

Coded  Second  Order  Polynonial  Response  Surface 
for  Number  of  Workers  Employed  in  Period  Five 

Coef  f  ii 

(R  = 

=  1.0000 

,  Adjusted  R-  »  0 

.9999) 

1 

Tc 

G 

w. 

Ng 

49.4412 

-3.4778 

-3.5475 

13.4030 

-1.2750 

1.4146 

-3.4778 

0.0000 

-0.0939 

0.4135 

-0. 1006 

0. 0000 

-3.5475 

-0.0939 

0.0000 

0,4108 

-0. 0934 

0.0750 

13.4030 

0.4135 

0.4180 

-1.2752 

0.3521  - 

0.2912 

-1.2750 

-0. 1006 

-0. 0934 

0.3521 

0.0000 

0.0000 

1.4164  0.0000  0.0750  -0.2912 


0.0000 


0. 0000 


..-w  ^  .  V  H ""  ^  ^  Ji 


Taole  E.lh.  Decoded  Second  Order  Polynomial  Response  Surface  Coefficients 
for  Number  of  Workers  Employed  in  Period  Five 

(R=  =  1.0000,  Adjusted  R=  =  0.9999) 


1 

T„ 

Tc 

G 

Ws 

Ns 

1 

10.3560 

-0.7353 

-0.7372 

0.9930 

-0.2682 

0.5433 

Tu, 

-0.7353 

0.0000 

-0.0030 

0.0051 

-0.0029 

0.0000 

Tr. 

-0.7350 

-0.0030 

0.0000 

0.0047 

-0.0025 

0.0030 

G 

0.9925 

0.0051 

0.0048  - 

0.0056 

0.0036 

-0.0057 

Wo, 

-0.2682 

-0.0029 

-0.0025 

0.0036 

0.0000 

0.0000 

Nr. 

0.5439 

0. 0000 

0.0030  - 

0.0057 

0.0000 

0.0000 

Table  E.  li . 

Coded  Second  Order 
for  Price 

Polynomial  Resp 
Index  in  Period 

onse  Surf 
Five 

ace  Coeff 

(R 

2  =  0.9990 

,  Adjusted 

R==  = 

0.9997) 

1 

Tm 

Tc 

S 

W, 

Ns 

1 

185.6951 

-14.8952 

-12.8819 

59.6286 

-18.4974 

15.2612 

T« 

-14.8952 

0. 0000 

0. 0000 

0.0000 

0.0000 

0.0000 

Tr 

tJ 

GO 

GO 

0.0000 

0.0000 

0. 8976 

0. 0000 

0.0000 

G 

59.6286 

0.0000 

0.8976 

0.0000 

-2.5389 

1.9044 

Wt 

-18.4974 

0.0000 

0.0000 

-2.5389 

4.0012 

-2.5593 

No 

15.0212 

0.0000 

0. 0000 

1.9044 

-2.5593 

2.3429 
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Table  E.lj.  Decoded  Second  Order  Polynoeial  Response  Sur'face  Coe^'ficient 

for  Price  Index  in  Period  Five 


(R  = 

=  0.9998 

,  Adjusted  R^  »  0. 

9997) 

1 

Tw 

Tc 

G 

W2 

Ns 

1 

23. 1609 

-2.7764 

-2.5197 

3.8950 

-2.3456 

3.2888 

Tw 

-2.7764 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Tc 

-2.5197 

0.0000 

0.0000 

0.0104 

0.0000 

0.0000 

G 

3.8950 

0.0000 

0.0104 

0.0000 

-0.0260 

0.0371 

Us> 

-2.3456 

0.0000 

0.0000 

-0.0260 

0.0957 

-0.1164 

Ng 

3.3947 

0.0000 

0.0000 

0.0371 

-0.  1164 

0.2027 

E.  Ik. 

Coded  Second  Order  Polynomial  Response  Surface  Coeffic 
for  Gross  National  Product  in  Period  Five 

(R= 

=  1.0000 

,  Adjusted  R=  =  1. 

0000) 

1 

T„ 

Tc 

6 

Nr. 

1 

141.3628 

-9.2409  - 

10.6323 

36.0425 

-0.3852 

0.7046 

T„ 

-9.2408 

0.0000 

-0.2276 

0.9903 

-0.2246 

0.  1707 

Tc 

-10.6323 

-0.2276 

0.0000 

0.9796 

-0.2342 

0. 1769 

G 

36.0425 

0.9903 

0.9796 

-3.0659 

0.8543 

-0.7019 

M:; 

-0.3852 

-0.2446 

-0.2342 

0.8543 

0.0000 

0.  1568 

No 

0.7046 

0. 1707 

0.1769 

-0.7019 

0.  1568 

0.0000 

Table  E.ll.  Decoded  Second  Order  Polynoeial  Reiponse  Sur'facc  Coe'f'f icients 
for  Gross  National  Product  in  Period  Five 


(R 

=  =  1.0000 

,  Adjusted  R^  *  1. 

0000) 

1 

Tw 

Tc 

6 

Us 

No 

1 

43.4207 

-1.9947 

-2.1269 

2.6404 

-0.2824 

0.4108 

Tw 

-2.0001 

0.0000 

-0.0074 

0.0122 

-0.0065 

0.0094 

Tc 

-2. 1269 

-0.0074 

0.0000 

0.0113 

-0.0063 

0.0091 

G 

2.6404 

0.0122 

0.0113 

-0.0134 

0.0088 

-0.0137 

Us 

-0.2805 

-0.0071 

-0.0063 

0.0088 

0.0000 

0.0071 

Ng 

0.4108 

0.0094 

0.0091 

-0.0137 

0.0071 

0.0000 

Tables 

E.2a-e  contain  decoded  coefficients 

for  first 

order 

ynomial 

response 

surfaces. 

They  may 

be  compared  directly  to 

multipliers  computed  by  Goldberger  (Goldberger,  1959). 

Table  E.2a.  First  Order  Response  Surface  Coefficients  for  a  Unit 

Increase  in  Tw. 

Period 


0 

1 

2 

3 

4 

5 

Nw 

-0.3367 

-0.7538 

-0. 1081 

-1.2729 

-1.3347 

1.2937 

P 

-0.8293 

-1.9435 

-3.0677 

-4. 1201 

.  0475 

-5.3117 

Q 

-0.7631 

-0. 1738 

-2.5663 

-3. 12342 

-3.4011 

-3.4391 

Table  E.2b.  First  Order  Response  Surface  Ooe'ff icients  ^or  a  Unit 

Increase  in  Tc. 

Period 

0 _ 1 _ 2 _ 3 _ 4 _ 5 

Nw  -0.6769  -0.5960  -0.1038  -1.2756  -1.3243  -1.2333 

p  -0.1671  -1.4552  -2.6269  -3.5301  -4.1661  -4.5498 

Q  -0.1532  -0.1418  -2.5847  -3.3451  -3.6914  -3.7003 

Table  E.2c.  First  Order  Response  Surface  Coefficients  for  a  Unit 

Increase  in  6. 

Period 

_ 0 _ 1 _ 2 _ 3 _ 4 _ 5 

Nu  0.6110  0.1214  0.1621  1.8268  1.8623  1.7715 

p  1.5050  3.1437  4.6727  6.0734  7.3022  8.3124 

Q  0.13840  0.2804  3.0707  4.5208  4.47952  4.7662 

Table  E.2d.  First  Order  Response  Surface  Coefficients  for  a  Unit 

Increase  in  Wa. 

Period 

0 _ 1 _ 2 _ 3 _ 4 _ 5 

Nw  -0.2705  -0.4465  -0.5115  -0.5069  -0.4599  -0.3919 

p  -2.9658  -4.5079  -5.3932  -5.9814  -6.4147  -6.7606 

Q  +0.3867  -0.3642  -0.2389  -0.2863  -0.2314  -0.1149 


Apendix  F.  Optimization  Problem  Solution 


Shown  below  is  the  output  file  of  the  linear  programming  package 
for  the  optimization  problem  formulated  in  Chapter  VI.  The  output 
includes  a  problem  specification,  optimal  basic  variable  values, 
shadow  prices,  and  the  objective  function  value. 


Problem  Specified  for  Solution 
hax i mi ze 


TW 

TC  G 

M2 

X 

1  X 

2  X  3 

X 

4  8 

19 

Qbj 

- 

3.11 

-3.33 

4.50 

-0.29 

54.40 

Constraint 

1  - 

max  TM 

type 

is  le 

1.00 

0. 

0. 

0. 

11.63 

Constraint 

2  - 

max  TC 

type 

is  le 

0. 

1.00 

0. 

0. 

13.  14 

Constraint 

3  - 

max  G 

type 

is  le 

0. 

0. 

1.00 

0. 

61.70 

Constraint 

4  - 

max  M2 

type 

is  le 

0. 

0. 

0. 

1 . 00 

21.82 

Constraint 

5  - 

inf  latn 

type 

is  le 

- 

4.61 

-4.  12 

6.68 

-5.66 

96.79 

Constraint 

6  - 

unemp lay 

type 

IS  eg 

- 

1.25 

-1.26 

1.80 

-0.50 

45.  77 

Constraint 

7  - 

ram  TM 

type 

is  gt 

l.OO 

0. 

0. 

0. 

5.63 

Constraint 

8  - 

min  TC 

type 

is  gt 

0. 

1.00 

0. 

0. 

7.  14 

Constraint 

9  - 

min  G 

type 

IS  gt 

0. 

0. 

1 . 00 

0. 

37.70 

Constraint 

10  - 

min  W2 

type 

IS  gt 

0. 

0. 

0. 

1.00 

13.82 

Activity  variables  I  through  4 

Slack  variables  (S)  5  through  9 

Surplus  variables  (Pi  10  through  13 
Artificial  variables  (A)  14  through  18 


Answers 


Basic  Variables 

Value 

X  1 

TW 

a 

11.6300 

X  2 

TC 

= 

7.1400 

X  3 

G 

s 

44.5489 

X  4 

M2 

s 

21.8200 

S  6 

max  TC 

3 

6.0000 

S  7 

max  6 

S 

17. 1511 

S  9 

inf latn 

a 

5.4738 

SIO 

unemploy 

3 

6.0000 

S12 

min  TC 

3 

6.8489 

S13 

mi n  6 

S 

8.0000 

Increase  in  Qbj.  Function  ior  unit  increase  in 
right  hand  side  of  constraints 
Shadow  Prices  Value 


Y  1  :  max  TU  -  0.0246 

Y  4  :  max  W2  =  0.9688 

Y  7  :  am  TW  =  0.  1925 

no  :  min  W2  =  2.4985 


The  value  o^  the  objective  function  is: 


185.8114 
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